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von Friedrich Möglich 


1. Beugungserscheinungen an Körpern ii 
von ellipsoidischer Gestalt; 

(ye 

Es ist verschiedentlich versucht worden, eine Theorie der 
Beugung an Ellipsoiden zu entwerfen. Schon Lord Rayleigh?) 
hat, unter der Voraussetzung, daß die Ellipsoide klein sind 
im Vergleich zur Wellenlänge der einfallenden Strahlung, 
Formeln angegeben, die die Störung eines elektromagnetischen 
Zustandes durch ein der erwähnten Bedingung genügendes 
Ellipsoid auszudrücken. R. Gans?) hat dann diese Methode 
benutzt, um die Mieschen Untersuchungen, über die Optik 
trüber Medien zu vervollkommnen. Abraham‘) und Ehren- 
haft?) beschäftigen sich mit der Untersuchung der elektrischen 
Schwingungen eines Rotationsellipsoides, eine Aufgabe, die 
wesentlich darauf hinausläuft, die Wellengleichung so zu inte- 
grieren, daß die elektromagnetischen Grenzbedingungen an der 
Oberfläche des die Schwingungen erregenden Ellipsoides er- 
füllt werden können. Diese Aufgabe ist mit den Mitteln der 
Analysis ohne weiteres zu bewältigen; allerdings werden bei 
dieser Aufgabe nur solche Funktionen zur Lösung gebraucht, 
die von der Azimutalkomponente unabhängig sind. In dieser 
Hinsicht sind die Probleme der Schwingungserregung wesent- 
lich weniger kompliziert als die Beugungsprobleme, bei denen 
Lösungen der Maxwellschen Gleichungen gesucht werden 
müssen, die auch von der Azimutalkomponente abhängig sind. 
Solche Funktionen hat für den Fall des Rotationsellipsoides 


> 


1) Lord Rayleigh, Phil. Mag. 43. 8.159; 44. 8.28; Scientific 
Papers 4. S. 283, 305. 

2) R. Gans, Ann. d. Phys. 37. S. 881. 1912; 47. S. 270. 1915. 

3) G. Mie, Ann. d. Phys. 25. S. 377. 1908. 

4) M. Abraham, Wied. Ann. 66. S. 435. 1898; Math. Ann. 52° 
8.81. 1899. 


5) F. Ehrenhaft, Wiener Ber. 113, S. 273. 1904. 
Annalen der Physik. IV. Folge. 83. 
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Niven!) aufgestellt, der zu ihnen geführt wurde, als er ver. 
suchte, die Wärmeströmung in einem homogenen Körper, dessen 
Oberfläche ein Rotationsellipsoid ist, zu berechnen, eine Auf. 
gabe, die noch Heine?) ungelöst zurücklassen mußte. Dem 
Wesen seines Zieles entspricht es, daß Niven nur solche 
Lösungen suchte, die im Endlichen endlich sind. Diese werden 
in der vorliegenden Arbeit als Grenzfälle wesentlich allge. 
meinerer Funktionen reproduziert werden auf einem Wege, 
der einwandfrei ist, im Gegensatz zu dem, auf dem Niven 
seine Funktionen findet. Gegen die Nivenschen Unter. 
suchungen kann man vom mathematischen Standpunkte aus 
dieselben Einwände machen, die seinerzeit Dannacher?) und 
nach ihm Mary Curtis‘) gegen den von Heine geführten 
Existenzbeweis für die Funktionen des elliptischen Zylinders 
erhoben. Aber nicht nur in dieser Hinsicht bedürfen die 
Nivenschen Untersuchungen einer Erweiterung, sondern auch 
die Integrale, die im Endlichen einen Pol haben, müssen ge- 
sucht werden. Auch diese werden sich als spezielle Fälle 
gewisser Funktionen im folgenden ergeben. 


a Besonders wichtig und vereinfachend ist es, wenn man 
die Integrale der Wellengleichung so in Reihenform darstellen 
Pg kann, daß ein jedes Glied dieser Reihe für sich der Wellen- 
gleichung genügt. Solche Funktionen, die selbst Integrale 
der Wellengleichung, außerdem gewissermaßen die Elemente 


der allgemeinen Integrale darstellen, nennen wir Normalintegrale. 

K.F. Herzfeld®) hat zum erstenmal den Versuch unter- 
nommen, die Beugung an Rotationsellipsoiden zu berechnen. 
Es gelingt ihm nicht, Normalintegrale der Maxwellschen 
Gleichungen zu finden. Die von ihm aufgestellten Reihen ge- 
nügen nur in ihrer Gesamtheit den Maxwellschen Gleichungen, 
nicht aber jedes einzelne Glied. Die Aufstellung der Reihen 
ist außerordentlich kompliziert, ihre Diskussion nahezu un- 
möglich und auch von Herzfeld nicht versucht. = 


1) C. Niven, Trans. of the Royal Soc. 171. S. 117. 1880. 

2) E. Heine, Handbuch der Kugelfunktionen II, Verlag von 
Reimer. Berlin 1881, S. 328 ff. 
8) K. Dannacher, Inaugural-Diss., Zürich 1906. 
_ 4) Mary Curtis, Annals of Math. 20. S. 23. 1918. 
ae 5) K. F. Herzfeld, Wiener Ber. 120. S. 1587. 1911. 
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Sudhansukumar Banerji?) beschäftigt sich mit der Inte- 
gration der Wellengleichung in elliptischen Koordinaten. Diese 
führt er ein durch die Gleichungen 


z=a0c08% z=c-osindsing. 


Dann sucht er Lösungen zu gewinnen, die denen in Polar- 
koordinaten ähnlich sind. Er muß jedoch zugeben, daß die 
von ihm gefundenen Integrale nur Näherungslösungen dar- 
stellen, die nicht einmal den Vorzug der Einfachheit haben. 

Es ist bis jetzt noch nie versucht worden, die Beugung 
an einem dreiachsigen Ellipsoid zu untersuchen. Diese Auf- 
gabe soll im folgenden durchgeführt werden. Dabei zeigt es 
sich, daß die Lösung für den Fall eines dreiachsigen Ellip- 
soides keineswegs theoretisch schwieriger zu behandeln ist, 
als der Fall des Rotationsellipsoides. In beiden Fällen führt 
die Auflösung der Wellengleichung auf lineare Differential- 
gleichungen, in denen über zwei verschiedene Parameter passend 
zu verfügen ist. Gerade diese Bestimmung der beiden Para- 
meter ist die wesentliche Schwierigkeit des Problems. Die 
Aufgabe wird dadurch gelöst, daß man jedem der beiden 
Parameter die Eigenwerte gewisser Integralgleichungen zu- 
ordnet und nachweist, daß jedem dieser Eigenwerte ein be- 
stimmter Parameterwert entspricht. Die Integrale der Diffe- 
rentialgleichung werden gewonnen als Lösungen dieser homo- 
genen Integralgleichungen, die durch Verallgemeinerung einer 
zuerst von Whittaker gebrauchten Methode?) aus einem all- 
gemeinen Integral der Wellengleichung abgeleitet werden. 
Eine Beschränkung auf Rotationsellipsoide wird durch die 
Theorie in keiner Weise nahegelegt; nur in der numerischen 
Auswertung werden wir uns auf diesen Spezialfall beschränken, 
weil allein dieser rechnerich durchführbar erscheint. Doch 
ist es nicht ausgeschlossen, daß nach Erledigung einiger mathe- 
mathischer Vorfragen auch der allgemeine Fall rechnerisch 
wird durchgeführt werden können. 

Die Maxwellschen Gleichungen befriedigen wir durch 
den Hertzschen Vektor. Dieser unterliegt nur der Be- 


1) Sudhansukumar Banerji, Bull. of Calcutta Math. Soc. 1915. 
2) E.T. Whittaker, Trans. 5. Int. Congress of Mathematicians, 
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dingung, daß seine rechtwinkligen Komponenten der Wellen- 
gleichung zu genügen haben. Aus ihm leitet man durch ein- 
malige bzw. zweimalige Anwendung der Operation „rot“ die 
elektrische bzw. magnetische Feldstärke her. Wir werden uns 
daher im folgenden mit der Wellengleichung befassen. Diese 
wollen wir nicht allgemein integrieren, sondern nur Lösungen 
suchen, die zeitlich periodisch sind. Diese Lösungen lassen 
sich darstellen als Produkte von zwei Faktoren, von denen 
der eine Faktor eine Funktion nur der Raumkoordinaten ist; 
der andere Faktor ist e'”'. Die Funktion der Raumkoordi- 
naten muß dann, wie man leicht sieht, der Differential- 
 gleichung 
AV+ = V=0 


der sogenannten Schwingungsgleichung, genügen. Es sollen 
Partikularlösungen dieser Differentialgleichung angegeben 
werden, die geeignet sind, die elektromagnetischen Grenz- 


bedingungen auf einem Ellipsoid zu erfüllen. = | 
L Teil 
Die Integration der Schwingungsgleichung in elliptischen 
Koordinaten 


§ 1. Einführung der elliptischen Koordinaten 
Unter elliptischen Koordinaten verstehen wir die drei 
Wurzeln der folgenden Gleichung dritten Grades in 9’: 
x? y? x? 


+ + 


. -—1=0; 0 
— p* <p<q 


q? 
x? (9 — p*) (9? — 9?) + (9? — 9°) + 2° — p’) 

— —p*) (9? — 9°) =9. 
Durch diese Gleichung wird o? als eine Funktion der drei 
räumlichen Koordinaten x, y, z dargestellt und zwar ist, den 
drei Wurzeln entsprechend, o? eine dreiwertige Funktion der 
Raumkoordinaten. Der Wert jeder Wurzel, oder wie wir 
auch sagen wollen, jeder der drei Zweige dieser Funktion, 
liegt zwischen genau angebbaren Grenzen. Einer der drei Zweige 
liegt zwischen Null und p?, wir bezeichnen diesen mit »*. 
Der zweite Zweig, den wir mit u? bezeichnen wollen, liegt 
zwischen p? und g?; der dritte Zweig schließlich ist immer 


oder 


größ 
i 
wir 
4 drei 
leitt 
oan 
> 
er; Hal 
+ 
De: 
Ge 
der 
du 
a 
El 
Da 
De 
‘ sh 
= 
ny 
4 Hs 
al 
_ 3 
di 
G 
| 
ge 
de 
d 
awe 


1 


Beugungserscheinungen an Körpern v. ellipsoidischer Gestalt 613 


größer als g? und hat keine obere Grenze. Diesen nennen 
wir 2°. 

Die Darstellung der kartesischen Koordinaten durch die 
drei Zweige 2, u, » wollen wir hier nur angeben, da die Ab- _ 
leitung derselben bekannt sein dürfte.‘) Es ergibt sich: = 


| | 


Es sei nun ein Ellipsoid mit den Halbachsen /, m, n ge- 
geben. In seinem Mittelpunkt befinde sich ein rechtwinkliges 
Koordinatensystem, dessen X-Achse mit der Richtung der 
Halbachse 7, dessen Y-Achse mit der Halbachse m und dessen ~~ 
Z-Achse mit der Halbachse n des Ellipsoides zusammenfallt. 
Der Polabstand 9 sei der Winkel zwischen einer beliebigen the Be; 
Geraden durch den Nullpunkt des Koordinatensystems und = 
der X-Achse. Zur Festlegung der Richtung einer Geraden 
durch den Nullpunkt bedienen wir uns dann noch des ndr 
Y-Z-Ebene befindlichen Azimutes gm. Ein Punkt auf dm 
Ellipsoid hat die rechtwinkligen Koordinaten ?) ete. 2 


r=lcoost; y=msndcosg; z=xnsindsing. 


Da nun andererseits nach den Gleichungen (1) bzw. nach der 
Definitionsgleichung der elliptischen Koordinaten ein Punkt 
mit den Koordinaten z, y, z auf einem Ellipsoid mit den 
Halbachsen 


1) Vgl. z. B. Heine, Handbuch der Kugelfunktionen I. 8. 352. 
Berlin 1878. 

2) Diese Winkelkoordinaten 9 und haben allerdings nicht mehr 
die einfache geometrische Bedeutung wie die gewöhnlichen Winkel- 
polarkoordinaten. Sie sind nicht mehr die Richtungswinkel derjenigen 
Geraden, die durch den Nullpunkt und den fraglichen Punkt hindurch- 
geht. Ihre Bedeutung ist vielmehr die folgende: Durch den Punkt auf 
dem Ellipsoid und die zu den konfokalen Ellipsoiden gehérigen beiden 
orthogonalen Hyperboloidenscharen ist ein Hyperbelast bestimmt, der 
das betrachtete Ellipsoid in dem betrachteten Punkte senkrecht durch- 
schneidet. Die Richtungswinkel des einen Asymptotenstrahles sind ge- 
rade die Größen $ und gy. 
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liegt, so erhalten wir aus (1) folgende Gleichungen: 


p> 

Die elliptinchen Koordinaten haben die Eigenschaft Pr 
Orthogonalität. Dies zeigen wir durch die Berechnung des 
Bogenelementes, dessen Quadrat sich als rein quadratische 
Form ergeben wird. Zu diesem Zweck differentiieren wir die 
Gleichungen (1) logarithmisch und erhalten: 


cr = sin Fcosg = 


we 


sin # sin yp = 


du-u dv.v 
i? — p? p® 
di-A | duru 


Hier quadrieren wir und addieren dann. So erhalten wir a 
Quadrat des Bogenelementes ds: 


(u* — p*) (u? — 99) 
(v? — 12) — u®) 
(»* — p*) (»* — 9) 
Denn die Koeffizienten der doppelten Produkte 


di-dp; di-dv; dp-d» 


a 3 verschwinden identisch, wie man leicht nachrechnet. 

Pets Von den Koordinaten A, u, » können wir nun zu anderen 
elliptischen Orthogonalkoordinaten übergehen. Zum Beispiel 
können wir solche durch folgende Gleichungen definieren: 

a=, (1); 7 = 

_ Diese Koordinaten müssen ebenfalls als elliptische angesprochen 
werden, denn aus «= const folgt 2 = const usw. Also auch 
- @=const stellt im Raume der zx, y, z Ellipsoide dar, ebenso 
= const und 7 = const ein- bzw. zweischalige Hyperboloide. 
cad as Wir verfügen nun über die Funktion g,, indem wir be- 
stimmen: y, sei das elektrostatische Potential, das von einem 
der Ellipsoide A = 2, hervorgebracht würde, wenn dies eine 
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elektrostatische Ladung hätte, Ein solches Potential ist durch 
ein elliptisches erster — dargestellt, namlich’): 


Indem wir über einen willkürlichen Faktor und durch Fest- 
legung einer der Grenzen des Integrales über eine an sich un- 
wesentliche additive Konstante verfügen, führen wir an Stelle 
yon 4 eine neue Variable & durch folgende Gleichung ein: 


do 


A 


q 


k’ 
(4) 


Suchen wir die zu « inverse Funktion, d. h. versuchen wir A/q 
durch « auszudrücken, so werden wir bei der von uns ge- 
troffenen Wahl der unteren Grenze des Integrales (4) auf die 
Jakobische Funktion dn («, k) geführt, die außer von «noch von 


V®-p R 


elliptischen Funktion, abhängt. Die elliptischen Funktionen 
sn(@, k) und, cn(e, k), die im folgenden auch gebraucht werden, 
hängen mit der eben eingeführten dn(«, k) durch folgende 
Gleichungen zusammen: 


dem Parameter 4 = dem sogenannten Modul der 


(4a) sn(a,k)=— VI — A); one, 


k nennen wir den komplementären Modul. 


In derselben Weise wie A und & miteinander verknüpft sind, 
sollen auch die neu einzuführenden $ und y mit bzw. u und » 
zusammenhängen. Wir erhalten: 

6) A=qdn(a,k); w=gdn(B,k); »=qdniy, &). 

Für spätere Entwicklungen ist es vorteilhaft, anzumerken, daß 
die durch (4) definierte Größe & als Funktion der kartesischen 
Koordinaten x, y, z der Laplaceschen Differentialgleichung 
genügt: 
(6) 

Die Gleichungen (5), durch die A, u, v einerseits mit den 
Größen a, ß, y andererseits verbunden werden, reichen zur 


1) Vgl. z.B. M. Planck, Einf. in die Theorie der El. und des 
Magn. S.56. Leipzig 1922. Doch beachte man, daß wir A? nannten, 


was bei Planck a? +4 ist. Ebenso entsprechen unsere p, q den 
Planckschen Größen a? — b? bzw. a? — 
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F. Moglich 
eindeutigen Bestimmung noch nicht aus; denn wir hatten ge- 
zeigt, daß die Größen A, u, v gemäß ihrer Definition nicht 
beliebiger Werte fähig sind, sondern auf bestimmte Intervalle 


beschränkt werden mußten. 


Dies hat zur Folge, daß auch 
die Größen a, 8, y gewissen Beschränkungen unterworfen sind, 
die sich aus denen für A, u, » eindeutig ermitteln lassen 
müssen. Betrachten wir z. B. die Größe u; diese war auf ein 
reelles Intervall beschränkt, dessen obere Schranke g, dessen 
untere Schranke p ist. Setzen wir nun in Gleichung (5) 8 =0, 
so nimmt dn(0) den Wert 1 an, und wir erhalten die Gleichung 
u=g. Da g die obere Schranke für u ist, so haben wir in 
ß=0 sicher eine Schranke für die Variable # gefunden. 
Setzen wir andererseits 6 =— K, wo K die reelle Viertel- 
periode der Jakobischen Funktionen darstellt, so ergibt sich 
die Gleichung: dn(— K)= Kk =p/q 


und daraus der Wert p für die Größe u. Wir haben also 


als untere Grenze für die Variable @ den Wert —Ä, als 
obere Grenze den Wert 0 gefunden. 


Weiter ist die Funktion 
dn (u) für reelle Werte von x reell und im Intewail 
—K<u<0 
monoton wachsend; daraus folgt, daß das Intervall für die 
Variable 3 die reelle Achse zwischen — X und 0 ist. 

— K ist offenbar auch eine Schranke für die Koordinate 7. 
Denn y korrespondiert der Koordinate v, und das Maximum 
von » ist gerade der Wert p. Gehen wir nun auf der reellen 
Achse von — K aus nach —2K, so nimmt dn (u) für eine 
solche Veränderung seines Argumentes dauernd zu, bis es 
für v=—2K den Wer 1 wieder erreicht hat. Aus den periodi- 
schen Eigenschaften der Funktion dn(u) erkennt man, daß dn(w) 
für reelle Werte seines Argumentes nur Werte zwischen A’ undl 
annehmen kann. Sollen daher die Gleichungen (5) richtig sein, 
so müssen wir für die Koordinate 7 komplexe Werte zulassen. 
Dabei werden wir das Gebiet der Variablen y aber so ab- 
grenzen, daß der Wert von dn(u) reell bleibt. 

Die einzigen Werte von u, für die die Funktion dn(u) 
verschwindet, sind die Werte v= (2m +1) K+i(2n +1)X, 
wo m und 2 ganze positive oder negative Zahlen bedeuten 
und iX’ die imagine Viertelperiode der Jakobischen Funk- 
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tionen ist. 2 ist die imaginäre Einheit. Unter den Punkten 
von u, die dn(u) zum Verschwinden bringen, müssen wir 
nun denjenigen aussuchen, der so mit dem Punkte u =— K 
verbunden werden kann, daß längs einer bestimmten Kurve, _ 
die die beiden Punkte verbindet, dz(u) erstens nur reelle 
Werte und zweitens, Werte zwischen 0 und A einmal und. 
nur einmal annimmt. Man sieht leicht mit Hilfe der Additions- 
theoreme der Jakobischen Funktionen und mit Hilfe der © 
imaginären Transformationen von Jakobi!), daß der einge 
Punkt, der allen beiden Bedingungen genügt, der Punkt 
—K+iK’ ist, daß ferner die Kurve, längs deren dn(u) reell 
ist, durch diejenige Gerade dargestellt ist, die durch die beiden 
Punkte — AK und — X hindurchgeht. Für y finden wir 
also als Intervall das Stück einer Geraden parallel zur imagi- __ 
nären Achse, das durch die Punkte -— X und —K+iK’ be- | 
grenzt wird. Durch genau die gleiche Überlegung finden wir — 
für & als Intervall das Stück der imaginären Achse, ds r. 
zwischen 0 und 7 KX’ liegt. N 
Wir werden nun auch die Gleichungen (1) und (2) auf | 
die elliptischen Koordinaten @, 8, y umrechnen, wobei ir 
eine bemerkenswert symmetrische Form erhalten werden. Mt 
Benutzung der Gleichungen (4a) u sich leicht: See ae 


r= dnadnB dny; y= enaenBeny; 


z= 


Die Gleichungen (8) legen eine Erweiterung derjenigen — 
Intervalle nahe, die wir für die Variablen 9, y gefunden haben. ro 
Denn bei der eben eingeführten Beschränkung sind die Größen 
dnßdny; icnBeny; snßsny entweder nur positiv oder 
nur negativ. Da aber der Polabstand # zwischen 0 und n, 
das Azimut g zwischen 0 und 2 variiert, so werden aufden 
linken Seiten von (8) Vorzeichenwechsel auftreten, die durch 
Veränderungen der $# und y innerhalb der gesetzten Grenzen 
nicht reproduziert werden könnten. 


1) Vgl. z. B.: E. T. Whittaker, G. N. Watson, A course " oe 
modern Analysis, S. 491ff. Cambridge 1920. thin 
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Betrachten wir zunächst die Gleichung 4 
cos} dnßdny. | 


Setzen wir hier B = 0, y =— K, so erhalten wir wegen 
dn(— K) =k’, dn(0)=1, cos =1; also entspricht dem Werte 
# = 0 das Wertepaar 8 = 0, y =—X. Es ist leicht zu sehen, 
daß beim Festhalten an den vorhin eingeführten Intervall- 
begrenzungen das Produkt dnßdny immer kleiner als X ist 
und nur dann verschwindet, wenn y = — K+iK’ ist, da dng 
nicht verschwinden kann. Es entspricht daher dem Wert 
y=—K+iXK’ in eindeutiger Weise der Wert 9 = > Um 
nun auch negative Werte darstellen zu können, erweitern wir 
das Intervall fiir y über —K +i KX’ hinaus bis zu —K-+ 2ik’. 
Dann müssen wir allerdings das Intervall für # so erweitern, 
daß in ihm dn(f) nie negativ wird. dn(u) ist nun nie negativ, 
sondern liegt immer zwischen 1 und 4’, wenn « reell ist. Wir 
dürfen daher das Intervall für $# nur längs der reellen Achse 
erweitern und haben zu prüfen, ob diese Beschränkung ge- 
stattet, auch sin#cosp und sin #singp vollständig durch 
und y darzustellen. 

Die Größen sin #cosyp und sin # sin verschwinden beide 
fir $=0,n. entsprach aber dem Wertepaar = 0, 
y =— K, und für dieses Wertepaar verschwinden in der Tat 


icnßeny und sn?sny. Dem Werte + = > entsprach der Wert 
y =— K+iX’, so daß wir erhalten: 


Gemäß den Gleichungen (4a) ist aber 

wie es sein muß. Weiter sieht man, daß icnßcny und sn sny, 
wenn y durch den Wert —K +X’ hindurchgeht, ihr Zeichen 
nicht wechseln, da 


enu— K+ik')= 


csp=—enß; inp=—snß, 

da N 
en(—K+ik)=—-, sn(— K+ ik’) =— 


kenu 
und en(u) = en(—u); dn(u) = dn(—u) 

ist. Die in ienßeny bzw. sn sny auftretenden Zeichenwechsel 
men also den Nullstellen von cos pm bzw. sin p entsprechen. 
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Da nun die Nullstellen von cnß an den Stellen (2xn+1)K 
liegen, (bei deren Überschreitung immer Zeichenwechsel ein- 
tritt), ferner die von sn 8 an den Stellen 2n X, so erhalten wir 
die notwendige Anzahl von Zeichenwechseln dann, wenn wir 
das Intervall von f wie folgt erweitern: —K <8 < 3K. Wir 
erhalten also für y ein Intervall der Länge einer halben ima- 

fir # die Länge einer ganzen reellen Periode: 

-K<B<+B3K. 

Es bleibt nun noch übrig, das Quadrat des Bogenelements 
auf die Koordinaten «, 8, y umzurechnen. Aus Gleich. (4) folgt 


da* (i? — — ). 
Analoges für du? und d»*. Unter Berücksichtigung von 
(3) erhalten wir: ae 


ds? = — q*{(dn*? — dn* 8) (dn? a — dn? 
(dn? B — dn*y) (dn? — a) dB 
+ (dn? y — dn?a)(dn?y — dn* B)dy*t. 
Es muß hier zunächst verwunderlich erscheinen, daß das 
Quadrat des Bogenelementes nicht, wie zu erwarten, positiv, 
sondern negativ definit ist. Der Grund hierfür ist der folgende: 
Die Größen @ und y ändern sich längs Parallelen zur imagi- 
nären Achse der a- bzw. y-Ebene. Daher sind da und dy 
rein imaginäre'Größen, deren Quadrate natürlich negativ sind, 
und die so die Kompensation des Minuszeichens bewirken. 
dß ist reell, jedoch ist der Koeffizient von df? negativ, weil 
nach Voraussetzung dn?« > dn? ist. Daher ist der Wert 
von ds? immer positiv, wie es auch sein muß. 

Wir hätten durch imaginäre Transformationen die Mög- 
lichkeit, zu reellen Koordinaten überzugehen, so daß sich in 
diesen auch die quadratische Form des Bogenelementes als 
positiv definite Form schreiben läßt. Diesen kleinen Vorteil 
würden wir jedoch durch den Verlust der bemerkenswerten 
Symmetrie in den Formeln (7) und (8) erkaufen, eine Symmetrie, 
die uns erst in einem der folgenden Paragraphen von großem 
Nutzen sein wird. Wir haben es daher vorgezogen, zunächst 
die Größen a, 8, y in der hier eingeführten Form beizubehalten. 


(9) 


WE 
‘a 
Im 
Ba 
7 
| 
| 


F. Möglich 


§ 2. Die Differentialgleichung 4V V=0 
in elliptischen Koordinaten 
ars Wir finden durch Differentiation: BT 
de da 98 da + dy du’ 
Op Oy 02? 0a08 \dx)\öx 


(55) (55) + (oe) 


Bilden wir dieselben Ausdrücke für die Differentiation nach y 


und z, so finden wir, daß der Ausdruck AY 1. die Differen- 
tialquotienten 14 a6 a nicht enthält; denn deren Koeffi- 
RR: zienten sind bzw. de, Aß, Ay, die gemäß Gleichung (6) ver- 
0V 

daöß’ Gady’ öß Bay 
2 fallen heraus wegen der Orthogonalität der «, ß, y; 3. die 


schwinden; 2. die Differentialquotienten 


0?:V 0?V 
Koeffizienten von dat? 


(ss) + + (32) 


Wir erinnern nun daran, daß wir die Größen a, 8,7 
als Potentiale eingeführt haben. In dieser Terminologie fort- 
fahrend können wir sagen, daß die Differentialquotienten, die 


da de du ay 

0a’ Oy’ O%’ 0% 

die rechtwinkligen Komponenten des Gradienten von bzw. «, 
ß, y sind, und so können wir sagen, daß die Koeffizienten von 


= nichts anderes sind als die Quadrate der 
absoluten Beträge dieser Gradienten. Nun läßt sich aber der 
absolute Betrag eines Gradienten sehr leicht durch einen 
Differentialquotienten ausdrücken; bezeichnet nämlich s, die 


Richtung der stärksten Änderung eines Potentiales g, so ist 


sind bzw. 


auftreten, nämlich 
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Fragen wir nun in unserem Falle, durch welche Eigenschaft 
die Richtung des stärksten Abfalls von z. B. & gekennzeichnet 
ist, so sehen wir sofort, daß in dieser Richtung sich die beiden 
anderen Koordinaten nicht ändern; es ist also längs s,: Re 
0 0 
Das Quadrat des Bogenelements in der Richtung s, ist also 
nach Gleichung (9), da d8 = dy = 0 sind, gegeben durch: 
ds? =—da*.q*-(dn*a@ — dn* f)(dn? a — dn’ y). 
Und daraus dann: 


Oa\2 6a\2 1 
(53) + (55) + (52) (dn? « — dn? 8) (dn? — dn? y) 
Analoges für @ und y. 
Nunmehr bereitet es keine Schwierigkeiten, die Differential- 
gleichung AV’ + =0 auf elliptische Koordinaten umzu- 
rechnen; indem wir aus dem Ausdruck für 4 V den Faktor 


1 
~ — dn? B) (dn* a — dn?y)(dn?$ — dn?y) 


herausziehen, erhalten wir 


1 
q?(dn?« — dn* B)(dn* a — dn* y)(dn* 8 — dn?y) 
. Kan: (dn?y — dn2a) 
4] du: Y 

5 V 


+ (dna — 0 


oder, indem wir mit diesem Faktor herauf ah CRUE 


2 2 25) 
(dn? — dn?® + (dn y — dn*a) 


(10) + (dn? 
= w? (dn* — dn 8)(dn? a — dn?y) 


wo wir x-g = gesetzt haben. 


$3. Die linearen gewöhnlichen Differentialgleichungen, 
die aus Gleichung (10) resultieren ; 
Wir suchen jetzt Lösungen der Differentialgleichung (10), 
die sich als Produkte von drei Funktionen darstellen lassen, 
derart, daß die eine Funktion nur von «, die zwitenurvonß — 
und die dritte nur von y abhängt, Sei 4-B-C ein solches 


> 
= 
2 
u 
| ty 
> 
Tee 
Dur 
| 
| B 
“4-9 
h 
er), 
ty 
cr | 
er 
ie 
st 


ist; ferner genüge A-B-C als Funktion von «, # und y der 
Differentialgleichung (10). Setzen wir jetzt # und y konstant, 


so erhalten wir fiir A folgende Differentialgleichung: N 
= (m, dn« + E,dn?a+ H,)- A 


sm, E,, H, Konstante, 
| 
dn®#+ H,)+ B, q 
dnty + E. dn? H C. 
at Wir multiplizieren die erste Differentialgleichung mit 
= er (dn? 8 — dn?y)B.C, die zweite mit (dn?y —dn?a)A-C, die 
er: = dritte mit (du®& — dn? 8) A- B; dann addieren wir. So erhalten 
wir auf der linken Seite, wenn wir noch A-B-C mit V be- 


= zeichnen: 
(dn? 8 — dn® (dn? y — dn2o)® + (dn?a — 

also gerade die linke Seite der Gleichung (10). Auf der he. 
2 m, dn* a(dn? 8 — dn?y) + m, dn* B (dn? y — dn? a) 


+ m, dn* y(dn?& — dn?ß) + dn®wdn?ß(E, — E,) 

+ dn* adn*y (EL, — E,) + dn?#dn?y(E, — E,) 

dn® B(H, — H,) + dn?y(H, — H,) + dn? a(H, — H,). 

: Soll das System der drei Differentialgleichungen mit (10) äqui- 

A valent sein, so müssen auch die rechten Seiten übereinstimmen, 

d.h. es müssen die mit Z, und H, multiplizierten Größen 

 herausfallen. Dies geschieht dann und immer nur, wenn 
E=Eb=E=E H=-H,=4,=H 

ist, wo E und H zwei zunächst willkürliche Konstante be- 

deuten. Rechts bleiben nur die Größen übrig, die mit bzw. 

m,, m,, m, multipliziert sind, und diese stimmen dann und 

nur dann mit der rechten Seite von (10) überein, wenn wir 

setzen 


Ws, 
DIR 


m =m, = m, = w?, 
Wir haben also ein Integral der Schwingungsgleichung 
durch eine einzige Funktion eines oT f, (w) in folgender 
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Der Unterschied im Verhalten der Funktionen wird durch die 
verschiedenen Intervalle der a, 8, y gewährleistet. 

Den weiteren Untersuchungen haben wir also die Diffe- 
rentialgleichung 


(11) — (w2dnty + Edn?y + H)f= 0 
zugrunde zu legen. 

Die Differentialgleichung (11) enthält zwei Parameter Z 
und H. Diese bestimmen wir so, daß die Integrale der Diffe- 
rentialgleichung (11) periodische Funktionen der Variablen w 
werden. 

Setzen wir » gleich Null, so geht (11) in die Differential- 
gleichung der Funktionen von Lamé über. Die Forderung 
periodischer Integrale der Differentialgleichung läßt sich dann 
nur erfüllen, wenn E=— n(n + 1) (n ganze positive Zahl)!) und 


H-—n(n+ 1) die Wurzel einer bestimmten algebraischen 
2 
Sei f(w) eine Lösung, der Differentialgleichung (11) und 


F(8,7) = f(8)- fly), so ersieht man sofort, daß F(#,y) folgender 
partiellen Differentialgleichung genügt: 
— [0X(dn*ß — dn*y) + — dn?y] F= 0. 
Diese Differentialgleichung hat den groBen Vorteil, nur noch 
eine Konstante, nämlich £, zu enthalten, während die Kon- 
stante H sich bei der Differenzbildung herausgehoben hat. 

Die Größen 8 und y sind durch die Gleichungen (8) mit 
dem Polabstand # und dem Azimut gp verknüpft. Die Funk- 
tion F kann daher auch als eine Funktion von 9 und @ an- 
gesehen werden, und wir wollen ermitteln, welcher Differential- 
gleichung sie als Funktion von 9 und zu genügen hat. Es 
ist also unsere Aufgabe, die Differentialgleichung (12) mit Hilfe 
der Gleichungen (8) umzurechnen. Hier sind die Umrechnungs- 
gleichungen: 


Gleichung vom Grade > oder wird. 


1) Nach Whittaker, Watson, a. a. 0. S. 564 lautet die Lamé- 

sche Differentialgleichung: 
aA 

dat 


was, wegen k?sn?’a = 1 — dn*a, mit unserer Behauptung übereinstimm 


= + ASA 


ay an 


3 
ag 
4 
<A 
| 
fi 
- 
4 
’ 
n Kor 
“an 
14 


Analoges für 


daraus dann: 


3 


testa 


Mödih 


er (35)' OF #9 OF dy 
| ) 68 
Aus der ersten der Gleichungen (8) folgern wir: 


sind; = snBenPdny 


und daher: 


0% 
Aus der zweiten und dritten der Gleichungen (8) ergibt sich: 

em 
tg ik enßenr’ cos*@m 66 eny 


dp 


tr anB(sn*® + en* 3) 
oder: 
dp O08. dp i OF 


sind Gy end Of 


Aus dieser Gleichung folgern wir sofort, daß 
a9 89 do 


ist und daß somit der gemischte Differentialquotient 


0? 
0389 
nicht auftritt. Aus denselben Gleichungen schließen wir auf 
ö® — sind sm?? 08 Oy Oy? 
Daher wird auch der Differentialquotient 07/Oq nicht auf- 
treten. Endlich können wir noch den Schluß ziehen, daß 


(38) (at) - 


Od Op . Od 


folgt ferner 


ist. Aus 
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0° 
ap =—isind 555, + cot ? 
024 
op ap ay + 


oder 


und damit haben wir gezeigt, daß 
or oF 


OF dy 
übergeht in 
Od oF OF 1 or 
- (3) + cot + sin? 


Man berechnet weiter leicht den Ausdruck ( 
und findet fiir ihn: 
(ae) 
SchlieBlich ergibt sich noch aus (8) 
3 + dn? y = cos? + sin? F cos* + *. 


0p OY 


2 
_ (5) = du? P— du?y. 


Damit haben wir auch das gesuchte Ergebnis gefunden: Die 
Funktion F, die in ihrer Abhängigkeit von den Koordinaten f 
und y der Differentialgleichung (12) genügt, muß als Funktion 
von + und ein Integral folgender Differentialgleichung sein: 


a ar 1 @F 
— [w? (cos? + k? sin? cos? gy) + =0, 


wo wir E fir E+ wm*k’* gesetzt haben. 
Für # fügen wir folgende Bedingung hinzu: Es sei E 
so bestimmt, daß die Integrale von (13) auf einer Kugelfläche 
endlich, eindeutig und stetig sind. Für w =0 würde daraus 
folgen: 
=—n(n+ 1). 
Denn für © = 0 geht die Differentialgleichung (13) in die be- 
kannte Differentialgleichung der Kugelflächenfunktionen über, A 
und es ist bekannt, daß nur dann Integrale dieser Differential- = 
gleichung mit den geforderten Eigenschaften vorhanden sind, __ 
wenn wir die Größe Z’ in der eben genannten Weise wählen. Er 
Annalen der Physik. IV. Folge. 88. 42 er 
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§ 4. Das allgemeine Integral der Wellengleichung 
und eine daraus folgende Integralgleichung fiir die Integrale 
der Differentialgleichung (13) 

Häufig wird in der Wellenoptik ein allgemeines Integral 
der Wellengleichung gebraucht, daß wir in folgender Form 
für unsere Zwecke benutzen wollen’): 


a 


(14) 2% = f + ysin + zsin iF, sin +ct) 
B=0 p-a d Go: 
Sucht man das allgemeine Integral der Wellengleichung fir 


zeitlich periodische Zustände, bzw. das allgemeine Integral der 
Schwingungsgleichung 4V+x2?/=0, so findet man leicht 


aus (14): 
(15) V= f fe + ysin cong, ; pdt, dy, 
a 


Bezeichne nun f eine Lösung der Differentialgleichung (11), 
dann ist #= f(f)- f(y) ein Integral von (12), bzw. als Funktion 
der Koordinaten # und „ ein Integral von (13). Als solches 
werde die Funktion mit Y(+,g) bezeichnet. Es ist dann 
sicher f(a)-F(,q) ein Integral der Differentialgleichung 
AV +x*?V =0, und wir können daher erwarten, daß sich 
f(@) p) in der Form 


(16) f(a) F(9,y) = f 
darstellen läßt. Nach den Gleichungen (7) und (8) rechnen wir 
jetzt das Integral auf die Koordinaten a, #, p um und finden: 


| 
1 6a) J 90060, + kenasin dcong sin 3) cos py + iksnasin dsin psin 


Po) d 


wo wiederum x-g gleich » gesetzt ist. Die zunächst will- 
kürliche Funktion setzen wir gleich sin 9, 
wo H($,%,) ebenfalls eine zunächst willkürliche Funktion be- 


1) E.T. Whittaker, Math. An 
Aun. d. Phys. 77. S. 685. 1925. 
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gleich, z.B. «=0. Es ist nun: 
dnO=1, cnO=1, snO=0. 


So erhalten wir fir (7, eine 
Form: 


(17) Fis, p) of farm 90039, + ksin p sin C08 qo) 
po) sin 9, d 9, 


Jetzt soll die willkürliche Funktion 4(9,,9,) so bestimmt | 
werden, daß die durch das Integral (17) dargestellte Funktion 
F(&%,g) der Differentialgleichung (13) genügt. Dabei wollen 
wir fordern, daB H(9,,9,) auf einer Kugelfläche (also für 
0<9,<r, -a<gp,<r) endlich, eindeutig und stetig ist. 
Wir sehen sofort, daß diese Annahme hinreicht, damit auch 
die durch (17) dargestellte Funktion F(#,g%) denselben Be- 
dingungen genügt. Wenn dann noch die durch (17) dar- 
gestellte Funktion die Differentialgleichung (13) befriedigt, so 
haben wir sicher damit den ersten Schritt zur Integration der 
Schwingungsgleichung in elliptischen Koordinaten getan. Er va 

Um die Bestimmung der Funktion Z(9,9,) durchführen 
zu können, beweisen wir einen Satz, der sich folgendermaßen 
formulieren läßt: 

Bezeichne der Operator D, „ die Operation 


1 
D;,,= sine OF (sin x 
1 2 
+ an? — (cos? sin? ao cos? 
so genügt die in (17) auftretende Funktion 
+ ksind cos psin — 


als Funktion der Variablen 9, gy, %,, gy, der Differential- 
gleichung 

Dy, = Da, le?) | 

wo D,,,, diejenige Operation bedeutet, die aus Dy, hervor- = 

geht, wenn man durch 9, und durch g, ersetzt. 

Wir führen den Beweis dieses Satzes durch Ausrechnen _ er og 

etwas allgemeiner, indem wir zeigen, daß die Funktion 


ei® (a 3 cos Hy +b sin Asin J cos p cos + ¢ sin 9 sin sin sin Go) 
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F. Möglich 


immer der Differentialgleichung Dy, ,,(f) = Do, ,,(f) genügt, 
er = wenn wir die Größen a, 5, c der Nebenbedingung 


— cc? = }? 


unterwerfen. 


2 


= eiwe = (ac0o8#cos + bsin & sin 9, CoSpCosp 
» @ 0 0 p fo 
+ csin sin 9, sin psin g,) 


= iw ee (— asin cos #, + bcos sin +, cos cos 
+ ccos sin 9, sin sin 
0° f 


= — w?eive(— asin + bcos sin cos cos y, 
+ ccos # sin d, sing sin p,)? — imo ee, 
ur iw e'»e(— bsin sin J, sin cos gy 
+ esin # sin 9, cos psin g,), 
y = — w?eive(— bsin J, sin 
+ esin #sin cos p sin y,)?—iw(o—acos cos 


sar — = — we're fa? (sin? F — sin? 
— [sin? — sin? - [? cos cos? + e? sin? sin? 
— 2b ce [sin? + — sin? sin p sin cos cos 


of p 908 Mo _ sin? 
cot —cot dt, e sin $ sin 9, [sin sın 
sin psin _ gain? 
+ [sin? 9 — sin 
1 of 


and, de) co? e'»e (sin? 9, cos? p, sin? 
— sin? 9 cos* sin* g,] + e?[sin? 9, sin? cos? 
sin* cos? g,] 

+ 2bcsin sin gy, cos cos p,[sin? — sin? 


11777) wie . 
sin $ sin 3, [sin? + — sin 
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Hier addieren wir und erhalten: 
ind 89 0g? 


1 of 1 
sind, OF, (sin ö + sin? 9, 


= — m? ja?[cos? — cos? #] — b? sin? cos? p, 
— b? sin? # cos? p sin? p, + b?sin? 9, cos? y, cos? gy 
+ 6? sin? 9, cos? p, sin? p — c? sin? sin? sin? 
— sin? cos’ + ce? sin? #, sin? sin? ¢ 

+ c?sin? sin? cos? 7} 

= — w?eive{a?[cos? — cos? +] + 5?[sin? F, cos? 
— sin? g] + ec? [sin? 9, sin? — sin? sin? 

oder wegen der Nebenbedingungen, die a, 5 und c zu er- 

füllen haben: 

(18) = Dam 

Mit Hilfe des Operators D,, schreibt sich die Differential- 

gleichung (13), der die durch (17) dargestellte Funktion F(:7, g) 

geniigen soll, in fog Form: 


Wegen (17) ist 


Ds, F(t g) = const ff D, yl” (cos cos +k sin sin cos cos 90] 


0 -z 


H(9,p,)sin 3, dH, dg, 
und wegen (18) I 
0 
p,)8in F, dt, 
Nun betrachten wir das Integral TE 
SS 0d, (sin 
1 ö2 ve 


Dies integrieren wir zweimal partiell nach 9, und zweimal 
partiell nach g,. So erhalten wir: 


Lat 
t -- 
> 
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re 
| 
= 
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Bar 
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N a 


x 
+S abs; H 0 Po 


erste der beiden einfachen Integrale verschwindet, weil 
sin fir %,=0 und z gleich Null ist, und nach 
a setzung H für reelle Werte von gy, zwischen —a und z 
er endlich, eindeutig und stetig ist. Das zweite einfache Integral 
4 verschwindet, weil der in der inneren eckigen Klammer stehende 
u Ausdruck in bezug auf 9, von der Periode 22 ist, und weil 
a2. % man die Differenz von zwei Funktionswerten zu bilden hat, 
deren Variable gw, den Unterschied 2” aufweist. Wir finden 

somit die bemerkenswerte Gleichung: 


x 
p 


Ds, Po (H (7, fo ) sin ad, 
Daraus und aus Gleichung (17) folgt: 


io (cos cosd, + k sin # sin cos ~ COS qo) 
1 { Do, iG 0 Po) — E. 0 sin ty dit, dy, . 
Soll nun Do, » — F(9, 9) = 0 
sein, d.h. soll F(:+,q) der Differentialgleichung (13) genügen, 
so ist offenbar dafür hinreichend, daß auch H(#q) der 
_ Differentialgleichung (13) genügt. Da wir forderten, 4 soll 
auf einer Kugelfläche endlich, eindeutig und stetig sein, so er- 
kennen wir, daß eine Lösung der Differentialgleichung (13), die 
den gestellten Bedingungen genügt, auch als Lösung folgender 
homogenen Integralgleichung gefunden werden kann: 


-n 
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} fo =0 (20) 
könı 
= wer‘ 
dit, 
(20: 
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4 (19) F(t, = iff ei (cos cos + k sin sin Dy cos Cos qo) I 
= he 
0 -: . 
F(t, Po) sin ty dit, dg, . 
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An die Stelle des in der Integralgleichung auftretenden Kernes 
(20) K, (iF, to, Po) = (cos # cos +k sin sin cos p COS Po) 
kénnen, wie man sofort sieht, auch zwei andere Kerne gesetzt 
werden, nämlich 
(20a) K,(9, 9; Po) om d,sin psin go) 
fo) = (k’ cos cos + ik sin dein sin ~ sin Po) 
denn man hat in der auf S. 627 auftretenden e-Funktion nur 
a=0,ah c=ib=h’ baw. b=0, dh c=ika=k zu 
setzen. 


§ 5. Reihenentwicklungen fiir die Kerne 
Die Weiterentwicklung der Theorie beruht auf einer 
Reihendarstellung der Kerne (20) und (20a), die nach Kugel- 
flichenfunktionen fortschreitet. Um zu dieser zu gelangen, 
leiten wir kurz einen Satz über die Form spezieller Integrale 
der Schwingungsgleichung ab. 
Wir suchen die Differentialgleichung 47 + x?V =0 da- 
durch zu integrieren, daß wir den Ansatz a alas 
V = Ole, y, 2) | 
machen, wo die Funktion # die Variablen z, y, z nur in der 
Verbindung x? + y?-+ z? enthält, und ® der Laplaceschen 
Differentialgleichung geniigen soll. Setzen wir 
so finden wir leicht, daß unter den gemachten Voraussetzungen 
A(R) + x*(R@) = 0 übergeht in 


@®R 
q 


Hier erinnern wir uns daran, daß die einfachsten Integrale der 
Laplaceschen Differentialgleichung gewisse homogene Funk- 
tionen der drei Variablen z, y,z sind. Für homogene Funktionen 
vom Grade n der Variablen z, y, z gilt nach dem Eulerschen 


Satz: ar ar ar 
+37, +25, n-F 


Nehmen wir also an, sei eine homogene, harmonische 
Funktion vom Grade n, so muß die Funktion R folgender ti Bi 
Differentialgleichung genügen: 


& 
+ a 
l 
l 
- 
ean 
l 
3 
‘ 
=e 
} tal 
| 
ie 


632 F. Möglich 
Hier setzen wir nun: ete 

| für 
= und finden dann fir f die Differentialgleichung: pil, 
2 (n + 3) 0 
Die Integrale dieser sind die Zylinderfunktionen der Ord- | ” 4 
nung + (x + 3) des Argumentes xo. Um die späteren Ent- die 
wicklungen einfacher darstellen zu können, ist es vorteilhaft, (23) 
gewisse Funktionen ,(w) durch folgende Gleichung zu de | ,, 
finieren: 
vor 
Bezeichnet “si weiter D, (x, y, z) eine harmonische, homogene 
Funktion vom Grade n, so ist sicher 
Wn (x V2?+y? (x, y z) Je 
Vary 
ein Integral der Schwingungsgleichung “a 
Nunmehr gehen wir aus von der Reihenentwicklung fir P 
eireosy, Ks bedeutet: 
cos 7 = cos cos + sin cos(p — GP). 
Für eir«sr ergibt sich die Reihe‘): 
COSY S (2n + 1)i" wy, (r) P_(cos 7). 
sQ 
Die Kugelfunktion P (cosy) läßt sich in 2n+ 1 spezielle 
zerlegen ?): F b 
P,(cosy) = P, (cos i) P, (cos #,) | m 
(22) 9 (n — m)! ” 9 p™ 6 
| +2 2 (cos #) (cos cos m (p — u 
wo 
= in™ dm P (e 
(cos = sın i} d(eos jm [ „(cos )] 
bedeutet. ( 
1) Vgl. z. B. Handbuch d. phys. Optik, Bd. I, 8.567 GI. (81b). 1926. 
Doch nennt der Verfasser dort vn, was er hier mit y,(w) be 
zeichnet hat. 
2) Z.B. Whittaker, Watson, a.a. O. S. 328. 
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Wir stellen uns die Aufgabe, aus der Reihenentwicklung 

für eireo®r eine Reihe für die Funktion 
ei cos du + y sin COS po + z sin SiN Yo) 

zu finden, die nach Kugelfunktionen der Variablen 9,%, fort- 
schreitet. Es müssen also die Größen r, +, m durch die 
z, y, z ersetzt werden. 

Verbinden wir die Größen r, +, g mit den z, y, z durch 
die Gleichungen 
(23) z=rc8% y=rsn#cosy; z=rsinge~ 
so ist: 

x cos it, + y sin cos p, + zsin sin yw, = rcos7. 
Wegen Gleichung (22) enthält die Reihe für ei”«r Glieder 


von folgender Form: ER 
) (cos 0) sin m @ 


Jedes einzelne dieser Glieder geniigt der Schwingungsgleichung. = 

Wir werden nunmehr versuchen, an die Stelle jedes der 
speziellen Integrale der Schwingungsgleichung für Polar- 
koordinaten ein ebensolches für die kartesischen Koordinaten 
zu setzen, das die Eigenschaft hat, für die spezielle Trans- 
formation (23) mit dem Integral in Polarkoordinaten überein- 
zustimmen. Mit anderen Worten, es soll sein: 

m cos m n (Va? + y? + 2?) 

sobald 7,9, p mit z, y, z durch Gleichung (23) verbunden sind. 

Eine solche Funktion herzustellen gelingt leicht, wenn wir 
berücksichtigen, daß jede Kugelflächenfunktion der Ordnung n, a 
multipliziert mit r", eine homogene harmonische Funktion vom yous 
Grade n in den zx, y, z wird, wenn Gleichung (23) erfüllt ist; 
und zwar gilt’): 
ime P™ (cos #) mp 


@, (x y 2), 


« 
2 
4 
“ 4 
« 


m +n)! 
2n-n! 


fet iycore cos m udu 


= (x + iycosv + izsinv)" sinmvdv. 
0 


B 
d- 
t- 
t, 
/ 
1e 
ir 
e 
1) Whittaker, Watson, a.a. O. S, 392. 


(2,9, 2)= — (+i cosmud»; 
0 
(25) 
S” (2, y,2)= x dv, 
ei (z cos Dy + ysin cos + z sin Do sin Yo) 
/ 2 2 
= 2” (2n + 1) + y 
) (xy 2)P, (cos 9) + 
r S™ (xyz) P™ (cos #,) sin m Go | | 
. ae A, ; Aus dieser Reihe läßt sich nun aber leicht die Reihe für 


Möglich 


Damit ist die Reihenentwicklung aber bereits gegeben, denn 


wir haben in der Reihe für eir«sr nur noch die Größen 


cos 


y, (7) P™ (cos 9) 
durch 


Bid +n)! f n C08 MV 
j - x iy COS v iz si - 
+ sin m © dv 
0 


zu ersetzen. Bezeichnen wir mit S”"(ryz); C™ (xyz) die 


n 
homogenen harmonischen Funktionen 


2x 


nur zu setzen: 
=wcot; y=wksindcosp; z=0. 


2a 


fe + iy cos v)" sin mvdv = 0 


den Kern A, 9, ableiten; wir haben zu diesem Zweck 


m in (24) eingeführt erkennt man sofort, daß CE 


ist. I 
die d 
Wir | 


Integ 


mit 


defin: 


(27) 


d 
zZ 
| 
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ist. Es fallen in der gesuchten Reihenentwicklung alle Glieder, 
die den sin eines Vielfachen von g zum Faktor haben, heraus. 
einem geeigneten Faktor multiplizierten 


Wir bezeichnen die mit 
Integrale 


> 


0 


mit @” (2, y), indem wir G, (x,y) durch folgende Reihen 7 
definieren: 


fe + iy cos v)'cosmvdv 


2m — 2m)! On 
Grey) = 
Qe! EN 
2 m)! (oe —m)! 
n 
1 2n 
= Der), 
o=0 
& 279-1 gt g! J 
n 
2m (2m — 2m + 1)! ofzr +1 
Gang (% 3) en 2 | 20 ) 
20! 


a?n-20+1 


(27) 


4m)! (om)! 


n 


0 1 2 Y 
Const (z, y) 2 (— 1) | 20 ) 
o=0 
20! 
_ g3n—2041 „2 
n 
(x a (en (en — 2m)! 1), 2n+1 
angi (2 + 1)! 20+1 
o=m 
2° (a +m+1)! (oem)! 
2m+1 (2m — 2m — — 2m — at 
o=m 
(20 + 1)! 


rt 
AF ta 
; 
Ber: 
< 
m | 
| 
> 


Mit Hilfe dieser Funktionen kénnen wir dann schreiben: 
K, (iF, » Po) = eiw (cos # cos + k sin sin cos p COS qo) 


(2n +1) y,(@ Veos? + + sin? IT 


Veos? + + k* sin? + cos? 


n=O 


G” (cos #, k sin «+ cos p) P™ (cos cos m 


ze ie Tu Übertragen wir diese Methode sinngemäß auf den Kern 
Ki, 9, 9, Po), 80 haben wir zu ‚setzen: 


| z=okcos% y=0 z=iwksin sing. 
Le aa in (24) eingetragen, erhalten wir die Integrale: 7 
2x 


Br Das erste verschwindet, wenn m ungerade ist, das zweite, 
_ wenn m gerade ist. So erhalten wir folgende Entwicklung: 


K, ($, p; Po) = cos cos Jy + ik sin sin d, sin sin qo) 


dren +1) wn (wo VK cos® & — sin? sin? 9) 
Vk * cos? $ — k? sin? $ sin? p 


n=0 


(29) {> (—1)" k cos it, tk sin sin p) 4 


- (cos 9,) cos 2m 
+ G2"*" ( cos 9, ik sin sin g) 


. (cos #,) sin (2m + 1) 


Eine ähnliche Entwicklung erhalten wir auch für den Kern 
K, (9 p, Fy Po): 


§ 6. Ein Exkurs auf das Gebiet der Integrale von Sonin 


Zu einer bemerkenswerten Integralformel, die als eine Er- 
weiterung gewisser Integrale von Sonin angesehen werden kann, 
gelangt man mit Hilfe der or des vorigen ii 


Beug 


De! 


unt 


D 
f 
e sın 
Wol 
nah 
nacl 
so € 
Inte 
| 
is 
in 
so 
In 
wi 
Ä 
2 
_ 
d 
Ar 
= 


Beugungserscheinungen an Körpern v. ellipsoidischer Gestalt 637 


Die Kugelfunktionen 
P* (cos 9,)ecosmy, und (cos 9,)sinmgp, 


haben die Eigenschaft der Orthogonalität, d. h. die Integrale 
a 2n 
f f P™ (cos 9,) P* (cos 
v 
sind nur dann nicht Null, wenn n=x und m=4 ist. 
Wollen wir eine Funktion /(9, g,), die noch gewissen, hier nicht 
näher zu präzisierenden Bedingungen genügen muß, in eine 
nach Kugelflächenfunktionen fortschreitende Reihe entwickeln, 
so erhalten wir die Koeffizienten dieser Entwicklung durch die 
Integrale 
a 2x 


Der Faktor a” ist, wenn m von Null verschieden, MD =e 


08 A Yo 
sin dG, sin dt, dy, 


. 
yi n (n + m) GE 
und wenn m gleich Null ist, 
22% 


Wollen wir die Funktion + k sin cos p sin d, cos 
in eine nach Kugelfunktionen fortschreitende Reihe entwickeln, 
so erhalten wir zur Berechnung der Koeffizienten folgende 
Integrale: 


a 2x 


SS ei@(sc0sd, + tsin 9, P™ (cos $,) ™ Pogin dit, dy, 
sin Mm Po 
0 


wo s= cosd,¢t = ksin #cosy gesetzt wurde. 
Die Integration über g, kann man sofort ausführen. 
Es ist 


eio(s cos 9, + t sin cos qo) — giwscost,, piwt sin 603 po, 


Den zweiten Faktor entwickelt man in bekannter Weise nach 
den cos der Vielfachen von g, und erhält: 


= J, (ot sin it) + (w ¢ sin 9) 
k=1 


Fx 
= 
4 
ER 
Dal 
3% 
"VE 
| 
+ 
vi 
| 
67, 
as 


vx 
| 


Ya 


638 F. Möglich 
Führt man dies in das Integral ein, so erkennt man, daß die 
Integrale, die sin der Vielfachen von 9, enthalten, ver- 
schwinden, während die anderen Integrale übergehen in: 

2in na” (w ¢sin #,) P™ (cos 9,)sind,d%,. 

0 

Reihen, die nach Kugelflächenfunktionen fortschreiten, sind 
eindeutig. Daher dirfen wir schlieBen, daB die in dem vorigen 
Paragraphen gefundenen Koeffizienten der Reihe (28) mit den 
hier gefundenen Integralen übereinstimmen. Also: 


f eioscos% J (wtsin P} (cos F,) sin F, d J, 
(n + m)! (s, Vs + t*) 
n—-m! 0 
Setzen wir hier in Gleichung (30) 2 = m, so ist wegen 
2 (cos #,) = sin” i+, 
und [gemäß Gleichung (27)] 
n ” . 
ef oscos J (wt sin dd, 
(30a) 0 
" 


Sonin und Gegenbauer finden folgende Gleichung'): 


f ei dy J _ (@ tsin sin’ a, d 


0 


(@ Vs? + 
= w,, ( @ (w Vs? - + [vgl. Gl.(21)], 


1) iN. Nielsen, Handbuch d. Theor. d. Zylinderfunktionen, S. 182. 
Leipzig 1904. af 


Beu 


so | 
das 
bat 
wei 
= “ Mm: 
A 
Eure unt 
ver 
(30 
Die 
tio. 
Wil 
. 
Wil 
vo) 
kö 
be 
un 
Di 
nu 
Me 
Hi 
un 
lie 
ge 
ei 
in 
= 2 - 2 2 
Jy (@ Vs +1 ). sc 
Setzen wir in dieser 
d 
ei 
| 
8] 


Beugungserscheinungen an Körpern v. ellipsoidischer Gestalt 639 


so geht sie in die Gleichung (30a) über. Wir können daher 
das Integral (30) als eine Verallgemeinerung des Sonin-Gegen- 
bauerschen betrachten. 

Setzen wir in (30) ¢=0 und m=0, so haben wir noch, 


weil a, nur halb so groß ist, wie a” sein müßte, wenn wir 


m = 0 setzen, mit zwei zu multiplizieren; dann erhalten wir aber 
unter Beachtung von (27): 
(800) P (cos sin = 2i" y, (ws). 

v 
Dies aber ist die bekannte Integraldarstellung für die w-Funk- BEE. 
tionen. 


n 


§ 7. Formale Auflösung der Integralgleichungen 

Wir gehen daran, die Integralgleichungen aufzulösen, die 
wir mit Hilfe der Kerne (20) und (20a) bilden können. Dabei 
wird sich ergeben, daß wir die Reihenentwicklungen für die 
von uns gesuchten Funktionen numerisch nicht aufstellen 
können — zu diesem Zweck werden wir eine andere Methode 
benutzen —, jedoch werden wir in der Lage sein, die Theorie 
und die Eigenschaften dieser Funktionen vollständig anzugeben. 
Die Auflösung der Integralgleichungen wird dagegen auch 
numerisch durchgeführt werden können, wenn der Wert des 
Moduls der elliptischen Funktionen gleich 1 oder Null wird. 
Hierauf jedoch werden wir erst wesentlich später eingehen 
und zunächst die Theorie der gesuchten Funktionen für be- 
liebige Werte von k durchführen, 

Die Eigenfunktionen der Kerne (20) und (20a) sind, wie wir 
gezeigt hatten, Funktionen, die auf einer Kugelfläche endlich, 
eindeutig und stetig sind. Solche Funktionen kann man nach 
einem Satz aus der Theorie der Kugelflächenfunktionen immer 
in Reihen entwickeln, die nach Kugelflächenfunktionen fort- 
schreiten. Wir werden mit Hilfe der Kerne (20) und (20a) 
zeigen, daß solche Reihen wirklich existieren. 

Eine genauere Untersuchung der Differentialgleichung (11) 
führte den Verfasser zu der Vermutung, daß die Integrale 
dieser Differentialgleichung sich in acht verschiedene Typen 
einteilen lassen müssen. Diese Untersuchungen sind hier, da 
sie nur vorläufigen Charakter tragen, nicht mitgeteilt; es schien 


f > 

ie = 
J 
d 

n = 
n 

Be 

) 
5 


- 


| 
TE 


P. Möglich — 


diese Unterlassung zweckmäßig, um den Umfang dieser Arbeit 
nicht allzusehr anschwellen zu lassen. Auch rechtfertigt sich 
dies im Hinblick darauf, daß diesen Untersuchungen keine 
Beweiskraft zugesprochen werden kann. 

Da ferner, wenn f(w) ein Integral von (11) mit den ge- 
forderten Periodizitätseigenschaften ist, das Produkt f(ß)- f(y) 
mit Hilfe der Gleichungen (8) auf die Koordinaten 9 und » 
umgerechnet, der Differentialgleichung (13) genügen muß, s0 
lag weiter die Vermutung nahe, daß sich auch die Inte- 
grale von (13) in acht verschiedene Typen einteilen lassen 
werden. Diese Typen können sich aber nur dadurch von- 
einander unterscheiden, daß die sie darstellenden Reihen nicht 
die gleichen Entwicklungsfunktionen aufweisen, was aber dann 
nur möglich ist, wenn in den Reihen nicht alle Kugel- 
funktionen auftreten, sondern nur solche, die z. B. gleiche 
Irrationalitäten aufweisen. Nun sind aber alle Zugeordneten 
der Kugelfunktionen, die einen geraden oberen Index haben, 
ganze rationale Funktionen, während die mit ungeradem oberen 
Index immer die Irrationalität Yl — x? aufweisen. Daher liegt 
es nahe, Reihen anzusetzen, bei denen die Zugeordneten ent- 
weder gerade oder ungerade obere Indizes haben. Ferner 
wird man auch gerade und ungerade untere Indizes unter- 
scheiden, schließlich auch Reihen als zu verschiedenen Typen 
gehörig betrachten, je nachdem sie nach 


P™ (cos #)cosmg oder nach P” (cos +) sin my 
fortschreiten. 
Wir legen die Integralgleichung 


= af je (cos # cos d, + k sin sin cos cos (iF, sin d d 
00 


& 


zugrunde und machen für f(7,q) den Ansatz: 


(31) £09, 9) = 4, +S (S 42" P2" (0089): cos2mq), 


n=1 m=0 


dessen Richtigkeit wir zu prüfen haben. Um dies durch- 
zuführen, benutzen wir die für den Kern 
ei (cos Pcosd,+ksind sindcospeosg) 
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gefundene Reihenentwicklung (28), indem wir in ihr aber # 
mit 9, und p mit g, vertauschen. Die Reihe hat also fol- 
gende Gestalt: 


@ 


Seng ™ eR) 


n=0 


m=0 
Diese Reihe als Kern in die Integralgleichung eingefihrt, 
ferner f (9,9) auf der linken Seite durch die Reihenentwick- 
lung (31) ersetzt, ergibt durch Vergleichung der Koeffizienten 
der Kugelflächenfunktionen auf der rechten und auf der linken 
Seite: 


Veos?3, + k*sin® 9, cos*g,)” 


2m n n(wV cos* 9, +k? sin? cos® 
(82) J ( 


. (cos Asin cos py) - Po) 9, dg, 


Auf der rechten Seite erscheinen aber zunichst auch Koeffi- 
zienten von Kugelfunktionen, deren untere bzw. obere Indizes _ 
nicht gerade sind. Diese sind: a 

Wr (@ Veos* +k sin? F,cos*g) 
SI (Veos? +, + k? sin? 9, cos? 9,)" 


c08%,; ksin cosy) + f(t, sind, dd 
0 Po o» Po oF Po 
oe (r gerade oder ungerade), 


az 2n 


0 0 


(Veos? 9, + k? sin? 9, cos? "+! 


(cost; ksin cos + f(F, Po) sin dH, d yp. 
Man zeigt sehr leicht, daß diese, wenn man für f(:+,, 9) 

den Ansatz (31) benutzt, verschwinden. Wir wollen den Be- 

weis nur für 4°"*' 


" geben, fiir das andere Integral verläuft — 
er analog. 


Die Funktion 2 enthält, wie man sich leicht aus der is Pr 
Definitionsgleichung (21) überzeugt, nur gerade Potenzen von z. 
Daher ist auch 

Annalen der Physik. IV. Folge. 88, 


> 


Ps 


it 
ch = 
ne 
> 
7) n an a 
ai". . 6. m 
{ > G (cos; Asin cosy) P,, (cos cos m yl 
- 
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4 
(Veos? 4, + k* sin? 9, cos! 
durch eine Reihe darstellbar, die nur gerade Potenzen von 
cos ‘+, und cos g, enthält. Die Funktion ~ ' 


(cos Asin ch, cos p,) 


aber enthält nur ungerade Potenzen von cos +, und von 
sin 7, cosg,. Das gleiche gilt demgemäß auch für das 
Produkt 


w,(o Veos? + sin? 9, cos? gp) 
(Veos? 4, + sin? 9, cos? 
Führen wir jetzt in (32) für f(#,,9,) den Ausdruck (31) 


ein, so würden wir immer zu Integralen der folgenden Form 
geführt werden: 


2m+1 ‘ 
(cos 9,; ksin 9, cos g,). 


a 2x 

sn 2n +1 2m+1 2 2u 
Sf sin?" +1 - cos?™+! - P,” (cos 
00 


-c082ugp, 
Diese Integrale sind aber immer Nall, und damit ist die 
aufgestellte Behauptung bereits erwiesen. In der gleichen 


Weise zeigt man das Verschwinden der Größe A” er 


Jetzt zeigen wir, daß die Reihe (31), in der die 4’. 
durch die Gleichung (22) bestimmt sind, gleichmäßig konvergiert 
für alle Werte von + zwischen Null und a und @ zwischen 
Null und 2x. Zu diesem Zweck machen wir eine Ab- 
schätzung der durch die Gleichung (32) dargestellten Größen, 
indem wir beachten, daß a, abgesehen vom Faktor 


(— 1)"(4n+ 
immer kleiner ist als das Maximum des Integranden, multipli- 
ziert mit der Maßzahl 4” der Oberfläche der Einheitskugel. 
Der Integrand ist ein Produkt aus drei Funktionen, und das 
Maximum eines solchen Produktes ist sicher nicht größer als 
das Produkt der Maxima jeder der drei Funktionen. Von 
der Funktion f(#,, 9,) wissen wir nur, daß sie als Lösung der 
Integralgleichung (19) endlich ist auf einer Kugelfläche und 
daher auf dieser ein Maximum erreicht, das wir ohne Be- 
schränkung der Allgemeinheit dieses Beweises gleich 1 setzen 
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können, da in f (+, g,) ein konstanter Faktor noch zu willkür- 


Yon nn erreicht für 


licher Verfügung frei ist. Die Funktion 
z=0 ihr größtes Maximum, nämlich den Wert 


Zur Bestimmung des Maximums von (cos #,; ksin 9, cos 


denken wir daran, daß nach den Entwicklungen des § 5 
die fragliche Funktion durch das Integral 


(2n—2m)! 

2 (cos 9, + ik sin F, cos p,cosu)?"cos2mudu 

0 

dargestellt wird, aus dem man unmittelbar abliest, daß 07 
2m)! 


(cos 9,; ksin 9, cp) | = 2. 
0 — 


(cos #,; ksin cos Py) | =1 


ist. Da weiter 


2a 
_ (—1)"(2n + 2m))) 
- [Teos + ésin cos(p — u)]?"cos 2m udu 
0 


ist, so folgt die ebenfalls gleich zu benutzende Ungleichung; 


=< (2n — 2m)! 
(2n)! 


| (cos #) cos 2m y 


Wir können jetzt folgende Ungleichung anschreiben, deren 
Richtigkeit man leicht mit Hilfe der eben gegebenen Formeln 
nachprüft: 

A” <= 4Ya 2 (4n41) 
+177 (2 +— 5) 

Das 2n-te Glied in der Reihe (31) ist nun aber eine 

Summe, nämlich: 


an : (2n — 2m)! 
(2m)! 
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F. Möglich 


Für diese Summe müssen wir ebenfalls eine Abschätzung 
finden. Eine solche ergibt sich aus der Ungleichung für die 
‘ und aus der fiir die Funktionen at „ (08 F) cos 2m wp. 


m=0 


Man findet so: , 


(2n — 2) — 4)... @n—2m) 


Ada +1) 


Jeden der Briiche, die unter dem Summenzeichen stehen, 


ersetzen wir durch einen größeren, nämlich durch in, 


£ Dann aber setzen wir noch m = n. So erhalten wir a: 


™ (cos ı#) cos 2m p | ihe 


is ot 


< 4 Va’ +1) wo" 
(ons +) 


In—1+2?r.n"t}, 


HN Offenbar brauchen wir nur noch zu zeigen, daß ~ Reihe 
g2n o-” 


konvergiert. Zunächst kürzen sich die Potenzen von 2 weg, 
bis auf eine im Nenner. Dann beachten wir, daß 


(2n + 5) - (2n +.) (2" -5) (22-2)... (n +5) 
16-3) 


ist. Im Zähler steht ein Produkt aus n + 1 Faktoren, von 
denen jeder gleich n ist. Im Nenner steht außer der Größe 
2II(a — 4) noch ein Produkt aus n +1 Faktoren, von denen 
jeder größer ist als n. Daher ist 
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a, 


konvergiert. Daraus folgern wir die Konvergenz der Maxima ~ 
der absoluten Beträge von 5 


denn der noch hinzukommende Faktor 


Ant stért die Kon- 
vergenz in keiner Weise, so daß die Behauptung, die Reihe (31) 
konvergiere gleichmäßig, zu Recht besteht. 

In derselben Art, wie wir es eben getan haben, können 
wir vorgehen, um zu anderen Lösungstypen der Integral- 
gleichung (19) zu gelangen, z. B. können wir zeigen, daß fol- 
gender Lösungstyp existiert: 


(83) 9) = > |D (os 


n=0 m=0 


Die hier auftretenden Koeffizienten ergeben sich analog 
den Überlegungen, die zur Gleichung (32) führten; das Er- 
gebnis ist 


= th(—1)"(4n +3) 


(Veos? 9, + sin? 9, 


0 0 
‘ G5”, , (cos 9; kein 9,6089) FI, 


Die Gleichungen (31) und (83) unterscheiden sich nur da- _ BER: 
durch voneinander, daß die Reihe (31) nach Kugelfunktionen ‘ 
mit geradem unteren Index, die Reihe (33) aber nach solchen 
mit ungeradem unteren Index fortschreiten. Für beide Typen 
führen wir ein gemeinschaftliches Symbol ein, indem wir als 
Funktionszeichen für (31) und (83) den Buchstaben ? wählen. 
Wir unterscheiden die beiden Typen voneinander, indem wir 
der der Reihe (31) entsprechenden Funktion einen geraden 


¢ 
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g Wir haben damit gezeigt, daB die Reihe a 
9 
“tem 
A 
> 
| 
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: 
x vn 
enor: 
an 
n 
e iis 
n 


Möglich 


Index, der der Reihe (33) entsprechenden Funktion einen un- 


geraden Index erteilen. Also: 
(Bla) > P2* (cos 
n=0 Lm=0 


n=0 

Aus derselben ee FE (19) finden wir noch zwei 

andere Typen, die beide als Kennzeichen den Buchstaben § 
erhalten sollen, nämlich: 


(35) 9) = 


n=1 


n—1 7 
> *' (cos #) cos(2m + 1)g|, 


m=1 J 


@ 


n 


(cos #) cos(2m + 


n=0 Lm=0 J 
Die Koeffizienten dieser sind durch ka BEN be- 


(cos ksind, COS Po) S, Po) sin F, did, dy,, 


2m+1 
V cos? + k* sin? 9, cos? ,) 
1)"(4n + 3)i ro COS” Po! 
es) Veos? 4, + sin? 9, cos? 
2m+1 
-G 4 1 (C08 cos p,)S, Po) SII, dH, dp, 


_ Dies sind alle Typen, die wir aus der Integralgleichung (19) 
| ableiten können. Es fehlen noch vier, von denen zwei sich 
ohne weiteres aus den bisher gefundenen unter Betrachtung 
Br der Differentialgleichung (13) herleiten lassen. Es fehlen gerade 
noch die Typen, die nach sin der Vielfachen von q fort- 
oa: _ schreiten. Diejenigen Reihen, die ungerade Vielfache des sin 
enthalten, können wir aber leicht aus den Funktionen 8, (9,9) 


und 8,,,,(9,%) herleiten. Um zu dieser Entwicklung zu ge- 
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langen, ersetzen wir in der Differentialgleichung (13) das Azi- 

mut 9 durch = — y. Der Differentialausdruck wird dadurch 

(cos? + Foos*g)+ 

Sie geht nämlich über in Fig 

| + (7 


89) ) sin? + cos? g) 


Hier müssen wir nun beachten, daß die von uns betrachteten 
Funktionen nicht allein Funktionen der Variablen # und 
sind, sondern noch von zwei Parametern abhängen, nämlich 
von » und k. Es wird manchmal wünschenswert sein, auch 
diese Abhängigkeit zum Ausdruck zu bringen, was wir dadurch 
tun wollen, daß wir das Parameterpaar », k hinter 9, y, 
durch ein Semikolon getrennt, in die Klammer einschreiben. 
Wenn wir also bemerklich machen wollen, daß z. B. die Funk- 
tion 8,,(9,g) auch von » und k abhängt, so schreiben wir 
9;0@,%). 

Wir denken uns nun in der Funktion 8, (9,9; 
Parameter A durch den Wert ik/k’ ersetzt, betrachten also oe 
ik 
sieht, der Differentialgleichung 

1 
sin? 09° 
(cos? (Fy sin? cos? 9) + f=0. 


die Funktion $,, (», p;o, Diese genügt, wie man sofort 


6° f of 


Weiter setzen wir in S, , (9 93 0, 7) noch an die Stelle von 


gy die Größe + — 9. Die Differentialgleichung, der diese 


Funktion genügt, ist nach dem Vorhergehenden auch ohne 
weiteres zu ermitteln. Wir brauchen dazu nur den Ausdruck (39) 


zu betrachten und in ihm & durch > und &’ durch V 1+ a 


zu ersetzen. Dieser Ausdruck ist aber wegen A?+h?=1 
gerade gleich 1/k’; daher genügt die Funktion 


Pace 
. 
l 
) 
as 
# 
4 


1 0? f 
sin?’ gq? 


die also, die aus S,, (9,9; dadurch daß man 


durch — gy, wk’ und k durch * od ersetzt, geniigt 


; a ebenfalls der Differentialgleichung (13). Die. so erhaltene Funk- 
tion ist mit keinem bisher bekannten Funktionstyp identisch. 
_ Wir bezeichnen diese Funktion mit 7, (%,q;,h). Es gilt 
die Gleichung 


y ‚ick 
Die Funktion 7, (9,9; @,k) hat folgenden Reihentyp: 3 
ae 
» +! +1 (cog F) sin (2m + 
> n=1 Lm=0 


Die Koeffizienten BS +! dieser Reihe hängen mit den Koeffi- 
Er der Reihe (85) durch folgende Gleichungen zusammen: 


(wk). 
Nr Hier ist durch die benutzte Schreibweise darauf hingewiesen, 
daß die Koeffizienten unserer Funktionen 


des Parameterpaares w,k bzw. »%,°- sind. 


T 
Dieselbe Überlegung, die wir eben für die Funktion 
(9, 9;®,k) angestellt haben, läßt sich auch mit genau dem 
gleichen Ergebnis auf die Funktionen &,41(9,g) anwenden; 


wir erhalten so die Gleichungen: Be: 


‚ik 
= - >90 


(=) (cos? + k? sin? cos? + -f=0. (45) 
we 
Be, Diese aber stimmt fast vollkommen mit (13) überein, nur steht Die 
dar: 
Tan 
ae, 
Hilf 
We 
hab 
daß 
leit 
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n=0 Lm=0 


(4) -> (cos9)sin (2 m + 1)p 


(45) (@, A) = (- (ok 7) 
Die gleiche Methode auf die Funktionen AR,,(9,y) und 
Ran+ı (9, p) angewandt, liefert keine neuen Typen. Es sei auch 
darauf hingewiesen, daß man die Funktionen 7,,(9,Y) und 
Tonz1(%, 9) aus einer Integralgleichung hätte herleiten können, 
die, nicht wie Gleichung (19) mit dem Kerne (20), sondern mit 
Hilfe des Kernes X, (9, 9; 9,,%,) [Gleichung (20a)] in derselben 
Weise gebildet ist. 

Von den vermuteten acht verschiedenen Lösungstypen 
haben wir sechs gefunden. Man wäre geneigt anzunehmen, 
daß sich die übrigen beiden Typen aus den Kernen (20a) ab- 
leiten lassen würden. Dies ist nur zum Teil richtig. Der 
Kern X, kommt überhaupt nicht in Frage, weil bei der Ent- 
wicklung dieses Kernes in eine Fourierreihe der Variablen » 
sin der geraden Vielfachen von 9 nicht auftreten, die beiden 
noch fehlenden Typen als Funktionen von p sich aber in eine 
Fourierreihe entwickeln lassen müssen, die nur sin der geraden 
Vielfachen von 9 enthält. Hingegen gestattet der Kern 
KR, (9,9; eine solche Entwicklung, doch sind, wie eine 
nähere Betrachtung zeigt, die Koeffizienten von sin2mq nur 
gerade Potenzreihen von cos 9, so daß wir aus K, (7, p; +, Po) 
wohl Reihen erhalten, die aus Kugelfunktionen 


P}" (cos 9) sin 2m p 


mit geradem unteren Index bestehen, nicht aber solche, bei 
denen die Kugelfunktionen einen ungeraden unteren Index 
haben. Es besteht daher hier noch eine Lücke, deren Auf- 
fillung der folgende Paragraph gewidmet ist. Dabei ist es 
zweckmäßig, den Existenzbeweis für den siebenten Funktions- 
typ, der an den Kern K, (it, 9; %,,Y,) angeknüpft werden 
könnte, erst im nächsten Paragraphen zu führen, weil uns die 
folgenden Entwicklungen einen neuen Kern liefern werden, 
der zur Durchführung des Existenzbeweises geeigneter ist als 
der Kern A, (9595 For 

ah 
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88. Eine neue Integralgleichung, 


der gewisse Integrale der Gleichung (13) zu genügen haben 
Wir die Schwingungsgleichung 


Es sei f(z,y,2, u, = f(xz,y,0) ein Integral dieser Differential- 
gleichung, das die Variablen z,u,v nur in der Verbindung 
enthält. Dann genügt f(z,y,0) als Funktion von z, y, 9, der 
2 2 


Diese aber können wir auch, da z, y und o voneinander un- 


ea sind, so schreiben: 


Damit haben wir = Satz: Ist f(z,y, 0), wo o die Wurzel 
Vz? + u? + v? bedeutet, ein Integral der fünfdimensionalen 
Schwingungsgleichung (46), so genügt die Funktion ®=z - f(z,y,z) 
in ihrer Abhängigkeit von z, y, z der dreidimensionalen 


0° 0° 
Angenommen, wir hitten ein Integral von (46) in der 
Funktion f(x,y,o) gefunden. Wir ersetzen dann x durch die 
Größe g-dna-cosi, y durch gken«sin %cosp und o durch 
iqgksnasin sing [vgl. Gleichungen (7) und (8)]. Dann ist: 


iqgksna> sin 
- f(qdnacos 3d; qkenasin cos p;igksnasin sin 


ein Integral der auf die Koordinaten &, +, gm umgerechneten 
‘Wellengleichung. Schließlich sei angenommen, die Funktion 
f(qdnacosd; gkenasin cos p;iqgksnasin sing) zerfalle 
in ein Produkt von zwei Funktionen, von denen die eine nur 
- von «, die andere nur von * und p abhängt: 


f(qdna.-cos F; qkenasin cos p; igksnasin sin 


= Y (a) p)- 


52 6? ( (o 9 


Differentialgleichung = 


Schwingungsgleichung 
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Dann ersieht man leicht, daB dn«FY(«) der Differential- 
gleichung (11), sin gp) aber der Differentialglei- 
chung (13) genügen muß. Daraus folgt auch die Differential- 
gleichung, der O(:+,q) allein zu genügen hat, nämlich: 


| — (cos? + Hsin? + Z]O=0; 


(w + ixcos + iysin #,cosp, + iz sin sin p, COST, 
+iusin sin p, Sin T, cos 0, + ivsin sin Sin T, sin o,)" 
ist in ihrer Abhängigkeit von den Variablen z, y, z, u, v, w 
ein Integral der sechsdimensionalen Laplacegleichung: 


Dies ist unschwer nachzurechnen, so daß wir auf eine Aus- 
führung des Beweises dieser Behauptung verzichten können. 
Denken wir uns eine Funktion f(z) durch eine konvergente 
Potenzreihe ihres Argumentes gegeben, so ersieht man leicht, 
daß auch 


+ ty sin COs Py 
(49) +izsin sin p, cost, + tu sin sin Sin T, COS a, 
+ ivsin +, sin p, sin 7, sin o,) 


der Differentialgleichung (48) geniigt. Daraus finden wir dann 
das allgemeine Integral von (48), indem wir mit einer beliebigen 
Funktion der Variablen 9,, 99), T,, %) multiplizieren und dann 
über diese integrieren. Die Grenzen der Integrationen sind so 
zu wählen, daß das Integral wirklich der Gleichung (48) ge- 
nügt, was sich immer erreichen läßt. Ist z. B. die Funktion, 
mit der man (49) multipliziert, in bezug auf @,, 7,, 0, peri- 
odisch mit der Periode 2% und so beschaffen, daß sie für 
= 0 und +, = einen von den Variablen 9,, 0, un- 
abhängigen Wert annimmt, so genügt folgendes Integral der 
Gleichung (48), wenn D(3,,,,7),0,) eine Funktion mitdiesen 


Kr 
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(50) f f f [rw + ixcos 9, + iysin 4,0089, +++) 
| Po, II, dp, dr,do,. 


Die Forderung, daB die Integrale von (48) die Variable w nur 
als Faktor e*” enthalten sollen, führt die Gleichung (48) in die 
fünfdimensionale Schwingungsgleichung (46) über, während sich 
aus (50) das allgemeine Integral dieser ergibt: 

a 222222 

Wie wir zeigten, kénnen aber nur die Integrale Bezug 

auf unser Problem haben, die die Variablen z, u, v lediglich 
in der Verbindung 9 = Vz? + u® + v? enthalten. Aus (50) leiten 
wir solche Ausdrücke in folgender Weise her: Zunächst 
setzen wir: 
D(H), Pos Ts = Py) sin? 9, sin? sin , 
so daß also p,, unabhängig von o, wird. Ferner 
schreiben wir: 

z=0C08T; v=osinrsine, 
und definieren schließlich: 


cos = costcost, + sintsint,cos(o—o,). 
Dies in (51) eingesetzt, erhalten wir: 1 
a 2x 22 2x 


f Dy + y sin Dy c08 yo + sin 9, sin qo cos y) 
3 2 
:2(9,,9,)sin? J, sin? gp, sin t, d dg, dt, do. 


Beachtet man ferner die von Heine!) abgeleitete Gleichung: 


0 


unter Vernachlässigung Konstanten: 


2x 


(la W, (* o sin F, sin 


2(9,, Po) sin? F, sin? dF, 


‚Heine, 2.2.0. I, 8. 310. 
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Hier ersetzen wir o durch z und rechnen dann die Größen mt 2 
z, y, z auf die Koordinaten «, 9, [Gleichungen (7) und (8)] 
um. Da auch die vorhin eingeführte Funktion Y(a)- 0 (9, 
ein Integral der Schwingungsgleichung (46) war, so liegt es ne: 
nahe, die Darstellung a 


(imksna sin sin 9, sin p sin | 


a Po) sin? F, dd, dg; wo=x-g 


zu vermuten. Setzen wir hier a = 0, so werden wir auf folgende 
Darstellung fiir die Integrale der Differentialgleichung (47) 
geführt: 
(51b) O(%, = const f f (08 cos + k sin sin dy cos cos 

+ Q(H,, po) sin® sin? p, dF, d 
Hierin ist dann noch die zunächst willkürliche Funktion 2, 
von der wir nur Endlichkeit, Eindeutigkeit und Stetigkeit für 
0< t, <2; 0< gy, < 2a fordern, so zu bestimmen, daß die 
durch (51b) dargestellte Funktion © (9, 9) der Differential- 
gleichung (47) genügt. 


Wir führen einen Operator D #,9 ein durch die Gleichung 


= OF (sin 955) + sinp dp (sin 
— w* (cos? + k? sin? cos? g). 


Man zeigt nun, daß die Funktion ei® (cos # cos 4) +k sin A sin 9, cos p cos) a 


der Differentialgleichung 
Do,of —Daaf =9 

genügt. Dazu beachte man, daß der jetzt eingeführte Operator ee 


mit dem früher auf S. 627 eingeführten Operator in folgender —_— 
Weise 


Fi 


da ferner ei® (cos # cos + ksin sin dy cos 9cos 9) der Gleichung 
Ds,of — Do of =9 


Genüge leistet, ist es nur noch nötig zu zeigen, daß 


ir 
16 
sh 
te 
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h 
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st 
ER 
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FR 


‘ cot 0 cot @, 0 


ist; dies aber rechnet man leicht mit Hilfe der auf S. 628 ge- 
gebenen Formeln nach. 

Nun verfährt man genau wie früher, indem man mit Hilfe 
von partiellen Integrationen zeigt, daß die durch (51b) dar- 
5 gestellte Funktion sicher der Differentialgleichung (47) genügt, 
wenn diese selbe Differentialgleichung befriedigt. 
Man schließt daraus, daß die Lösungen der homogenen In- 
tegralgleichung 


(52) it, y) 1 f £08 p COS 


Po) Sin? JF, sin? gp, dt, dy, 


Integrale der Differentialgleichung (47) sind. 
4 Uns interessieren aber nicht die Integrale von (47), sondern 
wir suchen die Integrale der Differentialgleichung (13). Wir 
hatten aber gezeigt, daß die Integrale von (47), wenn man sie 
mit sin # sin g multipliziert, Integrale von (13) sein müssen. 
_ Wir werden daher die Gleichung (52) mit sin + sin  multipli- 
zieren und erhalten dann, wenn wir mit A(9, g) die Funktion 
x sin 9 sin y - O(, g) bezeichnen, den Satz, daB gewisse In- 


befriedigen : 
a 2x 
A(t, = af f sin # sin sin sin Py 
vv 

i (cos d,+ksin Hsin cos pcos po) . A Po) sin d ty 


- gleichung genügen müssen, so erhalten wir zur Antwort: Die- 
__ jenigen Integrale, die, durch sin 9 sin p dividiert immer noch 
if auf einer Kugelfläche endlich, eindeutig und stetig bleiben. 
Er Nun sieht man aber sofort, daß diese oe gerade 


er! 


Ar 
(58 


un 


(58 


De 
r 
(56 
| (57 
” Fragen wir aber, welche Integrale von (13) dieser Integral- 
> 
Be. Es ist leicht zu verifizieren, daß der Kern der Integral- u 
= gleichung (53) aus dem der Integralgleichung (19) dadurch 
ss hervorgeht, daß man diesen nach gy, differentiiert und dann 
mit — multipliziert: 
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sin sin sin p sin Po ei @ (cos # cos 9, + k sin sin cos cos gy) 


1 0 


(ei® (cos # cos du + k sin sin cos COs qo) ) 


Dadurch, daß wir diese Operation an der Reihe (28) ausführen, 
erhalten wir die Reihenentwicklung für 


sin sin sin p sin Po ei (cos # cos + k sin sin cos C03 @,) : 


sin sin sin p sin Po* eiw (cos cos d, +k sin sin 9, cos ~ Cos qo) 


mo, 1 wn(w Veos? + k* sin? cos? 9) 
(54) (Veos? + 1? sin? cos? 9)" 


{te m+ (cost, ksin cosy) P” (cos 9,)sinm Fo} 
m=1 


Fiir die noch fehlenden beiden Typen machen wir die 
Ansätze: 


55) U, (4,9) = (> Bs" (cos 9) sin 2my) 
n=1 m=1 
58) 9) => (> Pr, (cos 9)sin 2m 
n=1 m=1 


und zeigen wie früher, daß 


B,, = (—1jt!(4n + 1) 
a 2n 
Wen (o V cos? 9, + k? sin? +, cos? ,) 
67) J (Veos? $, + sin? 9, cos? Fo 
ı B™ = 4 8)- A-2m 
2n+1 ak 4 
f Won41(@ Vcos® 9, sin® 9, cos® qo) 
(68) J J 9, + sin? 
| 
, (cos 9,; ksin J, cos gp) 


»+1 (0, Po) Sin A, dd, dg, 
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F.Miglich 


| _ ist. Die Richtigkeit der Ansätze (55) und (56) folgt dann noch 
daß die Integrale 


. G” (cos ksin COS U, » Po) sin ty dit, dg, 


az 2x 


f wn(wV cos? 9, + k?sin? 9, cos? 
(Veos? + k* sin® 3, 


Po 


0 0 


a 2x 
J f Wn (@ Veos* 9, + sin? F, cos* 


(Veos? 9, + sin? 3, cos Po) 


@,, (cos 9, ; Asin Po) U, Po)8in dg, 


verschwinden, wenn in (57a) n=2o+1 und m= 2r oder 
m=2r+1 bzw. n=2o und m=2r+1 ist, oder ent. 
sprechend in (58a) n=2o und m=2r oder m=2r+1 
bzw. n=20 +1 und m=2r +1 ist. Der Beweis für diese 
Behauptung ist ganz analog dem früher geführten, so daß er 
hier unterdrückt werden kann. 

Schließlich ist auch in ganz der gleichen Weise, wie 
früher, der Konvergenzbeweis für die Reihen (55) und (56) 
durchzuführen. Auch dies soll hier unterlassen werden. Es 
ist zu diesem Zwecke nur zu bemerken, daß das Produkt 


(59) tg G™ (cos 9,; ksin 


ar 


wenn m > 0) ist (und nur solche treten in den hier betrachteten 
Reihen auf), endlich bleibt, weil fir gy, = = die Größe 
G™ (cos #,; ksin ı#,c0sY,) mindestens wie cos, verschwindet 
 [vgl. die Gleichungen (27). Daraus folgt dann noch, daß (59) 
im Bereich der Variablen 0=9,=r; ein 
Maximum besitzt, daß niemals größer ist als das von = 


(cos #,; ksin 9, cos y,), 


d.h. als Abschätzung für (59) können wir die auf S. 643 ge- 
gebene, für @” (cos #,; Asin 9, cos p,) gültige benutzen. 

Wir haben in den beiden letzten Paragraphen den Beweis 
erbracht, daß es acht und nur acht voneinander verschiedene 
Typen für die Integrale der Differentialgleichung (13) geben 
kann. Wir worden uns einer Methode zuwenden, 
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die uns in die Lage versetzt, für die Integrale von (13) an- 
genäherte Ausdrücke zu geben, werden aber, um nachher prin- 
zipielle Ausführungen nicht durch Nebenrechnungen unter- ae 
brechen zu müssen, einen Paragraphen einschieben, in dem ~ 
einige Formeln entwickelt werden sollen, die dann zur nume- aS; Er 
rischen Auswertung unerläßlich sind. 


Nach Heine!) ist es möglich, die Differentialquotienten 
der Kugelfunktionen als Entwicklungskoeffizienten einer be- 


Paty @) 
quotienten ha die Entwicklungskoefüzienten folgender 
Funktion: 
Vı-2rz+ rt = dx” 
Nach der Definition von 7 ist sicher i) haus 


und daraus folgt durch itera nach r: 


Hier ersetzen wir nun 7 durch die iaiaidiisiaidiie ale 
ordnen nach steigenden Potenzen von r. Die Gleichung kann — 
nur erfüllt sein, wenn die Koeffizienten der Potenzreihe ver- 
schwinden. Dies gibt eine Rekursionsformel für die Funk- — 
eP.« 
tionen ie Wir ersetzen dann noch x durch n—» 


v 


und beachten, daß der gleichen Rekursionsformel auch die 


Funktionen 
da” 
geniigen miissen. 2 finden: 


(2n + 1)-2-P" (2) + 1) PL, (2) +¢ 


2v + 1)-x-P” (xz) = P” , (2), 2 

1) E. Heine, a.a. O. I, S. 451 ff. 


Annalen der Physik. IV. Folge, 88, 
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Beugı 
Durch nochmalige Anwendung dieser Formeln ergibt 
sich weiter: Ri 
icht 

+ ) P’ (2) n>v+i1 Aus | 
4n?—1° n-2 
mit 

Mit Hilfe dieser Gleichungen gelingt es dann, die Integrale (62) 
PL 
2x 
cos? ‘+ P* (cos 9) (cos cos jpsin dd ddyp 
Von allen méglichen Größen ath, a die durch Variation lie | 
die 
des zweiten unteren und des zweiten oberen Index entstehen 
(wobei allerdings immer j= m sein muß), sind nur die von Null 
verschieden, die erstens gleiche obere Indizes haben, und für 
deren untere Indizes eine der Gleichungen m = n, n — 2, n+% 
erfüllt ist. Wir erhalten’): 
+ 1)Qn+ 8)(2n+5)’ zu | 
(2% —3)(2m—1) (2m+1) (n—r—2)! sit, 
(29 +5) (2¥+7)’ % +1; 41 3) (2v+5)’ Ind 
| 
1) Unter der Gleichungen ‘ 
ant 


+ i)! 
ve ind gleic ; —, wenn * = u 
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In Heines Handbuch der Kugelfunktionen') wird die 
Richtigkeit folgender Gleichung bewiesen: 
_ 

dx dz = (2n + 1) P (2). 
Aus ihr folgt durch »-fache Differentiation und Multiplikation 
mit yi — wenn die Ungleichung n>» + 1 erfüllt ist: 
(61) (2m + 1) V1 —2? PP (x) = (2) — Pi}, (a). 
Ferner folgt, wenn x >» + 8 ist, die Gleichung: a 

vr 1 +2 
[| 


(62) 


— 


2 vy +2 
| 


Setzt man schließlich fest, man wolle unter P’ (x) dann die 
Null verstehen, wenn der obere Index der Funktion größer 
als der untere ist, so gilt die Gleichung (62) ohne die durch 
die Ungleichung n > v + 3 eingeführte Beschränkung. 
Jetzt sind wir in der Lage, die Integrale 


2r 
o” j 2 DY 
J sin? cos? p P” (cos cos 


- P! (cos #) cos j p sin 
zu berechnen. Dazu beachte man, daß die Gleichung 7 
cos’ p cos Cos j y = = {cos(v + + cos(v — 1) gq} 

-{cos(j + 1)p + cos(j — 1)y} 
gilt, woraus geschlossen werden kann, daß die Größen o “aa 
mur dann von Null verschieden sind, wenn für die oberen 
Indizes eine der Gleichungen j =v, v +2; v— 2 erfüllt ist. 


Indem man die Integration über y ausführt, gelangt man zu 
Integralen: 


(2) (1 — 2%) dz; 2) (2)(1 — x) dz; 
«t 


x: 
‚+2 
Ly 
i 
= 


Diese lassen sich mit Hilfe der Formeln (61) a: Pe: ant 
die Form 


a 


bringen, wo beide Kugelfunktionen die gleichen oberen Indizes 
haben. Solche Integrale wissen wir aber gemäß der Fußnote 
auf S. 658 zu berechnen. 


Wenn wir noch als Definition einführen: wir wollen 


eS A unter el die Null verstehen, sobald einer der beiden oberen 
ea Indizes (oder auch beide) größer wird als der unter ihm 
en stehende, so können wir die gesuchten Größen wm: ‚ soweit 
nt x eine der Gleichungen j =», » + 2, v— 2 erfüllt ist, nur unter 


Beachtung dieser beiden Beschränkungen allgemein durch 
folgende Ausdriicke geben: 
n (n + »)! 


(4n? — 1)(2n — — »)! 


(n + »)! 
(n + »)! 
nn-2 2 (4n®— 1) (Qn— —»—4)!? 
(n +»)! 
(4n® ni’ 
o”” =o (m+r+2)! 
mn  (2n+1)? + 3)(m — »)! 


nem |; 
(2n — 1)(n — » + 2)! 


(n +» + 2)! 
(an? — 1)@n + y— 2)! 
22 n+ 1)(2n + 3) (22 + 5) (rn — »! 
2)! 
(n +» + 2)! 


+ 3)? (n — »)! 


n+r+4)! 
nn42 2 


(6 


in 


= 
gi 
“. 
Fu 
Ex 
die 
(6: 
au 
lä 
du 
d 
“| 
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§ 10. Das Variationsproblem, das der Differentialgleichung (13) 
aquivalent ist 


Wir betrachten das Integral Ir = 


+ (fle, 08+ [2° + (1 — 2%) gl} dx dy 


und suchen die in diesem vorkommende, zunächst willkürliche 
Funktion f(z, p) so zu bestimmen, daß das Integral (64) einen 
Extremwert annimmt, indem wir gleichzeitig der Funktion f(x, p 
die Nebenbedingung 


fife (2, g))?dzdg = | 


auferlegen. Dieses Variationsproblem mit Nebenbedingung oo 
läßt sich bekanntlich durch ein gewöhnliches ersetzen, namlich 
durch das 


(66) + (f (2 @) )? a? (x? + k*(1 — x?) cos? 


dg=0, 


in dem #’ ein noch zu bestimmender Parameter ist. 


daß folgendem Ausdruck gleich wird: 
22 
a 
0 
-f faz. dy-öf(e, (a — x’) = + 
“10 


— f (2, y) (@?[2? + (1 — 2*)cos* + F’)| 


so daB wir für das Verschwinden von JJ zwei Bedingungen 
finden, die gleichzeitig erfüllt sein müssen: Erstens hat f(z, ¢) 
der Differentialgleichung 


In bekannter Weise zeigt man durch partielle Integration, er a 


268 

en 
1m 
eit ( ) aa 
ter 
ch 

4 
) 

a P +4 


(67) 

— f (@? (2? + R?(1 — 22)cos?y] + E’) = 
zu genügen. Dies ist aber nichts anderes als die Differential. 
gleichung (13), in der nur an die Stelle von cos #* die Größe x 
gesetzt wurde. Die zweite Bedingung ist die, daß die ein- 
Integrale 


+1 
» Of d: 
ant verschwinden. Dieser Bedingung wird aber Geniige geleistet, 
wenn die Funktion f(z, gy) endlich, eindeutig und stetig für 
;  —1l=zzr=+1 und periodisch in g von der Periode 2x ist, 
Daher ist die Gleichung (64) mit dem in dieser Arbeit an die 


Zn u Wir werden für f(x,y) einen der Ansätze machen, 
deren Existenz wir in den 88 7 und 8 erwiesen haben. Da 
Sf der hier einzuschlagende Weg wesentlich unabhängig davon 
ist, welcher der acht verschiedenen Typen der Betrachtung 
zugrunde gelegt wird, wählen wir zur Veranschaulichung der 
folgenden Methode die durch die Gleichung (31) dargestellte 
Funktion, der als Typ der Bezeichnung AR,, (9, zukommt, 
_ Wir werden jedoch eine kleine Abänderung vornehmen, indem 


wir schreiben: 


v=0 
re Hier haben wir an die Stelle der Grenze oo für die über r 
zu erstreckende Summe die Zahl N gesetzt. Über N wollen 
‘ wir nichts voraussetzen; zunächst sei N nur eine sehr große 
Ex Zahl, die nicht weiter beschränkt ist. Wir entgehen durch 
diese Festsetzung der Notwendigkeit von Konvergenzbetrach- 
tungen für eine im folgenden aufzustellende Determinante. 
Wir können hinterher zeigen, daß die erhaltenen Ergebnisse 
von der Größe der Zahl N vollkommen unabhängig sind, woraus 
dann hervorgeht, daß N über alle Grenzen wachsen darf, ohne 
daß dadurch die angestellten Überlegungen hinfällig werden. 
a Setzen wir jetzt (68) in (64) rye 80 wird ov Integral I 


-* 


(cos cos 279). a 
=0 
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zienten ee; Da (64) einen Extremwert haben soll, haben 
wir auf die Funktion /(A,°, A,°, 4,”, 4,°...) die Regeln der 
Differentialrechnung anzuwenden, also die Gleichungen 


01 


69 2 


=0 
zu bilden. Die Anzahl dieser Gleichungen ist gerade der An- a 
zahl aller vorhandenen A?” gleich, so daß die 4?” aus diesen hear 
Gleichungen sich berechnen lassen werden. Die Gleichungen (70) 5 

lassen sich in folgender Form darstellen: 


} 
c+1 v=j=0 


Für die Größen ergeben sich die Gleichungen 


N m 
f 25 2i of 2% 25 2.2%, 2j 

7 2m + @ [ 2m + h os, 2m 

(70) ~ 

wo § 

2a +1 

9 

0 

0 wenn tm » 

(71) 

(27 + 2m)! : 
wenn 
=) 4r+1 (2r- 2»)! 
4n = m 


= 


2a +1 2j 
(72) 0 
| Pal a) sin g sin 2g} drdg, 


2y, 2j 2», 2j 


während die x," und o,”,, durch die Gleichungen (60), 
(60a), (60 b), bzw. durch (63) gegeben sind. 

Indem wir an den Gleichungen (72) die Integration nach » 
ausführen, sehen wir, daß nur diejenigen Koeffizienten von 
Null verschieden sein können, deren obere Indizes einander 
gleich sind. Wir erhalten so die BE .' von Null 
verschiedenen 0; 


j 
- 
| 
ze 


Möglich 


i Das erste der beiden Integrale integrieren wir partiell und 
erhalten: +1 9 


ER = — (x) dz (a 2) dz. 
 Beachtet man nun, daß die Zugeordneten der Kugelfunktionen 
der Differentialgleichung 
d 4v? 72V 
EM. genügen, so sieht man leicht, daß sich die weg von den 
0°”? nur durch den Faktor 2r(2r + 1) unterscheiden '): 


2,22 
2y, 2», 2» 
(74) = 2r(2r+ 1) 237 
Bezeichnen wir ferner mit a2”) die Größen + %20°” 
27, 2m 2r, 2m 2r,2m? 


so führt die Lösung der linearen Gleichungen (71) auf die 
_  Determinante, die auf der folgenden Seite dargestellt ist. 
In dieser Determinante ist jede Zeile und jede Kolonne durch 
ein Indexduo gekennzeichnet. Die Ordnung der Reihen und 
_ Kolonnen in (75) ist so vorgenommen, daß alle Zeilen bzw. 

Kolonnen mit dem unteren Index 2r vor den Zeilen bzw. 
_ Kolonnen angeordnet sind, deren unterer Index größer ist als 
2x. Ebenso sind Zeilen bzw. Kolonnen, die gleiche untere 
Indizes haben, nach der Größe der oberen Indizes geordnet. 
Die Koeffizienten a sind symmetrisch insofern, als man 
das erste Indexduo mit dem zweiten vertauschen darf, ohne daß die 
Koeffizienten ihren Wert andern. Es ist also =(, 


woraus folgt, daß (75) eine symmetrische Determinante ist. 


1) Die 637° 3%, mit verschiedenen unteren Indizes verschwinden 


wegen der Orthogonalität der P*(x) mit gleichem oberen Index. Vgl. 
Fußnote 658. 


+1 2» 2» 
> dP dP 
2y _ 2r 2m 
(i — 2") dx dz dz 
f 1 — 
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Das Gleichungssystem (69) ist 
linear und homogen. Daher gibt es 
dann nur von Null verschiedene Werte 
der m, die die Gleichungen (71) 
lösen, wenn die Determinante (75) 
verschwindet. Die einzige in (75) 
vorkommende Größe, deren Wert 


zahlenmäßig noch nicht eindeutig + j 
festgelegt ist, ist die Größe Z. Da- « Se 
her haben wir # so zu bestimmen, 
daß die Determinante verschwindet. 23 22% 
Man sieht sofort, daß die Gleichung, 4 
durch die Z’ bestimmt wird, eineso- ° © ° = BEE 
genannte Säkulargleichung ist, deren on 
Wurzeln alle reell sind, dadie Deter- or 
minante (75) symmetrisch ist. 
Die Determinante (75) ist eine > 
Funktion von Die durch das “S “s" 
Verschwinden der Determinante bee & 
stimmten Werte von E’ werden da- 
her im allgemeinen auch Funktionen 
von w?; es liegt nahe, für Z’ einen 
Reihenansatz mit unbestimmten Koef- 
fizienten zu machen, indem wir .. .. 
ca 
Auch die Determinante wird eine _- 9% 
Reihenentwicklung nach Potenzen von 


w* gestatten, um deren Konvergenz 
wir uns nicht zu kümmern brauchen, } a 
da die Determinante gemäß den über 
die Zahl N gemachten Voraussetzungen 
endlich ist; aus der speziellen Form der 
Abhängigkeit der einzelnen Glieder der 
Determinante (75) von w? folgt dann, 
daB auch die Potenzreihe endlich 
sein muß. In bekannter Weise ist: 


. . . . 
. . . . 8 ma 
oo on aa Ow + 
a a an = q 
+ 
a 
| 
3 
’ 
D 
ds 


F. Möglich 


d(09)® 


wo Q eine sehr große Zahl a die wir hier nicht näher be- 
rechnen wollen, die sich aber in eindeutiger Weise aus der 
Zahl N herleiten läßt, ebenso wie N eindeutig aus Q berechnet 
werden kann. 

5 Die Determinante (75) soll für in gewissen Grenzen be- 
liebige Werte von w? verschwinden. Dies ist nur dann möglich, 


wenn in der Reihenentwicklung (77) jeder Koeffizient für sich 


_ verschwindet. Wir bilden daher zunächst (D),,-o und suchen 
die Nullstellen dieser Funktion. Aus (75) folgt sofort, daß 
alle Elemente, mit Ausnahme der in der Hauptdiagonale, ver- 
schwinden, wenn wir » gleich Null setzen. Man kann diese 


Determinante dann in Form eines Produktes schreiben: 


[78 


(v=r=0 t=0 

R Aus dieser ersieht man dann sofort, daß der von w? freie 
 Koeffizient in (76) den Wert 


(78a) r=0,1,2,3...N 


hat. Auch sieht man, daß diese Nullstelle (z + 1)-fach ist. 
Es gibt also zt + 1 linear unabhängige Lösungen von (71), und 
damit mindestens r+ 1 linear unabhängige Integrale der 
Differentialgleichung (13) vom Typus R,,(9,g), wofern der 
Koeffizient &, von (76) gemäß (78a) gewählt wird. 

Das jetzt erhaltene Ergebnis setzt uns in den Stand, 
nähere Festsetzungen zu trefien über die Indizes, die wir den 
Funktionen gy), S(t, gq), T(t), geben wollen. 
Bisher haben wir nur zwischen geraden und ungeraden Indizes 
unterschieden, ohne nähere Eigenschaften unserer Funktionen 
angeben zu können, mit deren Hilfe eine Klassifikation möglich 
wäre. Jetzt aber können wir zeigen, daß die unserer Funk- 
tionen A(t, gy), die dem Werte ¢, = — 2r(2r + 1) entsprechen, 
für verschwindendes » in Kugelflächenfunktionen vom Grade 2r 
übergehen. Man sieht dies sofort, wenn man in der Differential- 
 gleichung (13) » = 0 setzt. Auch von Z’ bleibt dann nur das 
von w? freie Glied übrig, das, wie wir hier zeigen konnten, 
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gerade den Wert hat, der nötig ist, damit die Integrale von 
(13) Kugelflächenfunktionen vom Grade 27 werden. Wir setzen 
nun fest, daß unsere Funktionen 2, S, 7, U dann den Index m 
haben sollen, wenn sie sich für w = 0 auf Kugelflächenfunk- 
tionen vom Grade m reduzieren. Die Mehrdeutigkeit dieser 
Festsetzung werden wir durch Einführung eines zweiten Index 
erst etwas später beseitigen. 

Wir setzen jetzt den für &, gefundenen Wert in die Deter- 
mivante D [Gleichung (75)] ein. Dies hat zur Folge, daß in 
den Hauptdiagonalgliedern, die den unteren Index 27 haben, 
sich die konstanten, d. h. von w? freien Glieder wegheben. 
Es gibt somit t+ 1 Zeilen bzw. Kolonnen in der Determinante, 
in denen jedes Glied mindestens mit ? multipliziert ist. Daher 
können wir sicher die Potenz »°:+? vor die Determinante 
ziehen. Weil aber die Determinante D, wenn &, nach (78a) 
bestimmt ist, eine Potenzreihe nach wm? gestattet, deren 
konstantes Glied Null ist, so ist zu schließen, daß wenn 
e=— 2r(2r + 1) ist, alle | um 

6 
t+ 1 ist, verschwinden müssen. 


erst die Ableitung |-—_ . Diese betrachten wir zunächst. 
d(a)?** o=0 
TE 

Unter Dr» ™,™---™ wollen wir die Determinante ver- 
stehen, die aus D hervorgeht, wenn wir in D die Zeilen bzw. 
Sis Jay Jg Ju ( 
My, Mz, Mz. ++ My 
mehrere oder alle gleich sein können) nach w? differentiieren. 


Mit dieser Bezeichnung läßt sich dann der 7 + 1-te Differential- 

J2 +++ Jr 
quotient von D nach w? durch die Summe S Dmsm mr 
ausdrücken, wo über alle möglichen Indexduokombinationen 
zu summieren ist. Unter allen diesen Summanden ist aber 
im wesentlichen nur ein einziger, der für w = 0 nicht ver- 
schwindet. Dies sieht man so ein: 

Wenn eine oder auch mehrere derjenigen Kolonnen, die 
den unteren Index 27 haben, nicht ein einziges Mal diffe- 
rentiiert ist, so enthält sie, wofern ¢, nach (78a) bestimmt 
wurde, nur Elemente, die mindestens mit »? multipliziert 


sind. Setzen wir aber w gleich Null, so enthalten alle diese 


| , wofern eo kleiner als 
w=0 


Von Null verschieden ist 


Kolonnen mit den Indexduos unter denen auch 


FR 
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Möglich 


Determinanten mindestens eine Kolonne, in der alle Glieder 
Null sind. Diese Determinanten verschwinden also mit w. 
Alle anderen Determinanten, bei denen alle Kolonnen mit dem 
unteren Index 27 nach w? differentiiert sind, sind aber unter- 
einander gleich. Setzen wir = 0, so haben diese Deter- 
-minanten folgendes Aussehen: Die Elemente der Hauptdiago- 
ee die den unteren Index 27 nicht haben, sind von o frei, 
und eindeutig bestimmt, die anderen Elemente, deren beide 
er unteren Indizes nicht 2r sind, sind Null. Die Hauptdiagonal- 
elemente mit den unteren Indizes 2r sind die Größen 


a, +8 ae Die anderen Elemente, deren beide untere 


2z, ot 


‚Indizes 2r sind, ‘sind die GréBen a 
die fir ©=0 von Null verschieden zu dem Werte von 


pe Alle die Summanden, 


Diese Determinante entwickeln wir nach dem Laplaceschen 
_Entwicklungssatz, und zwar nach den Determinanten, die aus 
den Zeilen mit dem unteren Index 2r gebildet werden können. 
Man sieht sofort, daß aus den betrachteten Zeilen dieser 
a Determinante nur eine Determinante gebildet werden kann, 

die von Null verschieden ist. Diese besteht gerade aus den 
Elementen, deren beide untere Indizes gleich sind. Damit 
reduziert sich also die Bedingung, daß Fr ver- 

=0 

an sollen, auf die, daß die Gleichung 


a’ 0 a’ 2 


beitragen, haben also die Form: 
100000... 
- 
000001 0 0 0 00 
00 00 0 2 0o 4 0 2r 
000000 %9,,: 2, S25 2,22 34,22 997,22 
000000a5,9, 
000000 0 


(19a) 


| 
| 
| 
| 
Ar 
A= 
= 
> 
BEN 
| by 
2z,2r 22,28 


(19a) 


gelten soll. Dies ist eine Säkulargleichung vom Grade 7 + 1 
zur Berechnung von ¢,. 


Wir zeigen jetzt, daß die Wurzeln der Gleichung (79) alle 


voneinander verschieden sind. Daß die Wurzeln alle reell 
sind, folgt einfach aus der Tatsache, daß die Determinante (79) 
symmetrisch ist; außerdem hatten wir schon nachweisen 
können, daß die Größe Z’ als Wurzel der Gleichung (75) reell 
sein muß. 


Wir bemerken, daß die Elemente Er , aus denen sich 
die Determinante (79) aufbaut, nur dann von Null verschieden 
sein können, wenn die Differenz der oberen Indizes gleich 
Null oder 2 ist. Diese Behauptung folgt sofort aus den 


Gleichungen (60), (60a), (60b) und (63). Zur eat wollen 


wir die Glieder der Diagonalreihe von (79), a,’ arte mit 9, 
die Größen gg = Pr mit ec, und >: mit p, be- 
zeichnen. Dann läßt sich (79) in folgender Form schreiben: 

0; C15 0; 0; 

0; 0; 0; 0; 641 

0; 0; 0; Cr¢13 +P, 


Diejenige Determinante, die aus dieser (rz + 1)-reihigen 
durch Streichung der m letzten Zeilen und m letzten Kolonnen 
entsteht, wollen wir mit W,_»4+;, bezeichnen. Man erkennt 
dann leicht, daB, wenn W, nach seiner letzten Zeile oder 
Kolonne entwickelt wird, sich folgende Rekursionsformel ergibt: 
(80) W = + &) — W,-2; 

s=1,23,3..r+1. 
Wir setzen fest, W, sei 1, W_, sei Null. 

Die Funktionen W, sind in bezug auf e, rational und ganz 
vom Grade s. Wir betrachten jetzt die Nullstellen dieser 
Funktionen und leiten aus der Darstellung (80) drei einfache 
Eigenschaften dieser Funktionenfolge her: Erstens kann W, 
und W,_, keine Nullstelle gemeinsam haben. Denn gäbe es 
einen Wert von ¢,, für den sowohl W, als auch W,_ı ver- 
schwinden würde, so müßte auch an dieser Stelle W,_>, 
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zwis: 
damit dann W,.;, W,-ı ... W, verschwinden; dies aber ver- Erge 
stößt gegen die Voraussetzung W,=1. Zweitens, da (c,_»)? stell 
eine stets positive Größe ist, haben W, und W,_sa verschie- on 


denes Vorzeichen, wenn W,_, verschwindet. Drittens folgt daß 
daraus, daß die Funktionen W, ganz und rational vom Grade s 


alle 

in &, sind: Für sehr große negative Werte von «, nehmen die 
Funktionen W, positive oder negative Werte an, je nachdem 
s gerade oder ungerade ist. Für sehr große positive Werte Wu 
von &, sind alle Funktionen W, positiv. der 
Wir wollen annehmen, W,_s habe Nullstellen, die alle Fun 


voneinander verschieden sind; sie seien, der Größe nach ge- Wir 
ordnet, mit a, <w,<a,<<«@,_-a bezeichnet. Auch W,_, 
habe nur voneinander verschiedene Nullstellen, fiir die gelten 
soll: 8, < 8, < 2s << P,-ı. Schließlich sei angenommen, daß die 
zwischen je zwei Nullstellen von W,_, eine, und nur eine Null- Kup 


stelle von W,_s liege. Es sollen also die Ungleichungen nic} 
erfüllt sein. Unter diesen Voraussetzungen hat dann auch W, eo 
nur voneinander verschiedene Nullstellen, und zwischen je £ 
zwei Nullstellen von W, liegt nur eine Nullstelle von W,_ı. W 
Die Richtigkeit dieses Satzes ist so einzusehen: Für negative ie 
sehr große Werte von ¢, haben W, und W,_» gleiches Vor- ee 
zeichen. Für &, = ß,, also an der ersten Nullstelle von W, _1, 
sollen W, und W,_s verschiedene Vorzeichen haben, was nur 
möglich ist, wenn W, eine Nullstelle hat, die kleiner ist als der 
die kleinste Nullstelle von W,_ı. Zwischen der ersten und von 
der zweiten Nullstelle von W,_, geht W,_s voraussetzungs- ow 
gemäß durch Null, wechselt also sein Zeichen; daher muß die 
auch W, in diesem Intervall sein Zeichen einmal wechseln usw. ob 
Auf diese Weise sieht man sofort die Richtigkeit der eben W 
aufgestellten Behauptung. sti 
ae Jetzt beachten wir nur noch, daß die Voraussetzungen 26 
des Satzes für s= 3 erfüllt sind. Denn W, ist eine lineare ge 
Bu 9; u Funktion, die für große negative Werte von ¢, negativ ist, für 
große positive Werte jedoch positiv. An der Nullstelle von ha 
Sun! W, muß W, negativ sein, weil W, = 1 positiv ist. Da ferner Wi 
nae PR für sehr große negative wie positive Werte von &, die Funk- de 
eas tion W, positiv wird, muß es für W, zwei Nullstellen geben, 
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zwischen denen die eine Nullstelle von W, liegt. Aus diesem 
Ergebnis folgt zunächst, daß auch W, nur verschiedene Null- 
stellen hat, und daß die Nullstellen von W, zwischen denen 
von W, liegen. So schließt man weiter und findet schließlich, 
daß die Nullstellen von W,,,, d.h. der Säkulargleichung (79) 
alle verschieden sind. 

Die Wurzeln der Gleichung (79) können wir der Größe 
nach ordnen und bezeichnen die an m-ter Stelle stehende 
Wurzel mit &,”-1(z). Durch diese Festsetzung sind wir in 
der Lage, die früher unterlassene nähere Klassifizierung unserer 
Funktionen R(9,9), 8(9, 9), T(9,q), U(9,g) durchzuführen. 
Wir hatten bisher nur festgesetzt, unsere Funktionen sollen 
als Index den haben, der den Grad der Kugelfunktion angibt, 
auf die sie sich für m = 0 reduzieren. Die Unbestimmtheit, 
die in der Existenz einer gewissen Anzahl linear unabhängiger 
Kugelfunktionen von gleichem Grade liegt, wußten wir zunächst 
nicht zu beheben. Jetzt setzen wir fest, daß die Integrale 
von (13) einen zweiten Index erhalten, der durch ein Komma 
getrennt, neben den ersten gesetzt wird. Dieser Index soll 
übereinstimmen mit dem Index, den wir der Größe ¢,"-1(r), zu 
der das betreffende Integral gehört, geben. Mit anderen 
Worten: Wir suchen unter den r + 1-Wurzeln der Gleich. (79) 
eine aus, z. B. &,”1(r), nehmen ferner an, wir hätten das Integral, 
das unter Benutzung der Größe 


= —2r(2r + 1) + + whe, +... 


der Differentialgleichung (13) geniigt; dann haben wir dem 
vorliegenden Integral als ersten Index die Zahl (27), als 
zweiten Index die Zahl m — 1 zu erteilen; das Integral hat also 
die Bezeichnung m - y). Allerdings wissen wir noch nicht, 
ob zu einem bestimmten Wert von z und zu einer bestimmten 
Wurzel &,”-1(r) wirklich nur je ein Wert der noch nicht be- 
stimmten Größen &,, &, &,... usf. gehört. Erst wenn wir 
zeigen können, daß dies der Fall ist, ist die von uns vor- 
genommene Indexbezifferung eindeutig. 

Um den Beweis, dessen Notwendigkeit wir eben angedeutet 
haben, möglichst übersichtlich durchführen zu können, wollen 
wir eine neue Determinante einführen. Wir erhalten sie aus 
der Determinante (75), indem wir dort alle Zeilen und Kolonnen 
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F. Möglich 
streichen, deren unterer Index kleiner oder gleich 2r — 4 ist, 
ebenso alle Zeilen und Kolonnen, deren unterer Index größer 
oder gleich 2r +4 ist. (Es sei noch einmal darauf verwiesen, 
daß in der Determinante eine Zeile oder Kolonne erst durch 
Angabe zweier Indizes bestimmt ist.) In der hierdurch ent- 
stehenden Determinante differentiieren wir noch die Kolonven, 
deren unterer Index 2r ist, nach w*. Die so erhaltene (37 + 3} 
reihige Determinante bezeichen wir mit A und nennen sie die 
reduzierte Determinante. 

Wir hatten bisher gezeigt, daß aus der Bedingung D =0 
die Gleichungen 


d?D \ 
Do=0 = 9; 


folgten, von denen die Gleichungen (2. : =0; ost 
d wo=0 

dann von selbst erfüllt waren, wenn wir ,=— 2r(?2r+]) 


+1 
-0 führte auf 
@=0 


Säkulargleichung (79) und lieferte uns die Bestimmung von &.. 


i 


gefolgert werden kann. 

Die Gleichung = 0 läßt sich mit Hilfe des 
a=0 

+ 


5 schreiben als 


setzten. Die Gleichung | 


92 


d’+“D 
d 


von uns früher eingeführten Symbols D 
Summe solcher Determinanten: 
Is D 0, 

wo über alle möglichen Indexduokombinationen zu summieren 
ist. Vorweg betrachten wir alle diejenigen Determinanten der 
Summe, in denen mindestens eine Kolonne mit dem unteren 
Index 27 nicht differentiiert ist. Man sieht sofort, daß, wenn &, 
nach (78a) bestimmt ist, diese Determinanten mit w vel- 
schwinden. Unter den übrig bleibenden wollen wir zwei Fille 
unterscheiden: Außer den Kolonnen mit dem unteren Index 21 
ist noch eine Kolonne differentiiert, deren unterer Index kleiner 
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oder gleich 27 — 4 oder größer oder gleich 27 +4 ist. Wir 
wollen diesen Fall charakterisieren, indem wir sagen, die Diffe- 
rentiation sei außerhalb des reduzierten Schemas vorgenommen. 
Der zweite Fall ist dann der, in dem die Differentiation inner- 
halb des reduzierten Schemas vorgenommen ist. Jetzt beachte 
man, daß die Elemente a,” hd der Determinante nur dann von 
Null verschieden sind, wenn die Differenz ihrer unteren Indizes 
Null oder zwei ist. (Genau das gleiche Gesetz gilt auch fir 
die oberen Indizes.) Wenn nun eine Kolonne außerhalb des 
reduzierten Schemas differentiiert ist, so sind die Elemente, 
die diese Kolonne mit den Zeilen vom unteren Index 27 ge- 
meinsam hat, sicher Null. Wenn wir nun w = setzen, so 
sind die einzigen Elemente der Zeilen mit dem Index 2r, die 
von Null verschieden sind, gerade die Elemente der Säkular- 
gleichung (79. Man kann dann die Determinante nach La- 
place entwickeln und findet als Ergebnis das Produkt zweier 
Determinanten, von denen die eine die Säkulardeterminante (79) 
ist. Die Größe e, war aber als Wurzel dieser Gleichung be- 
stimmt, so daß also alle Determinanten, die zum ersten Fall 
gehören, verschwinden müssen. 

Alle anderen Determinanten aber, die zum zweiten Fall 
gehören, haben die Eigenschaft gemeinsam, daß alle Elemente 
außerhalb des reduzierten Schemas Null werden, mit Ausnahme 
der Diagonalelemente, die bestimmte, hier nicht näher inter- 
essierende Werte annehmen. Man kann daher die Determinanten __ 
sehr leicht entwickeln und findet außer einem konstanten 
Faktor, der das Produkt aller nicht zum reduzierten Schema 
gehörenden Diagonalkonstanten ist, den ersten Differential- 
quotienten nach w? der reduzierten Determinante 4 für w = 0: 


(61) = 


Die Gleichung (81) zerfällt in eine Summe von Determinanten, 
von denen jede aus 4 durch Differentiation einer Zeile ent- 
steht. Wir unterscheiden wieder zwei Fälle 1. Fall: Die 
Differentiation ist an einer der Kolonnen ausgeführt, die lee — 
mente der Säkulargleichung enthält. Die Elemente dieser 
Kolonnen sind konstant, mit Ausnahme der Diagonalglieder, 
die noch von w? abhängen. Die Differentiation macht also 

Annalen der Physik. IV. Folge. 83. 
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F. Möglich 
alle Elemente zu Null, mit Ausnahme der Diagonalglieder, 
Wenn man w=0 gesetzt hat, ist hier allein die Größe s, 
vorhanden. Damit ist gezeigt, daß es + + 1 Summanden gibt, 
die die Größe «, einmal als Faktor enthalten. Die Differen- 
tiation an Kolonnen, die keine Glieder der Säkulardeterminante 
enthalten, und dies ist der 2. Fall, ergeben Determinanten, in 
denen, wenn wir w = 0 setzen, die Größen «,; &,; &,... über- 
haupt nicht vorkommen. Wir haben somit zur Bestimmung 
von &, eine lineare inhomogene Gleichung, deren Koeffizienten, 
wie man leicht sieht, eindeutigeFunktionen von z und &,”-1(7) 
sind. Daher existiert zu jedem Wertepaar tr; &,”-!(r) ein und 


nur ein Wert von 
Wir betrachten jetzt die Gleichung 4 


d (w? r+3 w=0 

Diese werden wir wieder mit Hilfe des früher eingeführten 
Symbols als eine Summe von Determinanten schreiben: 


Is D"' Ne; =(. 
Man sieht wie vorher ein, daß die Determinanten, in denen 
eine Kolonne mit dem unteren Index 2r nicht differentiiert 
ist, für @ = 0 verschwinden muß. Unter den übrig bleibenden 
werden drei Fälle zu unterscheiden sein: 

1. Fall: Zwei Kolonnen außerhalb des reduzierten Schemas 
sind differentiiert, wobei auch der Fall, daß eine Kolonne 
außerhalb des reduzierten Schemas zweimal differentiiert ist, 
hierher gehört. 

2. Fall: Eine Kolonne außerhalb, eine Kolonne innerhalb 
des reduzierten Schemas ist differentiiert. 

3. Fall: Zwei Kolonnen (oder eine Kolonne zweimal) inner- 
halb des reduzierten Schemas sind differentiiert (wobei die ein- 
malige Differentiation jeder Kolonne mit dem unteren Index 2r 
in dieser Zählung als selbstverständlich nicht mitgerechnet ist). 
Die Determinanten, die zum ersten Fall gehören, sind 
für © = 0 selbst Null. Beweis wie im ersten Fall der eben 
erledigten Untersuchung. 

Unter den zum Falle 2 gehörigen Determinanten be- 
trachten wir alle die, bei denen die gleiche Kolonne außerhalb 


# Ber 
= des 
für 
des 
Nu 
ren 
das 
der 
4 auf 
ele 
der 
q sta 
Sch 
die 
ma 
set 
au 
au 
Fa 
vol 
all 
ti 
set 
a so 
® da 
2 sc] 
M: 
lie 
ei 
de 
nt 
m 


Beugungserscheinungen an Körpern v. ellipsoidischer Gestalt 


des reduzierten Schemas differentiiert ist. Wenn wir dann 
fir © = 0 Null setzen, gibt es in allen Kolonnen außerhalb 
des reduzierten Schemas nur die Diagonalelemente, die von 
Null verschieden sind, mit Ausnahme der Elemente der diffe- 
rentiierten Kolonne. Jedoch enthält dann diejenige Zeile, die 
dasselbe Indexduo wie die differentiierte Kolonne hat, außer 
dem Diagonalelement nur Nullen, wie auch alle übrigen Zeilen 
außerhalb des reduzierten Schemas mit Ausnahme der Diagonal- 
elemente nur Nullen enthalten. Daher können wir immer nach 
den Zeilen entwickeln und erhalten, abgesehen von einer Kon- 
stanten, die das Produkt aller außerhalb des reduzierten 
Schemas befindlichen Diagonalkonstanten ist, eine Determinante, 
die aus der reduzierten Determinante 4 dadurch entsteht, daß 
man in ihr eine Kolonne nach »? differentiiert und dann w = 0 
setzt. Alle Determinanten vom Falle 2, die die Differentiation 
außerhalb des reduzierten Schemas an der gleichen Kolonne 
aufweisen, haben, wie man leicht sieht, den gleichen konstanten 
Faktor vor der Determinante, die aus 4 durch Differentiation 
einer Zeile nach w? entsteht. Diesen Faktor kann man daher 
vor die Klammer ziehen und hat in der Klammer die Summe 
aller verschiedenen Determinanten, die aus A durch Differen- 


tiation einer Kolonne nach wm? entstehen, wenn man w = 0 
setzt. Dies ist aber nichts anderes als (sor) og 8 ist aber 


so bestimmt, daß diese Größe verschwindet; damit ist erwiesen, 
daß auch die zum Falle 2 gehörigen Determinanten ver- 
schwinden. 

Es bleiben somit nur noch die Determinanten vom Falle 3. 
Man zeigt leicht wie früher, daß nach Entwicklung nach den 
Zeilen oder Kolonnen, die außerhalb des reduzierten Schemas 
liegen, als Faktor nur übrig bleibt der zweite Differential- 
quotient nach w? der reduzierten Determinante 4. Dann zeigt 
man noch genau wie vorher, daß die Gleichung 


d?4 

eine lineare inhomogene Gleichung zur Bestimmung von ¢, ist, aus 
der dann wie früher folgt, daß zu jedem Wertepaar z, &,”-1(r) 


nur ein Wert von «, gehört. In ganz der gleichen Weise zeigt 
man dann auch, daß sich e,; (vy > 1) aus der Gleichung 
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676 j 
ss 0 
d(@*) a=0 
berechnen läßt, wenn &,; &; ¢...¢,, bereits bekannt sind, 
ferner auch, daß diese Gleichung linear und inhomogen ist 


(82) 


mit Koeffizienten, die eindeutige Funktionen des Zahlenpaares 


a daB der Gleichung 


tT, &”(r) sind, so daß man in voller Allgemeinheit sagen kann: 
Zu einem einzigen Wert von rt gehören tr + 1 und nurr+1 
Werte der Größe E’. Jede der Wurzeln der Determinante (75) 
ist einwertig; daher gehört zu jeder Wurzel EZ’ nur ein bis 
auf einen konstanten Faktor bestimmtes Integral. Den zu- 
nächst willkürlichen Faktor bestimmt man durch die Forderung, 


a0 dy=1 
56 


geniigt werden soll. 

Die im vorhergehenden ausgesprochenen Sätze über die 
Nullstellen der Determinante (75) sind von der ursprünglich 
eingeführten Größe N, die wir brauchten, um einem Konvergenz- 
beweis aus dem Wege zu gehen, vollkommen unabhängig; 
höchstens ist der Forderung zu genügen, daß die Größe N 
größer ist als die Zahl r, eine Bedingung, der immer genügt 
werden kann. Hieraus ist, wie schon vorher angedeutet, zu 
schließen, daß die Ergebnisse, die wir für die Wurzeln der 
Determinante (75) abgeleitet haben, exakt gültig sind. 

Wir wollen jetzt als Beispiel denjenigen Wert von L’ 
berechnen, der zu dem mit 2, (+, y) bezeichneten Integral 
gehört. Wir können in diesem Falle den zweiten Index weg- 
lassen, weil es nur eine einzige Kugelfunktion der Ordnung 
Null gibt. Wir haben daher in (78a) für z Null einzusetzen 
und finden ¢,=0. Die Säkulargleichung (79) reduziert sich 
in unserem Falle auf ae 


(1+ #) +422, =0, 


woraus dann für «, sich folgender Wert ergibt: @ 


(83) 


Wir schreiten jetzt zur Aufstellung der reduzierten Determinante, 
die in diesem Falle dreireihig ist und folgende Form hat: 
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4m 16a 


Aus ihr finden wir nach Abspaltung aller unwesentlichen kon- 
stanten Faktoren die Determinante: 


ave 


2h? , 

>? 6+ H+ +242) 


&; &5 ... Wir haben die Rechnung für e, und Gurchgfihr eh 


und finden: 
= 13 k4); 
4 + 108 — 317 4229). 
Es folgt somit als Wert von Z’, der zum Integral R, (+, p) 
gehört: 


100 + 108 k? — 317 kt — 42 


70 875 


Diesen Wert von Z’ hätte man jetzt in die Gleichungen (71) = 
einzuführen und dann dieses System aufzulösen. Dabei würde as os 

man bequemlichkeitshalber dieses System inhomogen machen, 
indem man willkürlich 4,° gleich 1 setzt. Dadurch verstößt 
man zwar gegen die Bedingung (67), die jedoch hinterher leicht 
durch Hinzufügung eines geeigneten Faktors wieder erfüllt 
werden kann. Die approximative Auflösung des Gleichungs- 
systems (71) bietet theoretisch nicht die geringsten Schwierig- 
keiten, doch dürfte die numerische Auswertung zwar nicht 
unmöglich, aber doch recht kompliziert sein. Daher wollen 
wir hier darauf verzichten. 
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Möglich 


Aus der Tatsache, daß die Funktionen R(t, y); 8(9,9); 
U(9,g) die Lösenden eines bestimmten Randwert- 


oder v+j 


Dn, sin i} ditdg 


v0 
f sin dF dp 1 wenınn=m 
) | und 
a 2x 
SSR... Se, sin #49 49 = [ ‚sin = 0 
00 00 
a 2x | 


SS ‚sin 9d9dp= [ [8,,7, ‚sin =0 
0 0 
22 a 2x 


«Sin = [ [T,,U, dodgy =0. 
0 0 0 


Aus diesen erkennt man, daß die hier gefundenen Funktionen 
sich in gleicher Weise zur Vornahme von Reihenentwicklungen 
eignen, wie z. B. die Kugelflächenfunktionen. Hierauf werden 
wir in einem der folgenden Paragraphen noch zurückkommen. 


§ 11. Die Integrale der Differentialgleichung (11) 


Wir zeigen jetzt, daß die Integration der Differential- 

_ gleichung (13) und die der Differentialgleichung (11) äquiva- 
em lente Probleme sind, insofern als die Kenntnis der Integrale 
(13) die der Integrale von (11) nach sich zieht und um- 
ae a gekehrt. Diese Behauptung folgt aus der Tatsache, daß es zu 
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einen Wert von H in der Differentialgleichung von (11) geben _ 
kann, zu dem ein einziges Integral mit periodischen Eigen- 
schaften gehört. Wäre dem nicht so, gäbe es also zu einer — 
bestimmten Wurzel der Determinante 75 mehrere Werte von H, 
za denen mindestens je ein Integral mit periodischen Eigen- 
schaften gehört (oder auch nur ein Wert von H mit zwei linear 
unabhängigen Integralen der gleichen Eigenschaft) — diese seien _ 
mit fi; fa---/„ bezeichnet —, so sind, wie wir früher zeigten, 
die Produkte (8)- fly); f(x) Integrale 
von (12) oder, mit Hilfe von Gleichung (8) auf % und gg um- 
gerechnet, Integrale von (13), die sicher linear unabhängig, zu — 
dem gleichen Wert von E’ gehören. Im vorigen Paragraphen — 
zeigten wir aber, daß es zu jeder Wurzel von (75) nur ein 
Integral mit Kugelfunktionseigenschaften gibt. Damit ist der 
eben behauptete Satz bereits bewiesen. Wir gelangen also zu 
einem Integral der Differentialgleichung 


—(wrdnty + Ban’y + H)f =0, 


wenn wir in einem Integral von (13) 9 und m gemäß Glei- — 
chung (8) durch # und y ersetzen und dann für # einen kon- — 
stanten Wert A, substituieren. Nehmen wir z. B. das Integral Sean 


n,v 62 p) = (> (cos cos 2j 9) . 

o=0 j=0 

Hier setzen wir 8 =0, was dem Wert gm =O entspricht; = ; 
so kommen wir zu einem Integral der fraglichen Differential- My a 
gleichung, das folgendes Aussehen hat: ae 


=0 
Doch ist. die Reihenentwicklung nach Zugeordneten a 
Kugelfunktionen mit variablen unteren und oberen Indizes 
nicht eindeutig, so daß diese Reihen nur geringen Wert 
haben. 

Zu brauchbaren Entwicklungen aber kommen wir, wenn 
wir denselben Vorgang an einer der früher eingeführten Integral- 
gleichungen vollziehen, indem wir jetzt, da die Berechnung der BE Dee 
gesuchten Funktionen bereits durchgeführt ist, die Integral- me 
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gleichung als Integraldarstellung für unsere Funktionen be- 
trachten. Wir wollen z. B. den Kern (20) zugrunde legen. Indem 
wir dort cos # und sin I COs p nach (8) umrechnen und dann 
ß = K setzen, erhalten wir die Gleichung: 


a 2x 


(86) = ff a) SiN Fy dy d 
00 


Hier haben wir an die Stelle des Wertes A die Größe do,, 
geschrieben, um anzudeuten, daß es sich um den Eigenwert 
der Integralgleichung (19) handelt, der zu der speziellen Eigen- 
funktion R2,,,(%, gp) gehört. Wir ersetzen nun sowohl die 
e-Funktion als auch die Funktion R2,,,(%,,q,) durch ihre 
nach Kugelfunktionen fortschreitenden Reihen, von denen jede 
gleichmäßig konvergiert. Daher ist gliedweise Integration ge- 
stattet und liefert die nach den -Funktionen fortschreitende 
 Beibe, die wir mit c - Ra„,,(Y) bezeichnen wollen: 


| 
= c+ Ran,,(Y) = Wo o(adny). 


= exakten Definition der edt die durch die Reihe (87) 


dargestellt ist, ist es nötig, die Berechnung des Eigen- 
wertes Aa, nachzutragen. Zu diesem Zwecke denken wir uns 
in Gleichung (19) an die Stelle der Funktion F($, gy) die 
Eigenfunktion R>,„,($,g) eingetragen. Dann setzen wir 
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Die ia Gleichung (87) noch auftretende Konstante c ist, 
wie man sofort sieht, nichts anderes als Ro, ,(K), so daß wir 
als Definitionsgleichung für die Funktion 
Gleichung erhalten: 


=) Sn 
(88) Ron, (¥) = (— 
o=0 


Zu einer anderen Entwicklung gelangen wir, wenn wir ER mo 
in dem Kern (20) nicht 6 = K, sondern B= K+iK’' setzen. 2 
Auf diese Weise erhalten wir: ear 


hon,» iokenyeind con (Roy y(Fo; Po) 
Hier wir die Integration über gy, zunächst aus, 
sich dann folgende Reihe ergibt: a 


(- 4: | Jo; (— ho eny sin P}(cos 9,)sin 9, 


Be = 


Die Integration über ist mit Hilfe der Gleichungen (30) | 


und (27) leicht durchgeführt. Setzen wir zur Vania = as 


so ergibt sich die Reihe: ak a 
R, n,v (+ 
(90) Ron, = Ä 


=0 


2A} Ro, ,(K+iK’) 


Eine dritte Entwicklung ergibt sich mit Hilfe des Kernes 
K,(%,p; 9,,9,) aus (20a). Setzen wir hier K+ik’, 


so gelangen wir zu der Gleichung: 


x 
: 
a 
A 
ig 
0 
A, 
BR: 
23 
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z 22x 


= % Rony (Fo; p)sin t,dq,, 


WO fon,» ein des Kernes K, (:‘+,p; ist. Dieser 
läßt sich in gleicher Weise wie die Bigenwerte des Kernes 
K, (4, %,,9,) berechnen; wir finden: 


22x 


ergibt sich dann die 


ny = 


24° R,,, (K 


Hier ist unter Az, folgende Größe zu verstehen: 


o 


(20 + 27)! 2j 
A, = — 1)j 


Eine vollkommen analoge Reihenentwicklung wie Glei- 
chung (87) erhält man für die Funktionen 22,,41,,(y), nur daß an 
die Stelle der w-Funktionen mit geradem Index solche mit un- 
geradem treten und daß sich die Eigenwerte Aa„;ı,, ein wenig 
anders berechnen. Man differentiiere die Gleichung (19) ein- 


mal nach :# und setze dann ? = + p= + 


5 3 80 ergibt sich, 


wenn man für R2,41,,(:%); 9) die Reihenentwicklung einsetzt: 
f 5) a 2x 
a9 = fiw cos 
0 0 


* Rongs, Po) 8in dit, 
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Die zu den Reihen (90) und (91) analogen fiirdie Funktionen 
Rau+1,,(9, 9) sind nicht so leicht herzustellen. Wir benutzen 
wieder die Integralgleichung (19), ersetzen den Kern durch die — 
ihm äquivalente Reihe (28), setzen schließlich an dis Stelle 
von F(t, gp) die Funktion R2,,1,,(%, 9) bzw. die dieser ent- 
sprechende Reihe (33a). Schließlich rechnen wir noch cos + 
und sin cos nach (8) auf die Größen $ und 4 um. Soer-. 
halten wir nach Integration die Reihe: ; ’ 


See 


2o4+1 k’ 20+1 2 ts 


(20 +29 + 1)! dnßdny. 2j 
+> (20 —27+ 1)! k’ k’ A; 4 


Wir beachten nun, daß die Funktion Ron,1,,(%,g), wie 
man leicht aus der Reihe (33a) erkennt, durch cos # teilbar ist. Be 
Wegen (8) muß man schließen, daß daher R2,,1,,(8) durch 
teilbar ist und daher mit dxf verschwindet. Daher darf man En ds 
jetzt nicht mehr 6 = K+iX’ setzen, weil dann dnß Null vi 
werden wiirde, sondern wir haben vorher die Gleichung (93) noch 
durch dn zu dividieren. Dann erst kann man 8 = K+iK’ 
setzen, ohne auf beiden Seiten der Gleichung (93) die Null =| 
zu erhalten. Aus (27) rechnet man leicht nach, dab E. 


1 9; 1 _ th 


dny (cny)* 


(20 — 2j + 1)! 


ist. Schreiben wir dann noch 


(20 + 2j + 1)! 


j=0 


so erhalten wir folgende Reihe: 
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0 Ron41, 90) 
(94); 20 A 
1 d 


In ähnlicher Weise erhalten wir aus dem Kern 


od 


Rong, oY) = 


dn 
tke A2c+i 


sny 


20+ 27 +1)! 
= 1); 


bedeutet. Um diese Reihe zu verifizieren, beachte man, daß 
die Eigenwerte u>„;1ı,, des Kernes K,(#, p; 9,, 9,) sich so 
darstellen lassen: | 


(5; 0) 


Ot 
96 


Es würde zu weit führen, für jede der noch fehlenden 
sechs Typen alle möglichen nach den w-Funktionen fort- 
schreitenden Reihen anzugeben; wir wollen uns auf die hier 
angegebenen beschränken, da es nicht schwer ist, in genau 
der gleichen Weise die übrigen Reihenentwicklungen her- 
zustellen. 
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§ 12. Die Integrale der zweiten Art der Differential- 
gleichung (11) 

Durch die bisherigen Entwicklungen sind wir in die Lage 
versetzt worden, Integrale der Schwingungsgleichung (10) an- 
zugeben, die im Endlichen überall endlich sind. Es ist be- 
kannt, daß solche Integrale gerade die einfallenden Wellen 
darstellen. Zur Beherrschung der Beugungserscheinungen ge- 
hört aber auch die Kenntnis der von den beugenden Teilchen 
wegeilenden Wellen, deren analytische Ausdrücke wir noch 
nicht angegeben haben. Dieser Ableitung wenden wir uns 
jetzt zu. er 
Wir erinnern daran, daß das Integral (86), durch das die — 
Integrale der ersten Art dargestellt wurden, aus der Integral- 
gleichung (19) direkt folgte. Diese aber erhielten wir aus 
einem allgemeinen Integral der Wellengleichung, indem wir 
die Grenzen so wählten, daß gewisse bei der partielen 
Integration auftretende Integrale verschwanden. Dies Ver-- 
schwinden war die hinreichende Bedingung für die Richtig- 
keit unserer Integralgleichung. Sie führte auf die Grenzen 0 — 
und = für 9,, und 0 und 2r für y,. Man kann sich jedoch 
die Frage vorlegen, ob die so bestimmten Grenzen die einzig 
möglichen sind. Dies ist dann zu verneinen, wenn man auch 
komplexe Integrationswege zulassen will, und es gibt in der 
Tat keinen Grund, der eine Beschränkung auf nur reelle 
Integrationswege rechtfertigen würde. Doch wollen wir den 
Weg für die Integration über g, ebenso wählen wie früher, — 
also zwischen 0 und 22; nur der Integrationsweg für 9, soll 
geändert werden. ; 

Betrachten wir einmal die auf Seite 630 durchgeführte 
partielle Integration, die uns die Integralgleichung (19) lieferte, 
so erkennen wir, daß das eine der beiden einfachen Integrale, 
nämlich : 

toe 

f Rn» Po) — - = addy; 
= cos cos + ksin sin Cos COS 
verschwindet, weil der Integrand eine Differenz zweier Funktions- — 


werte für die Argumente 9, = 0 und gy, = 2x ist, die vor- 
kommenden Funktionen aber alle endlich und in bezug auf y, 


4 
; 
om 
5 
B 
0 
r 
1 


yon der Periode 2a sind. Wir haben die Grenzen für 9, nun 
so zu wählen, daß das andere Integral 

d iwe 

f | 35, Ra,» (Fo Po) — 


Po = ( 


OR, . d=b 
n, Po d 


8%, 

& = cos cos ‘+, + ksin sin cos Cos 
ebenfalls verschwindet. Rechnen wir cos # und sin # cos g 
auf die Koordinaten 8 und y um und setzen schließlich y = K, 
so erhalten wir nach Ausführung der Integration: 


a, 
[sin ty Rav (Fos Po) 
SER Dann sind die Stellen a und 5 so zu wählen, daß jeder 

eh der Summanden Null wird. Für den Punkt a wählen wir 


wieder Null, den Punkt 5 jedoch nehmen wir im Komplex- 
Unendlichen so an, daß der zweite Summand exponentiell 
N verschwindet. Mit dieser Wahl von 5 erhalten wir dann ein 

zweites Integral der Differentialgleichung (11) in der folgenden 
Form: 


097) (8) = const f 9, , (cos 9,)sin I, 
0 


wo die Funktion 9 durch ein Integral oder die diesem äquiva- 
lente Reihe definiert werden kann: 


On,» (cos 4) = 
0 


(98) 02n,,(608 4) = 2a >" P, (cos 9); 
2 | 
(C08 H) = 2a (Cos). 
o=0 


a Anstatt 5 so zu definieren wie eben, hätten wir auch 
fordern können, 5 sei im Komplex-Unendlicheu so bestimmt, 
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daß das Integral (97) konvergiert. Man sieht leicht, daß beide 
Forderungen auf das gleiche herauskommen, jedoch ist es 
zweckmäßiger, diese zweite Formulierung zn wählen. 

Doch auch mit dieser ‚läßt sich die gestellte Aufgabe 
nicht direkt lösen. Denn im Unendlichen ist nicht nur die 
e-Funktion, sondern auch die Funktion o,,, wesentlich sin- 
gulär. Man müßte also das asymptotische Verhalten der 
Funktion g,,, kennen, ehe man Exaktes über das Integral aus- 
sagen kann. Asymptotische Formeln für o,,, anzugeben ist 
dem Verfasser aber nicht gelungen. Aus diesem Grunde soll 
hier ein etwas anderes Verfahren eingeschlagen werden, 

Da die Reihe (98), wie man leicht zeigt, für alle end- 
lichen Werte des Argumentes gleichmäßig konvergiert, so 
können wir versuchen, diese Reihe in (97) einzusetzen und 
dann gliedweise zu integrieren. Tun wir dies aber, so ge- 
langen wir zu folgenden Integralen: 

b 
(99) Jf P,(cos 9,)sin d Dy. 
0 
Wir wollen nun einen Integrationsweg bestimmen, längs 
dessen die Integrale unabhängig vom Wert des Index der 
Kugelfunktionen konvergieren. Dies aber ist nicht mehr 
schwer, da die Kugelfunktionen als rationale Funktionen im 
Unendlichen nur einen Pol haben, dessen Ordnung der Ord- 
nungszahl der Kugelfunktion gleich ist. Ist dxf dem Be- 
trage nach größer als 1, so genügt es, den Integrationsweg bei 


= —ioo enden zu lassen, damit die gestellte Bedingung er- 


füllt wird. 


Die Integrale (99) mit dem so bestimmten Integrations- 
weg sind nach Heine!) Funktionen, die derselben Differential- 
gleichung genügen, wie die von uns durch Gleichung (21) de- 
finierten -Funktionen, ohne mit ihnen identisch zu sein; sie 
bilden vielmehr mit dem yw, ein Fundamentalsystem der 
Differentialgleichung 

d? w (w) 


dw* 


4 


1) E Heine, a.a.0.18.847 
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"nach den Kugelfunktionen und finden: 


F. Möglich 


Wir erhalten daher als Integrale der zweiten Art Reihen, die 
nach den Zylinderfunktionen der zweiten Art fortschreiten. 
Doch wollen wir als Entwicklungsfunktionen nicht die Funk- 
tionen (99) wählen, sondern eine ganz bestimmte Linearkombi- 
nation aus diesen und den w-Funktionen. Wir definieren: 


(100) W, (2) = — (2) + 

Wir erhalten dann, je nachdem der erste Index der Funktion 
hk, ,(9,p) gerade oder ungerade ist, eine Reihe, die nach 
geraden oder ungeraden W-Funktionen fortschreitet. Indem 
wir über eine willkürliche Konstante verfügen, schreiben wir; 


n 
Ry, = 7 
R® 4 > o 4 
Ron, ‚P) = RY (—1) A,, W, dn 8); 
a4 2n41v¥\ 9? 
Ro =3 —1)° (odn 
Ron iw AY RS) (K) ) 2041 +1 ( 6) 


Diese Reihen konvergieren sicher auBerhalb des Kreises 
|\dn8|=1, wie man aus der Betrachtung der Majorante 
erkennt, die durch Einsetzung der Maximalabschätzungen für 
die 4”, S. (643) gewonnen werden kann. 


Es ist übrigens nicht schwer, die analytische Fortsetzung 
für die Funktionen (101) zu finden, die über das Verhalten 
innerhalb des Einheitskreises Auskunft gibt. Zu diesem Zweck 


entwickeln wir die Funktion %, (w Yr? — 2rr, cosy + r?) 


= (—1)°y, „(or,) P,(cos 7); 


r<r,; cosy = cos cos 7, + sin sin cos(p — 
Diese Gleichung rechne man nun nach (23) auf kartesische 
Koordinaten um, indem man den drei Koordinaten r,, 9, 9, 
kartesische Koordinaten z,, y,, z, entsprechen läßt. Setzt 
man dann noch z, y,, z, gleich Null, so ergibt eine ähnliche 
Überlegung wie die im nee 5 durchgeführte, die 
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Y, (@ Vie —2,)?+ y’) 
= +11 10 VE 2) le: 
W (cw G2 (x,y). 


Setzen wir andererseits in die erste der Reihen (101) fiir die 
Größen 4) die Ausdrücke (32) ein, so ergibt sich: 


na 2n 


Veos? & + sin? $ cos® 


. „(cos +; ksin cosy) (wa Ro, p)sind dd dg. 


Ein Vergleich mit der eben abgeleiteten Reihe zeigt, daß 
beide, abgesehen vom Faktor R2,,,(+,q) und der Integration 
miteinander identisch wären, wenn auch noch die Glieder mit 
ungeraden Indizes vorhanden wären. Nun hatten wir aber 
gezeigt ($ 7), daß die Glieder mit ungeradem Index, wenn man 
sie mit einer Funktion %s,,(9,9) multipliziert. und dann 
über «+ von O bis a und über 9 von 0 bis 2” integriert 
verschwinden. Daher dürfen wir sie ohne weiteres hinzufügen 
und finden daher als on Fortsetzung der Reihen (101): 
(102) (8) = const J W, (@ - + sin? 
(A psinddddg. 


Setzen wir 8 gleich Null oder 2iK’, so sehen wir leicht 
daß die Funktion (102) an dieser Stelle eine logarithmische 
Singularität hat; denn dann verschwindet das Argument der 
Funktion für # = 0 oder a unabhängig von dem Werte von g. 


Da W(x) = 


-iz 


= ist, nimmt das er in der Nähe des 


Punktes oder = a die Form 
Behauptung sofort erfolgt. 


an, woraus die 


Doch auch stetig sind diese Funktionen nicht mehr. Dies 
liegt daran, daß die Wurzel Verzweigungspunkte hat. Ist 
nämlich dn reell und dem Betrage nach kleiner als 1, so 


sieht man leicht, daß die Punkte cos # = dnß und g = ~ 


solche sind. Doch geben nicht alle diese Anlaß zu Ver- 
zweigungspunkten der durch (102) dargestellten Funktionen, 
Annalen der Physik. IV. Folge. 83. 
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weil diese Pole im allgemeinen auf der Kugelfläche nur diskret 
liegen und von nicht höherer als der ersten Ordnung sind, 
In solchen Punkten konvergiert aber das Doppelintegral. Nur 
dort konvergiert es nicht, wo der Pol unabhängig vom Werte 
von » auftritt, also an den Stellen % = 0,”. Das sind aber 
gerade die Punkte, in denen die Funktionen der zweiten Art 
eine logarithmische Singularität haben. Prinzipiell könnte man 
von diesen beiden Punkten aus die Verzweigungsschnitte in 
beliebiger Weise verlaufen lassen, doch hatten wir bei der 
Definition durch Reihen dafür gesorgt, daß außerhalb des 
Einheitskreises Unstetigkeiten nicht auftreten können. Wir 
legen daher den Verzweigungsschnitt in die reelle Achse der 
dn ß-Ebene, zwischen die Punkte —1 und 1. Der Sprung der 
Funktion an diesem „Querschnitt“ ist gleich der Differenz der 
Funktionswerte für positive und negative Quadratwurzel, also 
unter Beachtung der Bedeutung der Funktion %#, (2): 


R® _ = {const f cos | (@ V(cos — dn)? +k? sin*® cos*g) 
| -V (cos? — dnp) + Bein’ cos'g 
V (eos# — + k’sin? cos*y 


Ru (4, gy) sin ddd dg 
Würde man übrigens in (102) die Funktion %, durch y, er- 
setzen, so würde man zu einer neuen Integraldarstellung für 
die Funktionen der ersten Art gelangen. Obgleich auch hier 
die Quadratwurzel verschiedene Vorzeichen haben kann, sind 
die Funktionen der ersten Art doch eindeutig, weil nämlich 
w, eine gerade Funktion ist. Die Differenz der Werte am 
Querschnitt ist also Null. 

Die eben gefundenen Funktionen der zweiten Art unter- 
scheiden wir von denen der ersten Art — die im Para- 
graphen 11 angegeben wurden — dadurch, daß wir ihnen als 
oberen Index eine 2 anfügen, während die Funktionen der 
ersten Art von jetzt ab als oberen Index eine 1 erhalten 
sollen. Die Flächenfunktionen aber, also die, die von zwei 
Argumenten abhängen, erhalten wie bisher keinen oberen 
Index, da hier eine Verwechslung ausgeschlossen ist: als 
Flächenfunktionen werden wir nämlich nur solche gebrauchen, 
die das Randwertproblem lösen. 

Wir haben gefunden, daß man aus den Funktionen der 
ersten Art die der zweiten Art herstellen kann, indem man in 
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den Reihenentwicklungen der ersten Art die Funktionen w, durch 
die Funktionen %, ersetzt. Doch gilt dies nicht allgemein, wenig- 
stens wurde es nicht allgemein bewiesen, sondern nur für die 
Reihen vom Typus (88), also für solche, bei denen das Argument _ 
der Zylinderfunktionen die Funktion w dn 8 war. Esistsehr wahr- _ 
scheinlich, daB fiir die Funktionen der zweiten Art auch Reihen- 
entwicklungen angegeben werden kénnen, die den Reihen (90) und 
(91) für die Funktionen der ersten Art entsprechen. Dazu müßte 
man in den Integralen auch der Integrationsvariablen einen 
komplexen Weg vorschreiben; dazu aber gehörte wiederum eine 
Kenntnis des Verhaltens der Funktionen 2,, ,(:+, p,) für sowohl 
große komplexe 9, als auch ebensolche g,. Da aber die 
wesentlichen Züge der Funktionen der zweiten Art sich bereits 
auf dem hier durchgeführten einfacheren Wege ergaben, glaubte 
der Verfasser, der eben ausgesprochenen Vermutung nicht 
mehr nachspüren zu brauchen. 

Doch bedürfen unsere Untersuchungen in anderer Hin- 
sicht der Vervollkommnung. Für die Funktionen der ersten 
Art genügte es, gewissermaßen als Beispiele den Fall der 
Funktionen R” (6) zu entwickeln; dann war es ohne weiteres 
klar, wie auch die anderen Funktionstypen zu behandeln seien. _ 
Hier aber ist insofern ein Unterschied, als die Methode ja 
wesentlich darauf beruht, eine Funktion von cos +, allein an- 
zugeben (in dem eben behandelten Falle die Funktion 
0,,,(cos 9,), mit der man eine gliedweise Integration aus- 
führen kann. Daß eine ähnliche Funktion auch für die anderen _ 
Typen angegeben werden kann, das bedarf noch eines Nach- _ 
weises, der aber leicht zu erbringen ist. 


Betrachten wir z. B. die Funktionen 7, ‚(9,%). Sie sind aes 


n,v 


Lösungen der Integralgleichung (53). Dort ersetzen wir wieder ’ I er 
und durch und 7; dann setze man y= A. So ergibt 


sich fiir die Funktionen der ersten Art: hea a 
az 
00 


Definieren wir nun in Analogie zu der Funktion 9, , (cos ,) 
die Funktion 
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(cos 9,) = sin 9, Po) Po Pos 


so erhalten wir unter Beachtung der Reihenentwicklungen (41) 
und (44): 


Tan, „(cos = sin ty > (cos 4 


= sin, Sot, +1 (608%) 
Man verifiziere ferner, daß die Gleichung 


P, (cos #,)sin? I, dit, = + 1) 
0 

gilt. So ergibt sich schließlich für die Funktionen der 
zweiten Art: 


(103) 7) = const - w 


Analoges für die Funktionen mit ungeradem ersten Index. 
uf Zwar könnte man leicht mit Hilfe einer der Rekursionsformeln 
der Zylinderfunktionen %# (z2)/z durch #,_,(x) und %#,,ı(e) 
ersetzen, doch ist es vorteilhaft gleich an die eben aufgestellte 
Reihe anzuknüpfen. Jetzt haben wir für die Größen B' 
gewisse Doppelintegrale einzusetzen; diese haben wir bisher 
explizit nicht angegeben, doch findet man sie leicht aus den 
Formeln (57) und (58), wenn man dort den geraden oberen 
Index 2m durch den ungeraden 2m + 1 und U, durch 7,, 
ersetzt. So erhalten wir: 


| 


= eonst- sn ß ff I enes+ 1)(40 + 1)te 
i 0 0 o=1 
Wy Veos* 9, + + k* sin? 9, ‘cos? ) (a@dnp) ib 
x Veos? 4, + k* sin* 9, cos? wo dnp 


(cos it, ; Asin, COS Tan,» (I, Fd Gy 


Andererseits ist: 
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> (20 + 3)(o + 1)(o + 2) Va 


v1 
9): 
Daraus schließen wir dann wie vorher auf die Richtigkeit 
folgender analytischen Fortsetzung: 


(2 


a 22 


9 sing YY, (oVieos dnp) + +h sin® + cos? 
3 V (cos + — dn + sin? cos? p 
Auch dies ist eine mehrwertige Funktion mit Polen in den 
Stellen daß =+ 1. Doch sind diese Pole nicht mehr loga- 
rithmisch, sondern von höherer Ordnung, 
Ganz analog kann man verfahren, um die Funktionen 
” (?) abzuleiten. Dazu beachte man nur, daß nach den auf 


(104) 


8.647 gemachten Bemerkungen die Funktionen 7, ,(:+,q) in 
die 8, ,(9,g) übergehen, wenn wir @ durch w/k’, k durch 


_ 7 und p durch oo gy ersetzen. Führen wir diese Sub- 
stitution an dem Kern der Integralgleichung (53) durch, so 
ergibt sich: 
ik 
Dieser Kern hat zu Kigenfunktionen sicher die 8, (9,9). 
Dann können wir genau so verfahren wie eben und erhalten 
schließlich: 
(o » Vıcos: + dn + k*sin* + cos? 9) 
sin 
-ffü V (cos $ — dn + k? sin? cos? 
-S§ sin dit do 


(105) 


00 


Mit Hilfe des Kernes pays Integralgleichung (53) und unter 
Beachtung der Gleichung “~~ = y,(z) verifiziert man mühe- 
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sin? + sin sin? iF, sin 2g, 08 cos 
Wo (@ ksin sin cos COs Y,) 
ein Kern ist, dessen Eigenfunktionen der Differentialgleich, (13) 
genügen. Dieser ist geeignet, um die Funktionen zweiter 
Art (0) herzustellen. Man findet: 
snß-en 
q 


(cos # — 6) + k* sin? cos? 


§ 13. Zwei Reihenentwicklungen 


Wie wir auf 8.630 gezeigt haben, hat die Integral- 
 gleichung mit dem Kern 


ei (dna cos 9008 d, + kcnasin sin 9, cos p cos py + iksna sin sin sin 


_ zu Eigenfunktionen die Integrale der Differentialgleichung (13), 

die dem gestellten Randwertproblem genügen. Denn die 
ar Größen dna, ken«, iksn« erfüllen die auf S. 628 ein- 
geführten Nebenbedingungn a?—c?=1 und 
Man kann daraus schließen, daß die eben genannte e-Funktion 
eine bilineare Entwicklung von folgender Gestalt zulassen wird: 


= = S Ey,» (@) Y) Po) 


Dabei bedeutet das Zeichen § eine Summation „über Buch- 
 staben“, d.h. man hat den Buchstaben R der Reihe nach 
durch 8, 7, U zu ersetzen, und auch an die Stelle der zu- 
nächst unbekannten Größen e, , irgendwelche anderen Ent- 
_ wicklungskoeffizienten gesetzt zu denken. Diese Koeffizienten 
sind natürlich Funktionen von «. Da aber die e-Funktion als 
Funktion von «, +, 9 der Schwingungsgleichung genügt, so 
müssen die Entwicklungskoeffizienten so bestimmt sein, daß 
jedes Glied der Reihenentwicklung der Schwingungsgleichung 
genügt. Die Koeffizienten von R,,(9,9) können daher, ab- 
gesehen von einem konstanten Faktor, nur die Funktionen 


RY (a) oder Re (2) sein. Da aber die e-Funktion eindeutig 
und überall endlich ist, so kommen nur die Funktionen der 
ersten Art in Frage, so daß sich als Entwicklung ergibt: 


Beug 


(107) 
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pio (dna cosd cosd, +k cn asin dein d, cos p Cos + ik sna sind sin sin y sin 
Hier bedeutet also S, daß man A, , durch 8, „7, 


n,v? 


und r,, 


U, und s,,, t,,. u,,, zu ersetzen habe. Die Reihenentwick- 
lung ist also bekannt, sobald wir die Koeffizienten r, ,, s, ,, 
t kennen. 


In»? Un,» 
”” Man beachte nun, daß die e-Funktion in den Kern der 
Integralgleichung (19) übergeht, wenn wir a = 0 setzen, Die 
Entwicklungskoeffizienten sind aber dann nichts anderes als 
die reziproken Eigenwerte dieser Integralgleichung. Mit anderen 
Worten: Für die Bestimmung z B. der Größen r,_, ergibt sich 


Da wir aber bereits auf S. 680, 682 und 684 die Berechnung 
dieser Eigenwerte auseinandergesetzt hatten, so dürfen wir r 
als bekannt ansehen. Analog verfahren wir zur Berechnung 
der s, ,, wobei wir auf die Durchführung der Berechnung der 
zu diesen Funktionen gehörenden Eigenwerte wohl verzichten 
dürfen, da eine Unklarheit hierüber nicht herrschen kann. 
Die Koeffizienten ¢,, sind auf diesem Wege nicht zu 


n,v 


erhalten, weil die Funktionen 7“) (a) für « = 0 verschwinden. 
Wir setzen deshalb jetzt «= K-+:iXK’ Dadurch geht die 
e-Funktion in den Kern e~ ©’ sin sin cos cos + sin sin yo) [Glei- 
chung (20a)] über; dann verfahren wir so wie früher. 

Ein wenig komplizierter ist die Berechnung der Koeffi- 
zienten uw, ,, Dazu betrachten wir nicht die eingangs an- 
irichens e-Funktion, sondern die Differenz aus dieser und 
der e-Funktion, in der g durch —p ersetzt wurde. Wir 
erhalten auf diese Weise die Funktion 


—2ioksnesin sin sin sin 


Nn,v 


(109) 


. + kcenasin sin cos po) 


Die dieser entsprechende Reihenentwicklung enthält natürlich 
nur noch solche Funktionen, die in bezug auf p ungerade 
sind; die Koeffizienten dieser Reihe sind daher entweder 24, , 
oder 2x, „ Die Reihe verschwindet für «= 0 und ist daher 


« 


| 
rar 
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0) 
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| 


F. Möglich 


b ae durch sn« dividierbar. Auch der Faktor 2 läßt sich auf 
beiden Seiten wegdividieren. Setzen wir dann «= 0, er- 
halten wir gerade aus der Funktion (113) den Kern der Inte. 
gralgleichung (53) (nur mit dem Faktor —wi%). Die Größe 
ak Mr a=0 
sind daher nichts anderes als die reziproken Eigenwerte der 
Integralgleichung (53), die sich in ähnlich einfacher Weise be. 
rechnen lassen, wie die Kerne (20) und (20a). Damit ist auch 
die Berechnung der u, | durchgeführt. 
Einer freundlichen persönlichen Mitteilung Hrn. Dr. 
a Hammersteins verdankt der Verfasser die Kenntnis der 
Tatsache, daß die bilineare Entwicklung eines Kernes K(z,7; §,1) 
wha im Fundamentalbereich der Variablenpaare x,y und &, 7 immer 
N ay dann konvergiert (mit Ausschluß des Randes sogar gleich- 
Br mäßig), wenn der Kern K stetige erste und zweite Ableitungen 
x Er, besitzt. Da diese Bedingung für unsere e-Funktion erfüllt ist, 
80 schließen wir auf die Konvergenz der Reihe (111). 
ER Die Gleichungen (85) lehren, daß, wenn man mit den 


2x 


) p)sind drags 


ay na 2n + 
f F(9,9) 8, g) sind dddg; 
4 


(A p)sinkdtdp 


bildet, diese im Falle der Konvergenz die Funktion /(?,¢) 
darstellt. 


Wir wollen uns die Aufgabe stellen, die Funktion 2, (7,7; ®) 
mit abgeändertem Parameterwert @ nach den Funktionen 


n,v 
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(9, 8, (9,93 0); (9, @); p; @) 


n, 
. vs. x 
zu entwickeln. Da w = q- — ist, so kann eine Änderung von o 


einmal dadurch entstehen, daß man g, d. h. die großen Achsen 
der konfokalen Ellipsoide ändert, oder dadurch, daß man die 
Schwingungszahl » ändert. Im ersten Falle würde sich nicht 
nur w, sondern auch k und #’, die elliptischen Moduln, ändern. __ 
Diesen Fall wollen wir ausschließen, da er kein physikalische _ 
Interesse bietet und nur den Fall betrachten, daß allein » 
und nicht & und A’ sich geändert haben. 
Setzen wir jetzt in (114) an die Stelle der Funktion #(3,4) 
die Funktion 2, „(9,9; ©), so sehen wir leicht, daß die Inte- 
grale (114 II); (114 III und IV) verschwinden. Die Integrale 
(114 I) sind nur von Null verschieden, wenn n gerade oder 
ungerade ist, je nachdem der erste Index der zu entwickelnden 
Funktion R, „(9,y; ©) gerade oder ungerade ist. 


Angenommen, der erste Index der zu entwickelnden " 


Funktion ist gerade, so finden wir unter Beachtung der Reihen- _ 
entwicklungen (31a): 


a 2x 


0 0 8 t 


m @ 

und unter Beachtung der Orthogonalität der Kugelfunktionen: ‘ BR 

20,* par NT 2 2s + 28)! 26,% 2t 2n,v 2 

8 t 


Damit ist die gesuchte Reihenentwicklung hergestellt. 
Die entsprechenden Koeffizienten fiir die Entwicklung der 


Funktionen.S, , p; @) usw. nach den 8, „(#, p; usw. nennen 


wir ); o); a; @). 
§ 14, Fälle der Ausartung 

In den bisherigen Entwicklungen hatten wir über den 
speziellen Wert des Moduls & nur die notwendige Voraus- 
setzung gemacht, daß er kleiner als 1 ist. Hier werden wir 
die spezielle Annahme machen, daß der Wert des Moduls & 
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P. Möglich. 


A 
: Wy entweder Null oder 1 ist. Dabei wird sich zeigen, daB sich 
er Pie in den Formeln große Vereinfachungen anbringen lassen, so 


können. Wir nehmen zunächst an, & sei gleich Null. 

j k=0 bzw. = 1 bedeutet nach (4), daß die Größen p 
und g einander gleich werden. Geometrisch heißt dies, daß 
die beiden kleineren Achsen der Ellipsoide übereinstimmen, 
80 daß diese in ein gestrecktes Rotationsellipsoid ausarten. 
Aus der Gleichung k?sn?u + dn?u = 1 entnehmen wir, daß 
für A=0 dn?u den Wert 1 annimmt, und dies in die Diffe- 
rentialgleichung (11) eingesetzt, ergibt: 


(112) +H+H-f=0. 


Dies aber ist nichts anderes als die Differentialgleichung fiir 
die trigonometischen Funktionen. Die Forderung der Ein- 
 deutigkeit im Raume verlangt, daß 
(113) ®+H+E=»% 
wird, wo » eine ganze Zahl bedeutet. 
= = Dies Ergebnis befriedigt einerseits, nimmt andererseits 
auch Wunder. Denn die Differentialgleichung (11) galt für 
A alle drei Koordinaten «, ß, y, und es ist ersichtlich, daß ein 
Produkt aus drei trigonometrischen Funktionen nicht ein Integral 
der Schwingungsgleichung in rotationselliptischen Koordinaten 
sein kann. Hier zeigt sich nun die Bedeutung der Tatsache, 
daß «, 8, y auf drei ganz verschiedene Intervalle beschränkt 
werden mußten. Auf S.619 hatten wir bereits darauf hin- 
gewiesen, daß wir auch reelle Koordinaten einführen könnten, 
allerdings unter Aufgabe der bemerkenswerten Symmetrie. 
Hier aber brauchen wir die Asymmetrie, und deshalb soll die 
damals angedeutete Transformation hier durchgeführt werden. 
Wir hatten auf 8. 619 für y das Intervall 
— K<y<—Kr2ik 
festgestellt. Zunächst bilden wir die Mitte dieses Intervalles 
auf den Nullpunkt ab, indem wir setzen: „=y+K-iK 
Für y’ ist das Intervall rein imaginär, so daß wir jetzt durch 
die imaginäre Transformation von Jakobi') auf ein rein reelles 
Intervall kommen. Es ist nämlich: 


€ 


si 
Pr 


af Whittaker, Watson, a.a.O. 8.505. 


Be 
Da 
dn 
gef 
jet 
= 
Inv 
» 
ha 
(11 
4 Di 
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ge 
en 
Gi 
AS fü 
ig da 
da 
va 
di 


a 


Daraus folgt dann für beliebige Werte von k der Übergang von 
dn(y,k) in — k'sn(n, k’), wo für 7 das Intervall — K’< 4 <+ K’ 
gefunden wird. Die Differentialgleichung (11) nimmt daher 


jetzt folgende Gestalt an: = 
+ (w? (yn; + HR? sn? (ns + Af =O. 


Hier kénnen wir nun den Rte vornehmen, indem 
wir A gleich 1 setzen. oad ist nun zu beachten, daß sn (n; k’) 
die zu 


inverse Funktion ist; diese aber geht für k’= 1 in . Area Tgy 
über, so daß wir jetzt die folgende Mimnellipeiikung er- 
(114) + (0? Tg!n + + H)-f=O0. Mere 


Die K’ läßt sich als Funktion des Moduls 4’ 
als hypergeometrische Reihe darstellen. Es ist?): 


(Übrigens gilt die gleiche Beziehung für X und k) Die hyper- 
geometrische Reihe ist innerhalb des Einheitskreises überall 
endlich, wächst aber mit Annäherung an diesen über alle 
Grenzen. Für k’=1 ist daher K’=oo und das Intervall 
für 7 ist daher - +m. 


Man sieht übrigens auch leicht auf dem gleichen Wege, 
daß X für k=O den Wert n/2 annimmt, und dies beweist, 
daß die Koordinate $ für k=0 in das Azimut 9 übergeht. 

Da der Wert von Tg innerhalb des betrachteten Inter- 
valles reell ist und zwischen —1 und 1 liegt, so kann man 


1) Whittaker, Watson, 2.2.0. S. 501. 
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durchführen. Bezeichnen wir cos? mit x, so können wir die 
Differentialgleichung (114) in folgende Form bringen: 

ge wir schließlich noch die Größe H durch F und » 


Lassen wir in Null gehen, so erhalten wir 


i aus (115) die Differentialgleichung der Zugeordneten der Kugel- 
funktionen. 

Um nun auch für die Koordinate & den entsprechenden 
Übergang zu finden, erinnern wir uns daran, daß der Mittel- 
punkt des Intervalles für vom Anfangspunkt des Intervalles 
für « um den Betrag X +:XK’ entfernt war, daß andererseits 


a k=0 K= und K’=oo unendlich wird. Daher wird 
dn(y +K+iK’) fir A=0 in den Ausdruck lim Cot(y + A) 


A>o 
_ übergehen; dieser aber ist gleich e-?”, Doch empfiehlt es 
sich, die Transformation 
e-?" = Sofa 
_ durchzuführen, weil wir dann durch die Gleichung (115) das 
gesamte Problem beherrschen. Für das gestreckte Rotations- 
ellipsoid gehen dann die Gleichungen (7) in folgende Form über: 


z = qGojacost; y= 


K, (0, 9; Po) = Kg (9, Fos = 

(l 17) | K, 07 p; Foy Po) = sin d sin % cos (p — Po), 

Der Kern (53) nimmt für diesen Grenzfall die Form 

(118) Ky Yo) = sin F sin p sin sin p, ei@ «os cos 
an. Schließlich ist auch leicht zu zeigen, daß der Kern 

f (119) Ky (9,9; Fos Po) = sin” sin n p sin” sin n gp, cos 8% 


_ Kigenfunktionen hat, die der Differentialgleichung (13) genügen, 
wenn man in ihr k = 0 setzt. 
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Die Integrale der Differentialgleichung (115) müssen sich, 
wie wir bemerkten, für w = 0 auf die Zugeordneten der Kugel- 
funktionen P”(z) reduzieren. Wir werden diese Integrale mit 


pe’ (z) bezeichnen und ihnen dieselben Indizes geben wie der 
Zugeordneten, auf die sie sich reduzieren, wenn ® Null wird. 
Nun ist andererseits sicher z. B. pe,” (cos #)c0os2v ein Inte- 
gral der Differentialgleichung (13), wenn wir in ihr k= 0 
setzen. Als solches aber bezeichnen wir es mit R,, | (i, p). 
Der erste Index der Funktion R(t, q) stimmt mit dem unteren 


Index von pe,” (2) überein; es fragt sich nur, in welcher Be- 


ziehung der obere Index 2» von p e,” (2) und der zweite Index o 
von R,„.(#,y) zueinander stehen. Dazu beachte man, daß 
der Index o von R,,, ,(,q) angab, zu welcher der nach der 
Größe geordneten Wurzeln der Säkulargleichung (79a) die 
Funktion R,,„(#,g) gehörte. Setzen wir nun in der Säkular- 


gleichung k = 0, so reduziert sie sich auf eine Form, in der 
nur die Diagonalglieder von Null verschieden sind. Die 


Säkulargleichung nimmt daher die Form | ar 
II 07 +3 P,) = 0 


e=0 
an, so daß die Wurzeln dieser Gleichungen gegeben 
‘ 
= 


. Po 

Andererseits sieht man sofort, daß zu der Funktion 
p e (cos F)cos2vp 

die Wurzel &,”’ gehört, so daß nur noch zu entscheiden ist, 

welche Reihenfolge die ¢,” annehmen, wenn man sie der Größe 


nach ordnet. Dies aber ist leicht zu entscheiden, wenn man 
beachtet, daß fir 4 = 0 die Größen it, den Wert 


\ „(an+2g! 2(8n* — 8" + 4n—1) 
(2% — 20)! (16n* — 1)(4n +3) 


die Größen p, den Wert =" : ae annehmen [nach 


8S. 664, 669 und den Gleichungen (60) und (72)]. So erhalten 
wir pämlich: 

=— (Su? — 8v? + 4n — 1) 


as 

| re 

| 
| 

‘Mele 

= 
| 
sinddurch 
A, 
; % 


pe,°(z) finden wir also die Integralgleichung: 
pe,*(z) =i, | eio2ype(y)dy. 


= (= (4 +1) (wy)p 


| 
pe (0) 1+ 4.5 
fied 
der Eigenwert %,: 
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und der Wert dieser Größe wird mit wachsendem » immer 
größer, da ja immer n=» ist. Dadurch ist gezeigt, daß die 
natürliche Ordnung der ¢,, d.h. die nach den oberen Indizes 
der Funktionen pe,” (2) zugleich auch die Ordnung nach der 
Größe ist, und daraus folgt dann sofort die gesuchte Beziehung 
zwischen dem zweiten Index der Funktion R,, (9,9) und 
dem oberen Index der Funktion pe,” (a): 


(120) pe,” (cos = 
Analoges er auch für die anderen Funktionstypen: 


(12 ) ) | | ) + m y ); 
( a 2n+1 


2n,¥ 


(7, p) usw. 


Wir wenden uns nunmehr den Integralgleichungen zu. 
Man sieht sofort, daß der Kern ei@«s#c0s% nur die Funk- 
tionen R, = pe,’ (cos #) zu Eigenfunktionen haben kann, 
weil dies die einzigen sind, die von g unabhängig sind. Für 


—1 


Schreiten wir nun zur Berechnung z. B. der Funktion pe,’ (z), 
d.h. derjenigen, die sich für m=O auf die Kugelfunktion 
(2) = 1 reduziert. Wir setzen also an 


Mit aiden. Ansatz finden wir sehr leicht: 


+1 


_ Ferner ergibt sich aus der Gleichung 
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(123) dy = + pe°(0) = 


oder da 

ist 


Da nun die Funktion y,, als kleinste Potenz von m die 
Potenz w?* hat, so werden wir schließen, daß auch die ge- 
suchten Koeffizienten A,, als kleinste Potenz von » die 
Größe w*° haben, und wir machen den Ansatz: 


indem wir Potenzen von höherer als der (2¢ + 4). Ordnung 
unterdrücken. Weiter werden wir in (121) A,; w,, (wy) und Br 
pe,°(y) durch ihre ersten drei Glieder annähern und erhalten so: fe as Ba 


(= 1)" + boo + C40") 77 
Vr(4s + 3) 
+2,0*+ c,0°) + is (a,o* +2,0°.. ) a 

+1 99 5) 
1,0 ,w' + ¢,@ 
+ + b, w* + e, 2 
+ [a, +b,0 +...) P,(y)}dy. 
Aus dieser Gleichung lassen sich dann in verhältnismäßig PR, ia 
einfacher Weise die Koeffizienten «a, ; b,,; «,, ausrechnen; 
denn die Gleichung (124) besagt natürlich, daß die Koeffizienten i. . 4 
gleicher Potenzen von w auf beiden Seiten gleich sein miissen. oa : Br 
Um diese Vergleichung durchzuführen, brauchen wir nur noch a 
die folgenden WORAN unter denen auch die ist, die man zur rs 
Berechnung der pe, „ı(2) mit ungeradem Index gebraucht: a 
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Jv? Pulydy = Ve ap)! 
1 k)! I(p+k+4) 
+1 0 wenn p<k 
dy = (2p + 1)! 
“1 


der Gleichung 
+ 1 
22” w,(o) =| eie:P (z)dz 


—1 


Die Richtigkeit dieser Formeln beweist man, indem — 


beide Seiten durch ihre nach Potenzen von g fortschreitenden 
Reihen ersetzt und dann die Koeffizienten gleicher Potenzen 
von o gleichsetzt. 


Ganz ähnlich verführt man zur Berechnung der Funktion 
pe(@). Diese reduziert sich definitionsgemäß auf die Kugel- 
funktion P,(z), wenn man w = 0 setzt. Daher ist also in der 
nach Kugelfunktionen fortschreitenden Reihe nicht mehr 4, 
sondern 4, gleich 1 zu setzen. Dann gestaltet sich die Be- 
rechnung des Eigenwertes A, in folgender Weise: Man bilde 
die Funktion 


-a!(n — a)! (n — 20)! \iw 24 
m= hn — 1), je nachdem n gerade oder ungerade. 
Da 


(2 - 3 )! 
(2) > (— 1) a! (n = = 20)! 


iwzy P 
22772777 

Setzen wir jetzt x = 0, so erhalten wir wegen der FERN 
der Kugelfunktionen: 


4 Im übrigen vollzieht sich die Rechnung in der gleichen Weise. 
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Als Beispiel einer auf diesem Wege gewonnenen Eigen- 
funktion führen wir an 


x 2 9 9 
= | + + 1) 3) 82. (20 + 3) 


(1296 -8+0° + 6778+8-04 — 9117-4-0°— 24663 0?— 16065 -0) 


22.39.52. 7-.(20 + 3)(20 + 5) + $) 


(2) 
Die Berechnung der Funktionen pe’ (x) gelingt mit Hilfe 
des Kernes: 
— »sindsin d,cos(p — qo) 


Aus der 


pe (x) cosyg = [fe — @ Vi- at Vi con(y - 

-1 0 “pe” cosy p, dy dg, 

folgt, wenn man die Integration über g, ausführt und dann 
noch durch cos» p dividiert: 


(127) pe” (x) [7,60 Vl — 2? Vl —y*)pe’ (yay. 


Um diese Integralgleichung in der gleichen Weise zu be- 
handeln wie die Gleichung (121), brauchen wir vor allem die 
Entwicklung des Kernes J, (iw — 2x? ¥1—y?%) nach den Zu- 
geordneten P” (x), wo oberer Index der Kugelfunktion und 


Index der Zylinderfunktion iibereinstimmen. Es handelt sich 
also um eine Entwicklung bei unveränderlichem oberen Index 
der Zugeordneten. Die Koeffizienten einer solchen Entwicklung 
finden wir durch das Integral 


2n+1 (n— fier a (4) dz, 


wo die zu entwickelnde Funktion Auf unseren 


Annalen der a IV. Folge. 83, 


f er 
> 
» 
vl 
7 
2 
72 V1 _ „a pr = 


dies aber haben wir — sogar in einer Verallgemeinerung — 
gu berechnen gelernt [Gleichung (30), in der s = 0 zu setzen ist] 
Wir können die eben angewandte Methode auch hier benutzen, 
so daß sich eine nochmalige explizite Darstellung erübrigt. Wir 
brauchen zur Durchführung der Rechnung die Formeln: 


Jom + 1 (dw yi — 2 yl —y’) 


= (20 — 2m)! (20 + 3) 2n+1 


 o=m 


#= (20 +3) yl — y?) (2) 


(a — m)!(o + m)! 


+1 


128 28 ..9 
—1 


(k — m)!(k + m)!(2n + 2m)! 
=o (n + m)!(n — m)!(k — n)! Wk +n + 4) wenn in 


— m)!(k + m + 1)!(2n + 2m + 2) 
+ m+ 1)!(n — m)! (k — n)! +) wenn k=n 


wenn mSk<n ist, sind die beiden letzten Integrale Null. 
Wir erhalten auf diesem Wege z.B.: 


| 80"? +?. a(4o+3) 


(1024944 -0° +6208 016-0% +179 925-0 +397 560) 
2°°+7.54 (20 + Hor + i) 


@)- 
Durch die Integralgleichung (127) erhalten wir nur die- 
jenigen der pe*(z), bei denen die Differenz zwischen dem 


unteren und dem oberen Index eine gerade Zahl ist. Zur 
Berechnung der noch fehlenden benutzt man z.B. den Kern 
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(130) z.y-J, (io Vi —ı? yi-y?) 
oder auch einen der Kerne (119). 

Zur Auswertung des eben genannten Kernes brauchen wir 
Formeln, die den Gleichungen (128) analog sind. Diese lassen 
sich leicht aus den Gleichungen (128) herstellen, wenn man 
die Gleichung (60) auf S. 657 beachtet. Wir führen nur eine 
der so zu gewinnenden Formeln hier an: 


yl—2? yl - 
IE >> (Qo — 2m + + 
wo Vi Am — m)! (0+ m + 1)! 


Bei Betrachtung der Formel (129) könnte man über die in 
dem dritten Summanden auftretenden sehr großen Zahlen 
erstaunt sein und Zweifel an der Konvergenz hegen. Es sei . 
daher hier eingeschaltet, daß c, von der Größenordnung 7}; ist, i 
während bereits c, kleiner als 4-10 ist. 

Die Funktionen pe’ (z) gestatten nicht nur eine Entwicklung | 
nach Kugelfunktionen bzw. nach den Zugeordneten, sondern 
auch nach den Zylinderfunktionen w,(az). Diese Reihen- 
entwicklung ist der vollkommen analog, die wir für die 
Funktionen 2°) («) usf. im Paragraphen 11 gefunden haben, 
Wir wollen die Entwicklung nach Zylinderfunktionen für die 
pe‘ («) kurz ableiten. Dazu setzen wir in die Integralgleichung (121) 
unter dem Integralzeichen für pe. (y) die bereits gefundene 
Entwicklung nach Kugelfunktionen an. Auch entwickeln wir 
den Kern nach Kugelfunktionen der Variablen y und integrieren 
dann gliedweise. So ergibt sich, wie man leicht nachrechnet: 


| pes, (2) 2A, , 1)’4,, W, , (@ 2); 


wo die Koeffizienten A, die gleichen sind, die auch bei der 


(132) 


w 


Entwicklung nach den Kugelfunktionen auftreten. Ferner gibt % = 
es Entwicklungen nach den -Funktionen mit dem Argument _ lal 
47* 

ae 


= 
ir 
= 
r) 
ink 
N 
. 
1 
1 


F. Möglich 
iw Yl—2*. Um diese zu erhalten, würde man die Kerne 
(iw bzw. z-y-J, (io V1 — x? V1 — y?) benutzen, 
und um die pe” (x) nach den w-Funktionen mit dem Argument wz 
re oa. zu entwickeln, verwendet man die aus (119) folgenden Kerne 
Diese ergeben die gesuchten Reihen- 
entwicklungen, in denen jede Zylinderfunktion noch mit dem 


In analoger Weise wie im Paragraphen 13 zeigt man dann, 
daß man zu Funktionen der zweiten Art gelangt, wenn man 
die Zylinderfunktionen y, durch die Funktionen ersetzt. 
Wir erhalten so Reihen, die genau wie früher überall außerhalb 
des Einheitskreises konvergieren und durch Summation in den 


bezeichnen sie mit ge’(z). Den im Paragraphen 13 durch 
 Summation gewonnenen Integraldarstellungen für die all- 
gemeineren Funktionen können wir hier die folgenden Aus- 
drücke an die Seite stellen: 


(188) ge” (z) = const. Yi — 2” m, dy. 
[o (a — ur? 

Die nach Zylinderfunktionen fortschreitenden Reihen 
=. wir im folgenden immer dann benutzen, wenn das 
Argument der Funktionen pe’ (x) größer als 1 ist. Denn die 
Zylinderfunktionen sind auch für unendlich große reelle Werte 
ihres Argumentes noch endlich, so daß die mit ihnen ge- 
bildeten Reihen auch noch für unendlich großes Argument die 
Diskussion über das Verhalten der Funktionen pe’ (x) gestatten. 

Es sei noch kurz darauf hingewiesen, daß ein Teil der 
in diesem Paragraphen abgeleiteten Ergebnisse bereits von 
Abraham und in größerer Allgemeinheit von Niven), aller- 
dings auf einem ganz anderen Wege gefunden worden sind, 
wie wir auch kurz in der Einleitung andeuteten. 

Ein Integral der Schwingungsgleichung ist nach den jetzt 
angestellten Überlegungen das Produkt = 


(134) pe™ (Cojn)pe” (cos 9)cosmy. 


Vgl. die Literaturangaben auf S. 609 und 610. 


Er Dal 
nat 
pe. 
fra, 
elli 
Gr 
Ku 
für 
da 
r ga 
fu: 
os Einheitskreis hinein analytisch fortgesetzt werden können. Diese ell 
oo Funktionen haben in den Punkten 1 und — 1 einen Pol. Wir sit 
di 
A 
wi 
d: 
sc 
W 
re 
is 
q 
I 
e 
( 
] 


Beugungserscheinungen an Körpern v. ellipsoidischer Gestalt 709 


Dabei denken wir uns die Funktion pe” (Cof7) durch ihre 
nach Zylinderfunktionen fortschreitende Reihe, die Funktion 
pe"(cos%) durch die Kugelfunktionsreihe dargestellt. Wir 
fragen, was aus dem Produkt wird, wenn wir vom Rotations- __ 
ellipsoid zur Kugel übergehen, d.h. wenn wir in o= i. die 
Größe g gegen Null gehen lassen. Offenbar wird dann auch 
w gleich Null, und damit verschwinden in der Reihe nach 
Kugelfunktionen alle Koeffizienten A” bis auf einen 4”, der 
für o = 0 gleich 1 wird. Das Produkt pe” (cos :#) cos m ıy geht 
daher in die Kugelflächenfunktion P” (cos +) cosmy über. Nicht 


ganz so einfach ist zu entscheiden, was aus der Zylinder- _ 
funktionsreihe wird. Zwar reduziert auch sie sich auf ein 
einziges Glied, da ja die Entwicklungskoeffizienten die gleichen 
sind wie in den Kugelfunktionsreihen. Es bleibt somit von a 


> 


z 


dieser Funktion nur w,(®@oj»), doch steht jetzt auch im 
Argument die gegen Null gehende Größe w, so daß untersucht 
werden muß, was lim (*# Cof n) bedeutet. Man beachte dazu, 

7 
daß g gewissermaßen das natürliche Längenmaß ist und 7 das 
Ergebnis einer Messung mit diesem. Wird das Maß kleiner 
so wird natürlich das Messungsergebnis größer, d.h. n wächst, 
wenn g gegen Null geht, über alle Grenzen. Andererseits 
rechnet man aus (116) leicht nach, daß 


r= yr? + y?+ 2? =g: — — cos? it 
ist. Für sehr große Werte von n unterscheidet sich daher 
q&ojn beliebig wenig von der Größe r, der Entfernung eines : 
Punktes vom Nullpunkt des Koordinatensystems, so daß sich _ 


ergibt: lim (74 Gof n) = =~, Damit finden wir also durch ~ 

q->0 
Grenziibergang aus dem gestreckten Rotationsellipsoid ein ead 
Integral der Schwingungsgleichung fiir Polarkoordinaten: 
17 r) P™ (cos #)cosm . 


rn 


Dies aber sind in der Tat die für diesen Spezialfall bekannten 
Integrale. 


ig 


Der andere Grenzfall, A = 1, ist ohne weitere Überlegung Es u 
auf den eben behandelten zurückführbar. Dazu betrachten Ry ne 


e 
1, 
r 
e 
- 
| 
| 
| 
2 7 
| 


von = bis = +ioo und © von -7 bis = wahrend die 


(116) 


Im Falle des gestreckten Rotationsellipsoides läuft dann a 


Größen g durch ig zu ersetzen sind. 


 Schwingungsgleichung anknüpft, vollständig durchgeführt und 


Angriff zu nehmen. 


F. Möglich 


wir eine feste Ebene im Raume, die wir zur zy-Ebene machen. 
In ihr zeichnen wir eine Ellipse, deren große Achse mit der 


z = y= gEojncos#+cosp; 
z = gCojncos 


z-Achse und deren kleine Achse mit der y-Achse zusammen- ell 
fällt. Ist q der Abstand eines Brennpunktes vom Mittelpunkt ms 
der Ellipse, so kann man, wie Sieger’) zeigte, z und y aus- ph 
drücken durch die folgenden elliptischen Koordinaten: h de 
z=glojncos y= gSinnsin #. q 

Dreht man jetzt die Ellipse um die z-Achse, so kann man in an 

der Ebene senkrecht zu dieser ein Azimut » einführen, mit au 
dessen Hilfe sich die rechtwinkligen Raumkoordinaten z, y, z gl 
gerade gemäß den Gleichungen (116) berechnen lassen. In je 
der Tat ist bei der Drehung der Ellipse ein gestrecktes ß 
Rotationsellipsoid beschrieben worden. Anders wird aber die g 
Sache, wenn wir die Drehung um die y-Achse ausführen. d 
Hier wird ein abgeplattetes Rotationsellipsoid beschrieben, g 

und die diesem entsprechenden Koordinaten sind die folgenden: d 

d 

y= qSinysind; x = q cosg; z = Cojy cos? sing 
oder, wenn man die z-Achse als Drehachse benutzt: 

| 


Gleichung (116) und (116a) können zu einer zusammengefaßt 
_ werden, durch folgende Gleichung: may 


z =— igsinasinOsinp. 


von 0 bis 100, © aber von 0 bis z, im anderen Falle läuft a 


Wir haben hiermit das Eigenwertproblem, das an die 


sind jetzt in der Lage, das physikalische Te | in 


1) B. Sieger, Ann. d. Phys. 27. S. 567. 1908. ic 
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-gemeinsten Integral der Schwingungsgleichung [Gleichung (15)] % 
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Zusammenfassung des mathematischen Teiles 

Im § 1 werden zunächst in der bekannten Weise die 
elliptischen Koordinaten A, u, » eingeführt. Von diesen geht 
man über zu einer anderen Art elliptischer Koordinaten, die 
physikalisch als Potentiale der Flächen zweiten Grades ge- 
deutet werden können, die also der Laplaceschen Differential- 
gleichung genügen. Die Umrechnungsformeln der kartesischen 
Koordinaten auf diese wird gegeben. Dann wird (§ 2) die Form 
angegeben, die die Schwingungsgleichung in diesen Koordinaten 
annimmt. Wir suchen spezielle Integrale der Schwingungs- _ 
gleichung, die ein Produkt dreier einfacher Funktionen sind; _ 
jede dieser Funktionen soll nur von einer der drei Koordinaten a, _ 
8, y abhängen. Durch diesen Ansatz wird man zu drei linearen 43 
gewöhnlichen Differentialgleichungen zweiter Ordnung fir die = 
drei Funktionen geführt, in denen zunächst je zwei, also ins- 
gesamt sechs, willkürliche Konstanten auftreten. Die Form © 
der drei Differentialgleichungen ist die gleiche. Die Bedingung, 
daß diese drei Differentialgleichungen wirklich der Schwingungs- _ 
gleichung äquivalent sein sollen, reduziert diese sechs Kon- Er | 
stanten auf zwei: je drei sind immer gleich. Daraus erkennt 
man: die drei Funktionen, deren Produkt ein Integral  . ae 
Schwingungsgleichung ist, genügen einer und derselben Diffe- fer Be. 
rentialgleichung [Gleichung (11)]. Das verschiedenartige Ver- 
halten dieser Funktionen wird dadurch gewährleistet, daß die 
Argumente «, f, y auf verschiedene Intervalle beschränkt sind. 

Zuerst betrachten wir das Produkt /(f)- f(y), wo die 
Argumente ß, y auf die im § 1 näher bestimmten Intervalle _ Fr 
beschränkt sind. Dieses Produkt muß, wenn anders es über- > ae 
haupt zur ane des Problems in F rage kommen soll, ge- Eich 


von und y, von den „Winkelkoordinaten“ be- 
trachtet. Als Funktion dieser Koordinaten ist das Produkt 
ein Integral der Differentialgleichung (13). Die hier noch auf- 
tretende Konstante Z’ soll so bestimmt werden, daß das In- _— 
tegral auf einer Kugelfläche endlich eindeutig und stetig ist. _ 
Die Differentialgleichung (13) in Verbindung mit dem all- _ 


führt zu mehreren Integralgleichungen, deren Eigenfunktionen RS A i 


n, 
er 
kt 
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Möglich 


die gesuchten Integrale der Differentialgleichung (13) sind, 
Aus diesen Integralgleichungen folgt dann die Existenz von 
Reihen, die gleichmäßig konvergent auf einer Kugelfläche, 
nach Kugelflächenfunktionen fortschreiten und die gesuchten 
Funktionen darstellen. Bei der Darstellung dieser Reihen 
ergibt sich eine Klassifikation der gesuchten Funktionen: Acht 
verschiedene Reihentypen werden für sie gefunden. 


Im § 10 wird das Ritzsche Verfahren auf unsere Diffe- 
rentialgleichung (13) angewandt. Durch geschickte Ordnung 
der dabei auftretenden Determinanten fällt der gewöhnlich 
diesem Verfahren anhaftende Übelstand, nur Näherungsresultate 
zu liefern, fort, und wir können den Satz exakt gewinnen, daß 
alle Eigenwerte E’ unserer Differentialgleichung einfach sind. 
Mit diesem Satze gelingt weiter leicht die Angabe der Integrale 

der Differentialgleichung (11). Für diese findet man Reihen, 
die nach Zylinderfunktionen fortschreiten. Diese konvergieren 
auch in dem Intervall noch gleichmäßig, das wir für das 
Argument @ angegeben haben. Daher liefern diese Reihen 
präzise Aussagen über das Verhalten unserer Funktionen im 
_ _Unendlichen. Der § 12 bringt eine Darstellung der Integrale 
a zweiter Art der Differentialgleichung (11); diese werden bei 
i om der Darstellung der Beugungswelle Verwendung finden. Der 
814 behandelt den Spezialfall, der durch das Vorhandensein 


einer Symmetrieachse im Koordinatensystem gekennzeichnet ist 


§ 15. Präzisierung der Aufgabe; Koordinatensysteme; 
Hertzscher Vektor 

Die Aufgabe, die in diesem zweiten Teil behandelt werden 
soll, können wir so aussprechen: In ein beliebig kompliziertes 
_ elektromagnetisches Feld, das im Endlichen keine Singularitäten 
enthält, werde ein materieller Körper von ellipsoidischer Ober- 
 fliche hineingebracht. Die Materialkonstanten, Dielektrizität 
und Leitfähigkeit seien bekannt. Die Größe des Körpers sei 
a mit der Wellenlänge des einfallenden Lichtes vergleichbar. Zu 
Yin untersuchen ist der Zustand, der von dem Körper hervor- 
Rn gerufen wird. 
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Ein Koordinatensystem, das sich als nächstes zur Fixierung 
der Raumpunkte darbietet, ist durch die drei großen Achsen 
des Ellipsoides bestimmt. Wir können z. B. festsetzen, daß 
das Koordinatensystem rechtshändig sein, daß die x-Achse mit 
der größten, die z-Achse mit der kleinsten Achse des Ellip- 
soides übereinstimmen soll. Zwar ist auch dadurch die Lage 
eines Koordinatensystems noch nicht eindeutig bestimmt, weil 
sich im allgemeinen mehrere solche finden lassen, die diesen 
Bedingungen genügen; doch spielt dies keine Rolle, da wir 
unter den möglichen uns eines willkürlich aussuchen dürfen, 
und wir wollen annehmen, daß dies geschehen sei. Die Punkte 
des Raumes haben in bezug auf dieses Koordinatensystem die 
Koordinaten z, y, z. Die Oberfläche des beugenden Ellipsoides 
ist gegeben durch die Gleichung: ee me 

wo a, b, c drei uns bekannte Größen sind. Suchen wir diese 


Gleichung mit der Definitionsgleichung der elliptischen Ko- 
ordinaten (§ 1) zu vergleichen, so finden wir: 


Mit p und g sind uns dann auch die Moduln unserer ellipti- 
schen Funktionen gegeben [vgl. Gleichung (4). In den Ko- 
ordinaten «, ß, y, wie sie auf S. 614 eingeführt wurden, schreibt 
sich dann die Gleichung für die Oberfläche des beugenden 
Ellipsoides in der einfachen Form: 


bestimmt [Gleichung (5)]: 


(135) a=g-dna,; - 


Der Wert des Quotienten eu ist immer größer als 1 


/ 


und auch reell, da voraussetzungsgemäß a größer als c ist. 
Daher wird durch diese Gleichung auf der imaginären «-Achse 
zwischen 0 und der Wert von a, eindeutig definiert. 

Die Maxwellschen Gleichungen 


&+4nx&=crotd; nH =—crov€; divE=0; divG=0 


73 
h 
2 
E 
4 
ax Er 2,2 
a A 
Ar: 


= 


Lee 
lo 


714 F. Möglich 


1 


gehen für zeitlich periodische Zustände (Zeitfaktor e’”‘) in 
folgende Form über: 


dsv 


(136) 


divE=0; divH=0. 


Ein Integral dieser Gleichungen erhält man aus einem Vektor ®, 
dem Hertzschen Vektor, dessen rechtwinklige Komponenten 


2 
der Schwingungsgleichung AB + Fe PB = 0 genügen, dadurch, 
daß man auf ihn die Operation rot bzw. rot rot anwendet. 
Es ist: 


c 
Ein anderes Integral erhält man allein durch Vertauschung 
von & und §. 
In den älteren Arbeiten zur mathematischen Theorie der 
Beugung wird der Übergang zum Hertzschen Vektor ver- 
mieden. Dort geht man so vor, daß man für die zur Ober- 


fläche des beugenden Körpers normalen Komponenten der 


elektrischen und magnetischen Feldstärke aus den Maxwell- 
schen Gleichungen durch Elimination der anderen Komponenten 
je eine partielle Differentialgleichung herleitet. Prinzipiell 
würde wohl auch in diesem Falle dieser Weg gangbar sein. 


Hertzfeld z. B., von dem der erste und bisher einzige Ver- 


such zur Theorie der Beugung an Ellipsoiden herrührt, ist ihn 
gegangen, wenigstens für den speziellen Fall des Rotations- 
ellipsoides, — der allgemeinste Fall ist nie in Angriff ge- 
nommen worden. In der Hertzfeldschen Arbeit macht sich 
aber außerordentlich störend bemerkbar, daß die auf diesem 
Wege erhaltene Differentialgleichung nicht mehr separierbar 
ist, d. h. ihre Integrale sind nicht mehr durch Produkte aus 
drei Funktionen je einer Raumkoordinate darstellbar. Bereits 
für den Fall unendlich großer Leitfähigkeit wird die Durch- 
führung der Rechnung außerordentlich kompliziert. Ein wirk- 
licher Fortschritt konnte nur dann erzielt werden, wenn es 
gelang, eine Differentialgleichung der Betrachtung zugrunde 
zu legen, die in dem eben erklärten Sinne separierbar ist. 
Eine solche ist aber die Schwingungsgleichung und Integrale 


dieser sind die rechtwinkligen Komponenten des Hertzschen 
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Vektors. Andererseits sind die elliptischen Koordinaten wichtig 
für die Grenzbedingungen. Diese sind ja die Forderung nach 
stetigen Tangentialkomponenten von © und § auf der Ober- 
fläche des beugenden Ellipsoides. Welch einfache Form die 
Gleichung der Oberfläche in elliptischen Koordinaten annimmt, 
hatten wir gesehen. ‘Dies ist der Grund für die Einführung 
der elliptischen Koordinaten. Somit handelt es sich darum, 
Beziehungen herzustellen zwischen den rechtwinkligen Kom- 
ponenten des Hertzschen Vektors und den elliptischen Kom- 
ponenten von © und 9. Diese Beziehungen müssen hin- 
reichen, um aus den elektromagnetischen Grenzbedingungen 
die rechtwinkligen Komponenten des Hertzschen Vektors ein- 
deutig zu bestimmen, wenn dasjenige Feld bekannt ist, das 
geherrscht hat, bevor wir den beugenden Körper hinein- 
brachten. Doch soll gleich hier bemerkt werden, daß sich 
gerade diese Forderung im allgemeinen nicht hat erfüllen 
lassen. ‚Die Grenzbedingungen wurden so kompliziert, daß 
kaum eine Aussicht auf Durchführung bestand. Nur für den 
Fall absoluter Leitfähigkeit erwies sich die Benutzung des 
Hertzschen Vektors als zweckmäßig. Der allgemeine Fall 
läßt sich dadurch behandeln, daß man direkt die rechtwinkligen 
Komponenten von € und § benutzt; dann hat man allerdings 
noch den Divergenzbedingungen zu genügen. Die Einzelheiten 
dieses Weges sind später ausgeführt. 

Der mathematische Ausdruck des ursprünglichen Feldes, 
das geherrscht hat, bevor wir den beugenden Körper 
hineinbrachten, soll singularitätenfrei sein. Die Abhängigkeit 
von der Zeit ist periodisch mit der Frequenz v. Den all- 
gemeinsten Ausdruck für ein solches Feld hatten wir bereits 
in der Gleichung (15) aufgeschrieben. Diesem wollen wir eine 
etwas andere Form geben: 


zx 2x 


/ i (z cos d, + y Sin Dp COs Po + z Bin Dy sin Po) 
F (x, y; z) = f f ec 0 0 


| 
Über die Winkelfunktion machen wir die Voraussetzung, daß rt oe 
sie sich in eine auf einer Kugelfläche gleichmäßig konvergente, 
nach den Funktionen (Hy); 


n,m 
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F. Möglich 


annehmen es sei: 


un m 


Dann ergibt sich durch gliedweise Integration eine Reihe für 
eine Komponente des Hertzschen Vektors; wenn man noch 
die Gleichung en beachtet: 


= 


n m 


U, ,m(#, Po) fortschreitende Reihe entwickeln läßt. Wir 


Po) 


F(x, y, z (4, @); 


~ 


187) (a) R,,m 


= 7, 


Durch eine einzige Komponente des Hertzschen Vektors ist 
bereits ein möglicher elektromagnetischer Zustand bestimmt. 
Die Hinzunahme einer zweiten und einer dritten Komponente 
würde je ein anderes mögliches Feld geben, die sich dem 
ersten überlagern, ohne es zu stören. Es ist daher keine Be- 
schränkung der Allgemeinheit, wenn wir annehmen, der Hertz- 
sche Vektor des ursprünglichen Feldes habe nur eine einzige 
Komponente; der allgemeine Fall geht dann durch Super- 
position dreier spezieller Felder hervor. Im folgenden soll 


n,m n,m* 


daher angenommen werden, daß das ursprüngliche Feld sich 
aus einer einzigen Komponente des Hertzschen Vektors an- 
geben läßt. Weiter nehmen wir an, daß das Feld durch (137) 
gegeben sei. 
r _ § 16. Die Erfüllung der Grenzbedingungen 
- Das Beugungsfeld muß im Unendlichen verschwinden und 
darf dort nur divergierende Kugelwellen darstellen. Der 
analytische Ausdruck für eine Kugelwelle ist, wie wir hier be- 
. See weislos bemerken wollen, der folgende: 
Daher dürfen für die NDR nur solche Ansätze ge- 
macht werden, die im Unendlichen in diesen Ausdruck über- 


gehen. Es ist nun unschwer einzusehen, daß die im § 12 
eingeführten Integrale der zweiten Art diese Eigenschaft 
haben. Wir wollen den Beweis für diese Behauptung für die in 
der Gleichung un) dar, Funktion RY (a) kurz durch- 
führen. 
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In der Theorie der Besselschen Funktionen leitet man 
für die Zylinderfunktionen eine asymptotische Darstellung her, 
die Auskunft über das Verhalten der Funktionen bei sehr 
großem Argument gibt. Diese Darstellungen sind im all- 
gemeinen semikonvergente Reihen, die an geeigneter Stelle 
abgebrochen werden müssen, dort, wo der Fehler klein ist. 
Diese Reihen brechen nun dann von selbst an endlicher 
Stelle ab und machen daher Fehlerabschätzungen über- 
flüssig, wenn der Index der Zylinderfunktionen ein ungerades 
Vielfaches von !/, ist. Dieser Fall aber liegt gerade hier vor; 
die Funktionen % (2) sind ja, wie wir gezeigt haben, ab- 


gesehen von dem Faktor Te ‚die Hankelschen Zylinder-* — 
x 
funktionen Ben (x), für die sich na Darstellung ergibt'): 


(2) = Bi + R, (2), 


bedeutet. Diese Darstellung der Funktionen ¥ (x) setzen wir 
dann in die Gleichung (101) ein. Es läßt sich nun zeigen, 
daß die Reihe der „Beste“ 


außerhalb des Einheitskreises gleichmäßig konvergiert und 
daher mit wachsendem Argument gegen Null geht, weil jeder 


der „Reste“ dies tut. Daraus folgt dann für sehr großes 


wo 


Argument für RS) „(«) die asymptotische Darstellung: 


R® (a) ~ const 
Zn, m dna 


Es soll nur noch darauf hingewiesen werden, daß für sehr 
große reelle dna das Produkt gdn« sich von der Größe r 
(Entfernung Aufpunkt—Nullpunkt) beliebig wenig unterscheidet. 


1) G. N. Watson, Treatise on Theory of Bessel Functions. S. 220. 
Cambridge 1922. 
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Da ferner m(@) R,, (+, g) ein Integral der Schwingungs- 
gleichung ist, 'so un wir, daß dies für Werte von «, 
die sich sehr wenig von 7X’ unterscheiden, in den Ausdruck 
für eine divergierende Kugelwelle übergeht. Weiter sieht man 
auch leicht, daß die Integrale der zweiten Art allein diese 
Eigenschaft haben; denn die Funktionen y, (zr), mit denen die 
Funktionen der ersten Art gebildet worden sind, haben die fol- 
gende asymptotische Darstellung: 
e x 
(1 + R, 


„Hier erkennt man, daß neben der divergierenden Welle auch 
eine einfallende vorhanden ist; daraus ergibt sich dann die 
Behauptung. 

Für die re werden wir daher zwangläufig zu 


m usw. geführt: 


n m 


Diese Gleichung stellt eine Komponente des Hertzschen 
ss“ Wektors der Beugungswelle dar. Es sei gleich hier darauf 
hingewiesen, daß die Beugungswelle nicht notwendig aus einer 
einzigen Komponente des Hertzschen Vektors hergeleitet 
werden kann, auch dann nicht, wenn dies, wie wir voraus- 
setzten, bei der einfallenden Welle der Fall ist. Daher 
werden wir gewärtig sein müssen, daß wir für den Hertz- 
schen Vektor drei verschiedene Gleichungen der Form (138) 
erhalten werden, deren Koeffizienten wir dann noch durch 
ein oben angeheftetes bzw. z, y, z voneinander unterscheiden, 
je nachdem es sich um die bzw. z-, y-, z-Komponente des 
Hertzschen Vektors handelt. Ob diese Komplikation wirklich 
notwendig ist, das werden die Grenzbedingungen lehren. 

Der Ansatz (1:3) für das ursprüngliche Feld entspricht 
der Annahme, daß der Außenraum materiefrei ist, d. h. die 
Leitfähigkeit x ist gleich Null, die Dielektrizität ¢ gleich 1. 
Im Innern des beugenden Körpers sind Leitfähigkeit und 
Dielektrizität zwar konstant, doch von Null bzw. 1 verschieden. 
Man überzeugt sich leicht, daß die Komponenten des Hertz- 
schen Vektors der etwas komplizierteren Schwingungsgleichung 
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(139) AV + 


geniigen, in der also der Schwingungsparameter eine komplexe 
Größe ist. Damit wird auch der Parameter = Baie 


—Anixy 


eine komplexe Größe, und damit sind auch die Funktionen 2 a: 
$), Sn,m(#, usw., Ri) usw. komplex. Dieshat 
an sich keine große Bedeutung, denn man überzeugt sich leicht, es 
daß in den Beweisen des ersten Teiles dieser Arbeit niemals die 
Voraussetzung, die Größe w* sei reell, eingeführt worden ist. 
Nur ein Schluß ist nicht mehr stichhaltig: Wenn w*? komplex 
ist, darf man nicht mehr schließen, daß das _— 


von Null verschieden ist; es könnte auch Null sein. Dies hat 
nun eine merkwürdige Folge. Haben wir eine Funktion z. B. 


welche Weise das geschah, bleibe 


daB 


n, m -ffr- m sin #d 
- 


ist, Shem gelten die anderen Orthogonalrelationen 


N Rum‘ Ry, sndddtdg = 0 
0 0 


unverändert, doch ist der fragliche Schluß an die Rees | 


SS (Ra, „)? sin dy + 0 


gekniipft, ie diese ist nicht mehr notwendig erfüllt. Wir 
dürfen also wohl Reihenentwicklungen mit den R,,„ (‘7, p) usw. 
bei komplexem w* vornehmen, nur die Art der Koeffizienten- 


 , 
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berechnung läßt sich nicht mit Sicherheit so durchführen wie 
im reellen Fall. Ob eine allgemeine Methode zur Koeffizienten- 
berechnung angegeben werden kann, ist eine offene Frage, die 
hier nicht entschieden werden soll. 

Der Hertzsche Vektor für das Innere wird durch die 
Bedingung gekennzeichnet sein, daß er überall im Innern des 
__ beugenden Körpers endlich sein muß. Wir haben bei der Be- 
trachtung der Funktionen zweiter Art beweisen können, daß 

diese für «= 0 einen logarithmischen oder einen Pol höherer 
Ordnung haben. Daher werden Reihenentwicklungen, die 
Funktionen zweiter Art enthalten, nicht in Betracht kommen. 
Wir gelangen so zu folgendem Ansatz: 


Dabei haben wir durch die gewählte Schreibweise zum Aus- 
druck gebracht, daß die Funktionen für das Innere von dem 
komplexen Parameter ö abhängen. Übrigens müssen wir auch 
hier darauf aufmerksam machen, daß möglicherweise erst drei 
verschiedene Ausdrücke vom Typ (140) den Zustand im Innern 
bestimmen, und daß wir analog wie früher die Koeffizienten 
dreier verschiedener Vektorkomponenten durch oben angeheftetes 
bzw. x, y, z voneinander unterscheiden. 

Die Komponenten eines Vektors bleiben Translationen 
des Koordinatensystems gegenüber ungeändert. Bei Drehungen 
des Koordinatensystems transformieren sie sich wie die Koordi- 
naten selbst. Die Koeffizienten dieser Transformation sind 
also die Cosinus der neun Winkel, die die drei ursprüng- 
lichen Achsen mit den drei neuen Achsen bilden. Zu dieser 
Art von Transformationen gehört auch der Übergang von 
einem kartesischen zu einem elliptischen Koordinatensystem. 
Denn auch die elliptischen Koordinaten sind orthogonal. In 
einem bestimmten Punkt des Raumes werden die Koordinaten- 
_ achsen bestimmt durch die drei Normalenrichtungen, die in 
diesem Punkte für die drei dort sich schneidenden Flächen 

_ (Ellipsoid, einschaliges Hyperboloid, zweischaliges Hyper- 
boloid) existieren. Der cos zwischen der z-Achse und der 
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Ellipsoid hat. Wenn man in dieser Richtung um die Strecke 
ds, fortschreitet, ändert sich nur « um den Betrag da. Be- 
zeichnen wir die Wurzeln aus den Koeffizienten des Bogen- 
elementes [Gleichung (9)] mit bzw. A, B, C, so ist 

dx 1 62 

dsa A da’ 
und wir erhalten die Transformationsformeln fiir die Vektor- 


komponenten: 
1 


+%,5 
(141) P, = 5 (57 +8, $4 


= 


Bedeutet ® den Hertzschen Vektor, so BEE wir durch A: 
ein- bzw. zweimalige Anwendung der Operation rot zu den j 


Feldstärken € bzw. §. Den Ausdruck für die Rotation in 
krummlinigen Koordinaten entnehmen wir der Literatur’): 


ö é = | 
(rot P), = (CB) - 5 )} 


(142) (rot P, = aa | (C 


3 
Ersetzen wir hier die krummlinigen Komponenten durch die 


kartesischen [nach Gleichung (141)], so ergibt sich: 
Ox 


0 
— (Bede +B, + BZ 
und durch Wiederholung der Operation: 


1) Z.B.J.Spielrein, Lehrb.d. Vektorrechnung. $. 177. Stuttgart 1916. 
Annalen der a IV. Folge. 83, 


| 
| 
| 
| 
| 
“i 
|... 
| öx dy dx 
04) | a0 lay + +35) 
| ö dx dy dx 
a 


0%, Oz OB. 
83 | BO OB ay 


wo das Summenzeichen andeutet, daß noch zwei weitere Aus- 
drücke hinzukommen, die aus dem hingeschriebenen entstehen, 
wenn man überall z durch y bzw. z ersetzt. 

aa Die elektromagnetischen Grenzbedingungen fordern Stetig- 
keit der 8- und y-Komponenten von € und § an der Ober. 
Bäche des Ellipsoides. Bei der Durchführung dieses all- 
az: gemeinen Falles empfiehlt es sich allerdings, wie schon vorher 


: bs i drücken läßt, wie durch den Hertzschen Vektor. Das be- 
deutet, daß, wenn man in der Gleichung (148) ® durch € er- 
‚setzt, sich nur die Komponenten der magnetischen Feldstärke 
mit — = multiplizieren. Daher läßt sich die Stetigkeit der 


a tangentialen Komponenten von § an der Oberfläche des Ellip- 
soides so formulieren: 


> 


GR 
& 

Q Q Q Q 


4 


(1. 


1 OR. da OR, ae) Als 
ö [ B 6B. dx df, 
1 a A > dx _ (14 
afc dB, dB, dx 
| 6a #25 op” all Sel 
1; a[B OP. da 
|= 
po 
lig 
au 
Gl 
angedeutet, aus den Grenzbedingungen nicht die Komponenten u 
des Hertzschen Vektors, sondern die Komponenten z. B. der 
oh eb elektrischen Feldstärke selbst zu berechnen. Man beachte dazu, 
ag “i daB sich nach den Maxwellschen Gleichungen die magne- 
tische Feldstärke, abgesehen vom Faktor — ebenso aus- 
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Als zweite Bedingung kommt hinzu: 


SchlieBlich ist dann noch zu beachten, daB im ganzen Raume 
die Divergenz von © verschwinden muß, d.h. die Divergenz- 
bedingung gilt sicher auch an der Oberfläche des beugenden 
Ellipsoides und zwar für die innere wie für die äußere Kom- 
ponente. Dies sind zwei weitere Gleichungen. In rechtwink- 
ligen Komponenten, aber mit Differentialoperationen, die sich 
auf die Koordinaten a, y beziehen, schreiben diese 
Gleichungen so: 


(145b) 


(145 c) 


Diese Gleichungen sind keineswegs homogen, wie es zunächst 
den Anschein hat, weil sie ja auch die Komponenten des ur- 
sprünglichen elektrischen Feldes enthalten, die uns bekannt sind 
bzw. die wir aus der einen uns gegebenen Komponente (136) 
durch Ausführung der Operation rot berechnen können. Diese — 
lassen wieder eine Darstellung der Form (136) zu; wir finden 
nämlich durch die Operation rot, wenn wir (136) als die 
z-Komponente des Hertzschen Vektors ansehen: - 


€,=0 
(2 008% + ysin cosy, + ssin sin gy) 


- sin Pol Po) 


i 
Vee 
= 
Br: 
| 
| 
) 
7 
3 
on : 


F. Möglich 
Daraus folgt dann für die y- und z-Komponenten je eine der 
Gleichung (137) analoge Reihendarstellung. 
Die sechs Gleichungen 145a, b,c reichen hin, um die dreiKom- 
_ ponenten der inneren und die drei Komponenten der Beugungs- 
welle zu bestimmen; zwar nicht eindeutig, denn ein Teil der 
Beziehungen ist differentieller Art. Doch wenn wir die Be- 
_ dingung hinzufügen, daß die zu findenden Ausdrücke in ihrem 
Geltungsgebiet überall endlich sein sollen, so wird damit die 
Eindeutigkeit der Lösung gewährleistet. Diese Bedingung sagt 
nämlich aus, daß die gesuchten Funktionen Reihen sein werden, 
die im Falle der Beugungswelle die Form (138), im Falle der 
zy inneren Welle die Form (140) haben. 


Die Berechnung der unbestimmten Koeffizienten in Beugungs- 
welle und dem Zustand für das Innere, läßt sich nun dadurch 
leicht vollziehen, daß man die sechs Gleichungen (145a, be) 
nach den Funktionen eines orthogonalen Flächenfunktions- 
systems entwickelt, z. B. nach den Funktionen R, (9,9); 

des Außenraumes, weil hier 
nach Voraussetzung w® reell ist, und daher alle Sätze über 
die Koeffizientenberechnung gültig sind. Dabei stehen wir 
dann vor der Aufgabe, die Funktionen des Innern R, „(9,9; 0) 
usw. mit komplexem ©? nach den Funktionen des Außen- 
- raumes zu entwickeln. Dies aber haben wir bereits im Para- 
_ graphen 13 durchgeführt. Wenn wir dann schließlich die 

Koeffizienten gleicher Entwicklungsfunktionen einander gleich- 
setzen, so gelangen wir zu einem System unendlich vieler 


Gleichungen mit den unendlich vielen unbekannten ra, |; 


ra’; ran; 
stellte Aufgabe im Prinzip gelöst. Nur könnte sich noch ein 
Bedenken erheben: Zu unseren Grenzbedingungen gesellten 
sich die Bedingungen für das Verschwinden der Divergenz, 
die wir aber nicht für den ganzen Raum, sondern nur für die 
Oberfläche des Ellipsoides erfüllten. Es fragt sich nun, ob 
das so bestimmte Feld wirklich im ganzen Raum der Divergenz- 

_ bedingung gehorcht, was sicher unbedingt erforderlich ist. 
Sicher ist offenbar nur, daß das gefundene Feld in unendlich 
| kleiner Entfernung von der Oberfläche der Bedingung genügt 
u und damit mit dem gesuchten Felde identisch ist. Damit 


ri ; ri. Dadurch ist die ge- 
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kennen wir aber von jeder der Komponenten von & den Wert 
selbst und den Wert der normalen Ableitung auf der Ober- 
fläche des Ellipsoides. Dann aber sind wir in der Lage, mit 
Hilfe der Kirchhoff-Huygensschen Formel jede der drei 
Komponenten analytisch in den Raum hinaus so fortzusetzen, 


daß die Werte von z.B. €, und Se stetig in die von uns 


gefundenen Randwerte übergehen. Bei der analytischen Fort- 
setzung mit Hilfe der Kirchhoffschen Formel ist allerdings 
die Frage, ob nicht auch die unendlich ferne Kugel einen 
Beitrag liefert, da ja die Komponenten im Unendlichfernen 
nur wie 1/r verschwinden. Dies ist jedenfalls dann zu ver- 
neinen, wenn die Wellen divergieren'), und das tun ja gerade 
die Beugungswellen, denn die Tatsache vom Nullpunkt weg- 
eilender Wellen hatten wir durch die Benutzung der Funk- 
tionen zweiter Art zum Ausdruck gebracht. Dann aber ist 
die analytische Fortsetzung eindeutig, d.h. die auf diesem 
Wege gewonnenen Raumfunktionen stimmen mit den von uns 
aufgestellten Reihen vollkommen überein, und da dies die 
einzige Fortsetzung des Komponententripels ist, so muß dies 
auch die gesuchte Wellenbewegung darstellen, d.h. die Diver- 
genzbedingung muß im ganzen Raume erfüllt sein. 

Wir betrachten nun noch kurz den Fall, daß die Leit- 
fähigkeit des beugenden Teilchens unendlich groß ist. Die 
Grenzbedingungen hierfür lauten: Verschwinden der tangen- 
tialen Komponenten von © an der Oberfläche des beugenden 
Teiles. Aus dem Hertzschen Vektor erhielten wir durch 
zweimalige Rotationsbildung die mit o multiplizierte Feldstärke 
Andererseits gibt es eine Rechenregel, die besagt: 


rot rot® = graddiv¥ — 4¥ 


oder da „ 
4$+—F=0 

c 

2 
rot = + grad 
Wenn wir hiervon Gebrauch machen, so finden wir an Stelle 
von (144): 


1) P. Debye, Ann. d. Phys. 30. S. 755. 1909 oder Handbuch d. 
phys. Optik. Band I. S. 559ff. Leipzig 1926. : 
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Jetzt nehmen wir an, nicht nur das urspriingliche Feld, sondern 
auch die Beugungswelle lieBe sich aus einer einzigen Kompo- 
nente des Hertzschen Vektors, z. B. ®, ableiten. Dann aber 
lieBen sich die Grenzbedingungen (147) nur befriedigen, wenn 


Pa=a, =, al = 0 


ist. Dies aber ist unmöglich. Denn das elektromagnetische 
Feld ist bestimmt allein durch den Wert der Feldstärke auf 
einer geschlossenen Fläche oder durch den Wert der zu dieser 
Fläche normalen Ableitung. Beide Werte willkürlich vorzu- 
geben ist nicht möglich. Daher haben wir zu schließen, daß 
das Beugungsfeld sich nicht aus nur einer Komponente des 
Hertzschen Vektors ableiten läßt, sondern daß dazu im all- 
gemeinen alle drei Komponenten nötig sind. Dann aber finden 
wir als hinreichende Bedingung für das Verschwinden von €, 


und €: 
(148) - 


Um diese Gleichungen aufzulösen, verfahren wir analog dem 
allgemeinen Falle. Für den Hertzschen Vektor der Beugungs- 
welle machen wir den Ansatz (138) und zwar für jede Kom- 
ponente einen. Wieder entwickeln wir dann jede Gleichung 
nach einem System orthogonaler Funktionen und vergleichen 
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die Koeffizienten gleicher Entwicklungsfunktionen. So erhalten 
wir wiederum ein System unendlich vieler Gleichungen zur 
Bestimmung der unbekannten Koeffizienten. 

Hiermit ist das Beugungsproblem vollständig gelöst. An 
die Gleichungen (148) bzw. (145a, b,c) haben dann die mehr 
oder minder komplizierten Auswertungen anzuknüpfen. N} 


$ 17. Ein numerisches Beispiel 


Wir hatten im § 14 ausführlich den Fall des gestreckten Pee 
Rotationsellipsoides besprochen und ganz kurz gezeigt, wie auf 
diesen der Fall des abgeplatteten Rotationsellipsoides zuriick- _ 
geführt werden kann. In diesem Paragraphen soll als Beispiel §= 
fir die Auswertung der physikalischen Grenzbedingungen die aX 
Beugung an einem speziellen abgeplatteten Rotationsellipsoid 
durchgeführt werden, weil unter der Schar solcher Ellipsoide ein 
besonders einfaches ist, nämlich die Kreisscheibe. Die Beugung 
an einer Kreisscheibe soll hier untersucht werden. Dies ist 
eine Aufgabe, die bereits Lommel?) mit den Hilfsmitteln der _ aa 
Fresnel-Kirchhoffschen Beugungstheorie durchführte, eine a ie 
Aufgabe, die sich dort nur mit einem sehr umfangreichen 
mathematischen Apparat erledigen ließ. 

In der zy-Ebene sei eine Schar konfokaler Ellipsen ge- 
zeichnet, die Brennlinie liege auf der y-Achse, ihr Mittelpunkt 
sei der Nullpunkt des Koordinatensystems, der Abstand der 
beiden Brennpunkte sei 2g. Durch Rotation dieser Figur um 
die z-Achse entsteht dann eine Schar abgeplatteter Rotations- 
ellipsoide, unter denen auch die Kreisscheibe ist. Diese ent- 
steht durch die Rotation der Strecke zwischen den beiden 
Brennpunkten um ihr Mittellot. Ihr Durchmesser ist 2g. An 
ihrer Stelle befinde sich eine unendlich diinne, unendlich gut 
leitende Metallfolie der gleichen Größe; die Beugungserschei- 
nung an dieser Metallfolie ist zu untersuchen. Man sieht 
sofort, daß ihre Gleichung durch Te se 
gegeben ist. ae 

Die einfallende Welle soll eben sein und in der Richtung 
der z-Achse, von positiven zu negativen z fortschreiten. Die © 
Welle sei linear polarisiert, der magnetische Vektor schwinge 


1) E. vy. Lommel, Abh. d. bayer. Akad. d. W. 15. S. 531. 1886. 


Pin, 


r 


J 


‘dn 
Pie 
A 
4 


Möglich 


in der Richtung der z-Achse. Der Hertzsche Vektor, de 
diese Welle darstellt, ist dann: A. 


(150) P%=0; P= 3; 


i+ — 
e 


Unter Benutzung der Gleichungen (116a) auf S. 710 und (149) 
schreiben sich die Grenzbedingungen (148) in folgender Form: 


B,), cos p +(P.),sing =0, F 
4 Vine sin g— (B.' cos 0, 


(P,), cos cos + (P,), cos Fsin p + (Fe). sin = 0. 


Aus den ersten beiden Gleichungen ließt man sofort ab, dab ®, 
und ®, fir 7=0 zu verschwinden haben. Da es aber nach 
dem Ansatz (150) nur für ®, im Außenraum eine Anregungs- 
schwingung gibt, so ist zu schließen, daß ®, im ganzen Raume 
identisch verschwinden muß. Ebenso ergibt die dritte Gleichung 
das identische Verschwinden von ®,, so daß wir in unserem 
Falle nur eine nicht verschwindende Beugungskomponente oo 
haben. Diese, zusammen mit der Komponente der einfallenden 
hat der Bedingung = 0 zu genügen. 
Die Beugungserscheinungen treten besonders deutlich auf, 
er wenn die Wellenlänge der einfallenden Strahlung mit den 
_ Lineardimensionen des beugenden Körpers vergleichbar wird. 
1% Im Rahmen dieser Bedingung haben wir mit der Auswahl der 
Wellenlänge noch ziemliche Freiheit und werden daher die 


. as Wellenlänge so wählen, daß die Eigenfunktionen pe” (2), 
ae nach denen wir die Reihenentwicklungen vorzunehmen haben, 
is) sich möglichst einfach berechnen lassen. Wir setzen daher 


fest, die Wellenlänge sei 7 mal dem Durchmesser der beugenden 
Kreisscheibe. Dies hat zur Folge, daß o=—-g=1 wird. 


5 Wir hatten auf 3. 710 gezeigt, wie der Fall des gestreckten 
Me und des abgeplatteten Rotationsellipsoides formal durch die- 
selben Gleichungen behandelt werden können, wenn man den 
Koordinaten auch komplexe Intervalle gestattet. Dann hat 
af man allerdings die Größe g durch iq zu ersetzen. In unserem 
Falle wird daber der Parameter © =i. Die mit diesem be- 
__-rechneten Eigenfunktionen sind dann: 


Be 
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pe(z)= 0,707 + 0,081P,(~) + 0,001 P, (x)... 


152 | e,(2) = 1,281 P,@) + 0,051 P,(a).. 
( pe, =—0,036 + 1,579P,(~) + 0,036P, (x) + 0,0047P, (a)... 
pes(z)= -0,027P,(@) + 1,582P,(z) + 0,028P,(x)... 


Die vorstehenden Reihen sind so hingeschrieben, daB sie Er- 
gebnisse liefern, die auf zwei Stellen hinter dem Komma genau 
sind. Es ist dies die Genauigkeitsgrenze, die wir im folgenden 
einhalten wollen. Wir erhöhen die zweite Stelle, wenn in der 
dritten eine Zahl steht, die größer oder gleich fünf ist, im 
anderen Falle bleibt die zweite Stelle ungeändert. 

Wir haben jetzt die Entwicklung der einfallenden Welle 
nach den Eigenfunktionen vorzunehmen. Dabei kommt uns 
nun sehr zu statten, daß der Ausdruck für die einfallende 
Welle, abgesehen von einem konstanten Faktor, vollkommen 
mit dem Kern Gleichung (117) übereinstimmt, so daß die 
Reihenentwicklung nach den Eigenfunktionen nichts anderes 
ist, als die bilineare Entwicklung eines Kernes nach seinen 
Eigenfunktionen. Bezeichnen wir jetzt die Eigenwerte des 
Kernes e’®=y mit e,, so ergibt sich: ED 


Aus” ~~ Integralgleichung (121) rechnet man leicht mit — 
Hilfe der Reihenentwicklung für durch gliedweise 


Integration nach, daß wir aus pe, (zx) gerade erhalten, 


wenn wir die Kugelfunktionen ?,,(z) mit geradem Index durch 
2-w,,(@y) und die Kugelfunktionen mit ungeradem Index ~~ 
durch 27-y,,,,(@y) ersetzen. Bezeichnen wir diese Reihen eee : 

mit pe,(y), so erhalten wir die Reihenentwicklung: is pes 


(153) = Spe, (y)pe, (2). 


Das Argument der Zylinderfunktionen y, ist, wie man aus 
den Bemerkungen von S. 710 und unter Beachtung der Glei- 
chung w =i sofort erkennt, die Größe Sinn, die fir 7 = 0, 
also an der Oberfläche der beugenden Folie, verschwindet. Da die . 
Zylinderfunktionen yw, Potenzreihen sind, deren kleinste Potenz = 
gleich n ist, so verschwinden die Funktionen P 2n41(0), die 
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einen ungeraden Index haben, während die Funktionen pe,, 
mit geradem Index sich für 7 = 0 auf 24,2» reduzieren, wo 
4A," das konstante Glied in der Kugelfunktionsreihe der 
Funktion pe,(z) bedeutet, also z. B. 4,?= — 0,036 gemäß Glei- 
chung (152). Es zeigt sich nun, daß bereits 4, kleiner als 
4 10” ist, ferner, daß pe,(z) kleiner als drei bleibt, solange z 
_ -reell ist und nicht größer wird als 1. Diese Bedingung ist 
also erfüllt für eine Winkelfunktion, denn in unserem Falle 
ist z durch sin % zu ersetzen. In der von uns angestrebten 
_  Näherung ist also an der Oberfläche der Metallfolie bereits 
das Glied pe,(y)pe,(z) zu vernachlässigen. Da die ungeraden 


(134 )o= — ig {1,414 pe, (sin — 0,072 pe, (sin #)}. 


Oberfläche der Folie. 
ie Um nun der Grenzbedingung (P,), = 0 zu genügen, haben 
- wir eine Beugungswelle zu konstruieren, die an der Oberfläche 
der Folie gerade dem Negativen von (154) gleich wird. Dies 
läßt sich immer erreichen durch einen zweigliedrigen Beugungs- 
satz, nämlich: 
(155) = B, 7 ¢, (Sin 1) pe, (sin F) + B, ge, (Gin n)pe, (sin 9). 


- Man sieht dann leicht, daß der Bedingung (P,), = 0 genügt 

Das einzige, was noch Schwierigkeiten bereiten könnte, 
ist die Berechnung des Wertes ge,, (0). Wir haben zwar Reihen 
für ge,,(y). Die eine entsteht aus der Reihe pe,,(y), wenn 
wir ,, durch %,, ersetzen und die Funktion mit dem Eigen- 
_ wert e,, multiplizieren. Doch hört diese Reihe bereits zu 
pe auf, wenn y = 1 wird, ist also sicher nicht mehr 
für 7 = 0 zu gebrauchen. Aber ebenso, wie wir für die Funk- 
tionen erster Art pe,(z) verschiedene Entwicklungen nach den 
Zylinderfunktionen erhielten, einmal nach y,(oz), dann aber 
: ; auch nach w,(wiY1— z,), so ergibt sich eine analoge Reihe 


er 
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wo Ae”) die Koeffizienten der Funktion pe,, (x) in der Kugel- 

fiaktionsreihe (153) bedeuten. Unter j,, seien die Higenwerte 
des Kernes J, (iw V1 — z? — y? verstanden. ‘Die Reihe (157) 
konvergiert ebenfalls, solange das Argument der Zylinderfunk- 
tionen größer als 1 ist. In unserem Falle ist das Argument 
Goj7. Daher konvergieren die Reihen (157) bis an die Folie 
heran, aber nicht auf der Folie. Die Reihen (157) werden wir 
also immer anwenden müssen, wenn wir das Feld in der Nähe 
des beugenden Schirmes berechnen wollen, über die Werte 
ge,„(0) erhalten wir aus ihnen jedoch keinen Aufschluß. Daher 
haben wir auf die Integraldarstellung (133) zurückzugehen, in 
der allerdings noch die unbestimmt gelassene Konstante zu 
berechnen ist. Es zeigt sich, daß diese Darstellung mit der 
Reihe (157) in ihrem Konvergenzbereich identisch wird, wenn 
wir die Konstante gleich 2(e,,)? wählen, wo e,, wie vorher 
definiert, ein Eigenwert des Kernes e’ = ist. Wir erhalten so: 


+1 h 
-1 
Hier beachte man nun, daß YW, (x) = —_— ist. So ergibt sich 
schließlich: 
(159) q¢,, (0) =—2—" ~p ey, (2)dz. be | 


Hier setze man nun für die e-Funktion ihre Potenzreihe ein, Ex 


verschwinden, weil die pe,,(z) mit geradem Index gerade a ie ft 
Funktionen sind. Die e-Reihe breche man dann an einer Stelle 

ab, an der der Fehler kleiner wird als 5-10-*. Es zeigt 
sich, daß folgende Integral hinreichende 
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x 
(1598) ge,,(0) =— 2ike, (t+ art 


Dann ersetze man nur noch pe,,(z) durch seine Reihe (153) 
e bh und erhält dann eine Reihe von Integralen der Form: 
z zZ. 
Diese ii aber in der Gleichung (125) bereits ausgewertet. 
Auf diese Weise erhält man dann aus (156): 
(160) B, = 9:21,758; B, = — 9:0,001. 
Der Hertzsche Vektor für das gesamte elektromagnetische 
Feld ist dann: 
| P, =— Sn _ 21,753 ge, (Sin 7) p e,(sin #) 
= + 0,001 ge, (Sin 7) pe, (sin 9)}. 
Mit Hilfe der Gleichung (143) läßt sich dann die magnetische 
Feldstärke berechnen, und das Quadrat des absoluten Betrages 
dieser ergibt die Intensität. 

Führen wir die Rechnung durch, so zeigt sich eine be- 
merkenswerte Schwierigkeit. Der Ausdruck für die Feldstärken 
nimmt auf der z-Achse und direkt an der beugenden Folie 
den Wert oo an. Dies liegt daran, daß die Reihen nicht 
überall gliedweise differentiiert werden dürfen, weil auf ge- 
ae wissen Linien die Entwicklungsfunktion alle unendlich groß 

werden. Man könnte z. B. daran denken, mit Hilfe der Green- 
schen Formel den Ausdruck für die Feldstärken, der außerhalb 
der in Frage stehenden Geraden überall gültig ist, in diese 
hinein fortzusetzen. Doch ist es noch einfacher, nicht die 
magnetische, sondern die elektrische Feldstirke') zu berechnen. 
Diese ergibt sich aus ® durch zweimalige Anwendung der 
Operation „rot“. Da aber 

rotrot $ = grad div® — AB 

ist und ® der Schwingungsgleichung genügt, so ergibt sich: 


| E- RP + grad div®. 
en Der Vektor $ war so bestimmt, daß die Tangentialkompo- 


sup: 1) Diese ist natürlich nicht gleich der magnetischen Feldstärke, 
EN N da in der Nähe der Metallfolie stehende Wellen zu erwarten sind. 
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nenten von © an der Oberfläche der beugenden Scheibe Null 
werden. ‘§ selber ist aber ebenfalls dort Null. Daraus folgt, 


y 


Hertzsche Vektor für die einfallende ebene Welle divergenz- — 
frei, daher entsteht 8 allein durch Differentiation des aes j 
Teiles von ®, der vom Nullpunkt wegeilende Wellen darstellt, ist 
also wiederum eine solche und gestattet daher die Darstellung © 

= >> ge" (Sinn)pe” (sin #) cos m(p — 


Da nun Eu Funktionen ge” (0), wie ein Blick auf die a 
Gleichung (133) lehrt, niemals Null sind, müssen die Koeffi- 
. 
Oy ist 
also im ganzen Raume identisch Null. Ebenso zeigt man auc 
daß oat 
von Null verschieden, weil der Bedingung divE =0 genügt 
werden muß. Diese Komponente hätten wir durch Differen- 
tiation zu berechnen und würden finden, daß der analytische 
Ausdruck auf der x-Achse zu gelten aufhört, weil wieder dort _ 
die Entwicklungsfunktionen unendlich groß werden. Doch ist | 
der Wert von €, auf der z-Achse leicht anzugeben; denn aus 
Symmetriegründen muß für y=z=0 r 
der Energiestrom die Richtung der 
z-Achse haben, der elektrische und : 
magnetische Vektor in einer Ebene + sine 
senkrecht zur z-Achse liegen, woraus ’ 
sofort folgt: 
(€,)y=2=0 =. Z 
Auf der z-Achse berechnet sich daher ae 
die Intensität allein aus der y-Kom- 
ponente des elektrischen Vektors; 6 : 
diese stimmt mit $8, durchaus überein und wird durch Glei- 
chung (161) dargestellt. 
Für g=1 haben wir die Intensitäten auf der z-Achse ~ 
unmittelbar vor und hinter dem beugenden Schirm berechnet a, ; 
und in einer graphischen Zeichnung veranschaulicht. Die 
Kurve I stellt die Intensität vor die Kurve II hinter = 


h 


dort verschwinden muß. Nun ist aber dr 


zienten der Reihe verschwinden. Die Komponente 


identisch verschwinden muß, jedoch ist 


“= 
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vive Schirm dar. Die Abszisse ist die Entfernung des Punktes von 
Bu u der beugenden Scheibe, die Ordinate stellt die Intensität in 

BER er Prozenten der einfallenden Intensität dar; die Wellenlänge 
ist 2a. Die beiden Kurven berühren sich bei x == und dies 
ist nur natürlich, da der Ausdruck für R, vor und hinter dem 
_ $Schirm auf der z-Achse sich nur durch verschiedenes Vor- 
zeichen im Exponenten der e-Funktion unterscheidet. Er liefert 
also die gleichen Werte dort, wo der imaginäre Teil der 
e-Funktion verschwindet, und dies tritt bei jedem ganzen Viel- 
fachen von ein. 


$ ei F. Möglich. Beugungserscheinungen an Körpern usw. 


Zusammenfassung des physikalischen Teiles 
Seer Im § 15 wird gezeigt, wie sich aus der Gestalt des 
beugenden Ellipsoides das Koordinatensystem herleiten läßt 
N a und welcher Gleichung das beugende Ellipsoid in diesem Ko- 
ordinatensystem genügt. Der Hertzsche Vektor als Integral 
3 i der Maxwellschen Gleichungen wird eingeführt und das all- 
_ gemeine Integral der Wellengleichung angegeben, auf elliptische 
Koordinaten umgerechnet, und die ihm äquivalente Reihen- 
entwicklung gegeben [Gleichung (186—137)]. Der nächste 
Paragraph zeigt, wie die elektromagnetischen Grenzbedingungen 
benutzt werden können, um ein spezielles elektromagnetisches 
Feld wirklich zu berechnen. Der § 17 enthält die Illustration 
der Methode an einem Beispiel. 


2 Zum Schluß möchte ich nicht verfehlen, meinem hoch- 
verehrten Lehrer, Herrn Prof. Dr. M. von Laue, meinen er- 
= gebenen Dank auszusprechen für das ständige Interesse, das 
er dem Fortschreiten der Arbeit entgegen gebracht hat, sowie 
für manchen wertvollen Rat und Aufmunterung bei den oft 
J % schwierigen mathematischen Untersuchungen. 
Dank schulde ich ferner dem Kaiser Wilhelm-Institut fir 
Physik, das mich durch Gewährung von Geldmitteln in die 
a Lage versetzte, numerische Rechnungen von anderer Seite aus- 
fahren zu lassen. 
Fräulein L. Wohlleben unterstützte mich in freundlicher 
Weise bei der Auswertung der Integralgleichungen und bei 
anderen orientierenden Rechnungen. EN 


(Eingegangen 16. Mai 1927) 
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Uber die 
der Brownschen Bewegung durch das Licht; 


von W. Pospisii 


I. Experimenteller Teil 


1. Der Zweck der vorliegenden Arbeit war durch Mes- 
sungen der Beweglichkeit mikroskopischer Teilchen bei ver- 
schieden intensiver Beleuchtung quantitativ zu konstatieren, ob 
und und in welchem Maße das benutzte Licht diese Bewegung 
beeinflussen kann. 


2. Übersicht der älteren, ähnliches Ziel verfolgenden 
Arbeiten.!) Regnauld (1858) beobachtete auf kugelförmigen 
Partikeln von 0,25 Durchmesser, daß ihre Bewegung durch 
Verschwächung des einfallenden Lichtes verkleinert werde, 
und er findet deswegen die Ursache der Brownschen Be- 
wegung in den Flüssigkeitsströmungen, welche durch die 
Wärme des von den Partikeln absorbierten Lichtes hervor- 
gerufen werden. Zu derselben Anschauung ist 8. Exner (1867) 
geführt worden. Exner hatte die Bewegung verschiedener 
Substanzen mittels eines Gitters im Mikroskopokular gemessen 
und fand die Bewegung bedeutend größer bei einer inten- 
siveren Beleuchtung. Beide Arbeiten leiden an dem Umstand, 


daß die mit der Temperatur veränderliche Viskosität der — 


Suspension nicht berücksichtigt wird. 


Nachdem schon von Carbonelle (1874) die kinetische _ 


Ursache nach der Gastheorie ausgesprochen wurde, untersucht 
Gouy (1888) qualitativ verschiedene äußere und innere Ein- 
flüsse auf die Brownsche Bewegung und berichtet, daß er 


bei einer Lichtintensitätsänderung im Verhältnis 1:1000 keine _ 


Änderung in der Größe der Bewegung wahrnehmen konnte. 
Gouy sagt: „La lumiére ne joue done aucun role sensible 


1) Bezüglich der Literaturangaben sei auf die Monographie de Haas- 
Lorentz, „Die Brownsche Bewegung“, Sammlung Wissenschaft, ver- 
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dans le mouvement brownien“ (C.R. 109. p. 104. 1889), 
Gouy hatte aber die Bewegung nicht gemessen und suchte 
vielmehr die Nichtexistenz und auffälligere Änderungen der 
Bewegung. 

Auch Meade Bache (1894) hatte keine Einwirkung 
des Lichtes auf die Brownsche Bewegung gefunden. 

Trotzdem referiert Kolätek (1889) über die Gouysche 
Arbeit folgendermaßen: „Referent meint, hierin (d.h. in der 
 Brownschen Bewegung) nichtsdestoweniger eine Art radio- 
metrischen Wirkung zu erblicken, zu deren Entfernung bei 
der äußersten Feinheit und Beweglichkeit der Teilchen Be- 
lichtungen gehören würden, die eine Beobachtung der Er- 
scheinung ausschließen“. (Beiblätter 13. 8. 877. 1889). 

Quincke (1898) sieht die Ursache der Bewegung in der 
Beleuchtung, indem die das Partikel einhüllende Schicht 
ungleichmäßig erwärmt wird. 

Aus der Zeit nach den entscheidenden theoretischen 
Arbeiten von Einstein (1905) und Smoluchowski (1906) 
sei erwähnt die auch nur qualitative Untersuchung Sved- 
bergs („Existenz der Moleküle“ §.98) über den Einfluß 
der «-,ß- und kurzwelligen Strahlen auf die Bewegung der 
Kolloidpartikelchen, mit negativem Erfolg, und daß V.Henri 
(1908) bei einer sehr intensiven Beleuchtung, welche zu den 
kinematographischen Aufnahmen nötig wurde, das mittlere 
Quadrat der Verschiebungen der Partikel viermal größer 
gefunden hatte, als nach der Einsteinschen Formel zu er- 
warten war. 

Nach der glänzenden Bestätigung der Einsteinschen 
Formel für das mittlere Quadrat der Verschiebungen 4,? der 
Partikeln in einer fest gewählten Richtung x 

RTt 

(I) = 

wo R die Gaskonstante, T absolute Temperatur, ¢ die Zeit, 
N die Loschmidtsche Zahl, £ die Viskosität der Flüssigkeit, 
r den Radius der kugelförmigen Partikeln bedeutet, durch 
die Arbeiten von Perrin und Svedberg ist zu erwarten, 
daß eine intensive Beleuchtung die Brownsche Bewegung 
nur in dem Maße vergrößern kann, soweit sie die Temperatur 
der erhöht und ihre innere 
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3. Die Arbeitsmethode des Autors. Als MaB der Beweg- 
lichkeit der Brownschen Teilchen diente nach der Einstein- 
schen Formel (1) das mittlere Quadrat der Verschiebungen. 
Die Positionen der Teilchen wurde in Zeitintervallen von 
10 Sek. mittels des Abbeschen Zeichenapparates auf das 
Millimeterpapier aufgezeichnet, dann wurden die Verschie- 
bungen A, und A, in zwei zueinander senkrechten Rich- 


tungen abgelesen, quadriert und das Mittel von ihnen 4,,,? 
berechnet. Die beiden Richtungen x, y sind nur wegen der 
besseren Kontrolle und Verläßlichkeit berücksichtigt worden. 
Dies geschah für eine sehr schwache Beleuchtung, welche 
noch eine verläßliche Beobachtung zuließ und für eine inten- 
sive Beleuchtung. Die dazu gehörigen Verschiebungen be- 
zeichnen wir 4,2 und 4,2. 

Die gemessenen Verschiebungen wurden dann umge- 
rechnet auf dieselbe Temperatur (293° absol.) und auf die 
zu dieser gehörige Viskosität nach der Formel 

— 298 298 &; 

(wo Ty, Co, T, und Z, die gemessenen Temperaturen und Vis- 
kositäten bedeuten), und erst dann sind diese reduzierten 
Mittel verglichen worden. Die Gefahr einer zu großen Extra- 
polation droht hier nicht, weil es sich nur um wenige Zehntel 
von Grad handelt. 


Etwaigen Überschuß von 4,24 über 4,2,.4, der eben 
nach dem oben gesagten nicht zu erwarten ist, nennen wir 
die Vergrößerung der Brownschen Bewegung durch das Licht 
und es ist dann eine Erklärung für ihn zu suchen. 


4. Die Apparatur. Nach längeren mißerfolgten Messungen 
an Toluolemulsionen im Wasser verfertigte ich eine für meinen 
Zweck sich gut bewährende Beobachtungsküvette mit den 
Dimensionen 1 x 1 x 0,1 mm, welche eine Menge von etwas 
über 0,1 mg Emulsion und ein Thermoelement aufnehmen 
konnte.!) (Fig. 1.) Die sonstige Apparatur zeigt die Fig. 2. 


bzw. Arca = 


1) Die Kiivette verfertigt man auf folgende Weise. Auf das Unter- hire Be; 
lageglas wird mit dem Balsam ein glattes, dünnes Deckgläschen auf- Kegel v A 
geklebt und nach dem Eintrocknen wird dieses mit dem Diamanten zwei- ees ae 
mal länglich in einer Entfernung von 1 mm eingeritzt, mittels Xylol frei- er 
gemacht und in drei Teile gebrochen. Der mittlere enge Streifen wird = 


Annalen der Physik, IV. Folge. 83, oe Che 
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W. Pospisil 


Als Lichtquelle diente bei den endgiiltigen Messungen eine 
Nernstsche Glihlampe, weil sich die Bogenlampe wegen 
ihrer unvollkommenen Stabilität zu diesem Zwecke nicht 


gut eignete. 


4 


\ 
Paraffin‘ Ölimmersion \ 


\£mulsion mit Thermoel. 
! 


Fig. 1 


rie 


Zeichen- Wasserbad mit 
apparat Thermometer 
Berußte Glasplatte Wasserfilter 


‘ N. 


— — 


Fig. 2 wee) 


Es wurde mit einer Olimmersion und 3560maligen 
gréBerung gearbeitet, welche einer Bildentfernung von 47 cm 
entsprach. Die Größe der suspendierten Teilchen von Mastiz, 


in der Halfte gebrochen und alle vier Stiicke werden wieder auf ein Glas, 
wie die Fig. 1 zeigt, so aufgeklebt, daß sich in den einen Ritz noch die 
dünnen Thermoelementdrähte einbetten lassen. Das neue Deckgläschen 
über der Emulsion wird mit Paraffin abgedichtet. 
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Gummiguit und Ruß war nach der Art der endgültigen Mes- 
sungen nicht unbedingt notwendig zu kennen; sie war geschätzt, 
mittels des Zeichenapparates und mehrerer auf ein Papier 
gezeichneten Kreise von 1, 2, 3, ... mm Radius, welche einer 
wirklichen Größe von 0,28, 0,56, 0,84, ... u entsprachen, 


Die Temperaturänderungen in der Emulsion konnten 
mit dem Thermoelemente Konstantan-Nickelin bis auf Tausend- 
stel von Grad bestimmt werden, aber es genügte vollkommen, 
ihren Wert nur auf Zehntel von Grad zu messen. Die Vis- 
kosität der benutzten Flüssigkeit und ihre Veränderlichkeit 
mit der Temperatur wurde mit einem Kapillarviskosimeter 
durchgemessen. 


5. Die Hauptmessungen und ihre Resultate. Die ersten 
einigen Tausende von Punktierungen, welche bei verschiedenen 
Lichtquellen und Intensitäten an Toluol-, Mastix-, Gummi- 
gutt-, Tusche-, Zimmerstaub und Rußsuspensionen im Wasser 
vorgenommen wurden, brachten nicht genug sichere und ein- 
wandfreie Resultate, sie waren nur eine notwendige Schulung 


in der Handhabung der endlichen Messungen. Aus diesen sei a 


nur angeführt: 


a) die Messung der Bewegung der Gummigutteilchen im | 


Lichte, welches sie nicht absorbieren. Es wurde an einer wäß- 
rigen Gummiguttemulsion mit Teilchenradius von 0,22 bis 
0,35 u eine Reihe von 500 einzelnen Punktierungen durch- 
geführt. Eine Hälfte von ihnen bei einer intensiven Be- 
leuchtung mittels der Nernstschen Glühlampe und die zweite 
Hälfte bei demselben Licht, welches noch durch eine 8cm 


dicke Schicht von starker alkoholischer Gummiguttlösung a 


filtriert wurde, welche alles Licht, die rechte Seite des op- 
tischen Spektrums, das von den Teilchen absorbiert werden 
konnte, entfernte. Die Temperaturunterschiede in der Emulsion 
betrugen 0,3—1° C. 


Nach der Reduktion folgt für das Verhältnis der Ver- = 


schiebungen 
A*sitr, Licht : A?weißes Lient = 100 : 107,6 


b) Messungen an denselben Kohlenteilchen. Die Orien- 
tierungsmessungen an Kohlenteilchen haben eine größere 
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W. Pospisil 


_ Geschwindigkeitssteigerung der Bewegung mit der Licht- 
 intensitaitszunahme gezeigt.!) 
Weil es sich um nichtkugelférmige Teilchen handelt, fiir 


Formel zum Stokesschen Fallgesetz) nicht kennen kann 
fa ee damit auch der Fehler aus der Wahl von nicht gleich- 
aus 4 groBen Teilchen bei verschiedener Beleuchtung ausgeschieden 

wire, wurden die Messungen auf folgende Art durchgeführt. 
Das Sehfeld wird intensiv beleuchtet, dann wird in den 


stellt (Fig. 2), welche nur soviel Licht durchläßt, damit noch 
die Teilchen beobachtet werden können. In diesem Lichte 
machte ich an einem Teilchen 10—60 Punktierungen, und 
ohne das Teilchen zu verlieren, entfernte ich die RuBplatte, 

and an denselben Teilchen zeichnete ich eine beiläufig gleich- 
große Anzahl von seinen Positionen ein. Aus beiden diesen 
_ Reihen werden nur soviel erste Punktierungen in Rechnung 
gezogen, wieviel die kürzere von ihnen enthält. 

E Bei dieser Methode fallen viele Punktierungen leer aus, 
hauptsächlich wegen des vorzeitigen Verlustes des Teilchens, 
Weil das Wasser infolge seiner kleinen Viskosität eine zu 
lebhafte Bewegung erlaubt, wurde als Suspensionsflüssigkeit 
eine Mischung von Wasser und Glyzerin im Volumenverhältnis 
2:1 benutzt, welche bei 20°C das spezifische Gewicht 1,088 hat. 


Auf diese Art wurden an 13 Rußpartikeln zwei Reihen 
von Punktierungen, jede mit 250 zählbaren Punktierungen, 
gewonnen. Die Temperaturunterschiede betrugen 0,7° C. Nach 
der Reduktion bekommt man 


A,?: 4,2 = 100:147. 


: c) Vergleichende Messungen. Mit der sub b) beschrie- 
 benen Methode wurden über Zwölfhundert zählbaren Punk- 
tierungen an Mastix-, Gummigutt- und Rußteilchen durch- 


1) Die Kohlenteilchensuspension wurde nach verschiedenen miß- 
erfolgten Versuchen auf folgende einfache Weise vorbereitet. Der käuf- 
liche Ruß wird mit Wasser geschüttelt und nach seinem Herabsetzen 


ae a Rußteilchen verschiedener Form und Größe. 


zeigt das Mikroskop in dem obenstehenden klaren Wasser pi 


740 
a welche man die veränderliche Widerstandsformel (analoge 
° 
Dtraniengang des Lichtes eine stark berußte Glasplatt Je- 
> 
Ihre Viskosität beträgt bei den Temperaturen 18, 19, 20, 21°C 
die absoluten Werte: 0,0286, 0,0276, 0,0266, 0,0256 
5 
8 
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geführt. Es wurden immer gleichzeitig in derselben Emulsion > 
Mastix- und RuBteilchen und Gummigutt- und Rußteilchen 
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suspendiert. Die Bewegung der einzelnen Substanzen wurde 
dann abwechselnd immer auf denselben Partikeln gemessen. 

Ein Beispiel aus diesen Messungen bringt die Fig. 3 
welche die 128. Protokollartabelle reproduziert. Sie gibt die 


Wer. sich die 4 Kurventeile A zur sehr schwachen und die 
Kurventeile B zu einer intensiven (2,2 cal/em? min) Be- 
 Jeuehtung.!) 
Bat . Die zahlenmäßigen Angaben dieser Messungen enthält 
die Tab. 1. 

; ; Der Mastix zeigt bei der intensiven Beleuchtung das 
mittlere Quadrat um 2 Proz. kleiner. Diesen Unterschied 
kann man wohl als eine statistische Schwankung ansehen. 
Das Gummigutt zeigt eine ähnliche Vergrößerung wie sub a). 
Wenn wir bei den Kohlenteilchen auch das Fesultat sub b) 
(mit Rücksicht auf seine größere Reihe von Punktierungen) 
in Rechnung nehmen, so bekommen wir eine Vergrößerung 
um 45 Proz. 

Das Resultat der Messung ist also folgendes: 

1. Die Brownsche Bewegung der Mastixkügelchen bleibt 
_ riicksichtslos auf die Intensität der Beleuchtung in den sta- 
tistischen Schwankungsgrenzen konstant. 

2. Das Gummigutt vergrößerte sein mittleres Quadrat 
der Verschiebungen durch die benutzte intensive Beleuchtung 

von 100 auf 107 Proz. und 

3 Fi = 8. Die Kohlenteilchen von 100 auf 145 Proz. 


a’ 


6. Die Messung der Intensität der Beleuchtung wurde 
auf folgende Weise durchgeführt. Anstatt des Präparates wurde 
auf das Mikroskopkondensor ein Diaphragma von 2,24 qmm 
Öffnung so eingestellt, daß es den intensivsten Lichtstrom 
_ durehlie8. Darüber war ein Silberscheibehen untergebracht. 


ami .. Dasselbe war r auf der unteren Seite leicht beruBt und es wurde 


1) Wie man sieht, wurde die Punktierung einigemal (auf 2—3 Sek.) 
unterbrochen und das Zeichenpapier verschoben, damit nicht zu viele 
Aufzeichnungen auf eine kleine Fläche kämen. 


; 

Positionen eines Rußteilchens von äquivalentem Halbmesser 
0,4 u in Zeitintervallen von 10 Sek. wieder, und zwar beziehen 

a die obere Seite des Scheibchens angelötet waren. Das Scheib- 

hen wurde auf 2 Min. der Einwirkung des Lichtstromes be- 


743 


ößerung der Brownschen Bewegung durch Licht 


Vergr 


| | | 9'6r | Ist / 
os‘s | | 02208 | 81200) 181 
| | 01984 | 1/08 | OFZ / 
901 : 001 | OTT : 001 | 
005 | 9'308 | 0298F a120'0| / L 
| 
Li‘g | | 01098 | 033 | / _ L 
86 : 00T | TOT : 001 | | co xusem 
ri‘s | | 00808 | 8820/00 | | 988 | _ | / L 
dan 
onynpey uonyyn uezue [2 
youu JOA 008/, u qSUION |. w0,_01 
yeyıs | amyeı | -[yRZ UG zugjsqug 
-oyst, |-odmoy) | 148g 
“Se ıqez | Zunjqonojog 


I 


STR 

= 
= 
. 
: 
2 
| 
% 
| 
ur 
4 
i 


W. Pospisil 


Er lassen und aus der graphischen Darstellung des Galvano- 
meterausschlages (Fig. 4, Kurve A) wurden dann die Wärme- 
 verluste abgemessen. So bekommt man einen linearen Zu- 
+ wachs der Temperatur des Scheibehens während der Be- 
leuchtung (Fig. 4, Kurve B).!) 

Für die mittlere Intensität I des Lichtstromes wurde 
aus sieben Messungen 
I= 2,15 + 0,15 cal/qem min = 2,2 eal/qem min 


150 - 
AG 


70H 


7408 6=7°C 


Fig. 4 


gefunden. Von dieser strahlenden Energie kann man voraus- 
setzen, daß sie ausschließlich den optischen Wellenlängen an- 
gehört, denn die anderen Wellenlängen wurden von dem Wasser- 
filter und zahlreichen Glasschichten der Apparatur beseitigt.?) 

Die bei den Punktierungen benutzte sehr schwache Be- 
leuchtung konnte mit dieser Methode nicht gemessen werden. 


1) Bezieht sich zu einem Scheibchen von 1,011 g Gewicht. 

Mr 2) Mit derselben Methode wurde konstatiert, daß nur das Wasser- 
filter (10cm Wasserschicht) aus dem gesamten Energiestrom 79 Proz. 

absorbiert. 


pt 


\ 
> 
40. 
/ 
. 


Vergrößerung der Brownschen Bewegung durch Licht 745 


Weil der absoluten Größe dieser Beleuchtung keine wesent- 
liche Bedeutung zukommt, wurde sie bloß abgeschätzt. Sie 
betrug bedeutend weniger als 0,001 der intensiven Beleuchtung. 


7. Die Avogadro-Loschmidtsche ZahlN habe ich wegen 
le der Kontrolle der absoluten Größe der Verschiebungen aus 
1678 Punktierungen, welche auf den Mastix und Gummigutt- 
emulsionen bei sehr schwacher Beleuchtung gemacht wurden, 
berechnet und bekam für sie den Wert N = 6,22-1023, 


II. Besprechung der Resultate 


8. Über die Verläßlichkeit der Resultate. Die Tempe- 
ratur und Viskosität wurden so gemessen und berücksichtigt, 
daß in dieser Hinsicht kein Zweifel entstehen kann. Auch die 
Größe der Teilchen und die verschiedene Form der Rußteil- 
chen bleiben auf die Resultate ohne Einfluß, weil sich die auf 
sie bezogenen Größen bei der Bestimmung der Verhältnisse 
eliminieren. 

Die Hauptaufmerksamkeit müssen wir den möglichen 
Strömungen in der Flüssigkeit widmen, sie sind die wundeste 
Stelle der Arbeiten dieser Art. Daß die vorliegende Arbeit 
auch in dieser Hinsicht einwandfrei ist, geht aus folgendem 
hervor: 

a) Die Beobachtungsküvette wurde deshalb so klein ge- 
macht, damit die ganze Flüssigkeit gleichmäßig erwärmt wäre 
und keine Ursache zu Strömungen entstehen möchte. 

N b) Das gleichzeitige Suspendieren von Mastix und Kohlen- 
teilchen hatte eben den Zweck, zu beweisen, daß in der intensiv 
beleuchteten Flüssigkeit wirklich keine Strömung vorhanden 
war. Denn sonst müßten alle Substanzen dieselbe Vergröße- 
rung zeigen, wenn sie durch Strömungen verursacht sein sollte. 
Das Resultat bestätigt die obere Behauptung. 

E c) Daß die Bewegung der Kohlenteilchen auch während 
“) der intensiven Beleuchtung den Charakter einer reinen Brown- 
schen Bewegung trug, kann man folgendermaßen nachprüfen. 


Nach der Einsteinschen Formel sollen die 42, welche den 
Zeiten t, 2t, 31... entsprechen, im Verhältnis 1:2:3... stehen. 
Aus den Diagrammen kann man jetzt auch die Verschiebungen 
ay zwischen der 1. und 8., der 2. und 4. usw,, und zwischen der 


1. und 4., der 2. und 5. Position usw. ablesen und die 4? berechnen. 


0- 
e- 
u Ber. 
| 
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Ne 


So wurde gefunden, daß die 4?, welche den Zeiten 10, 20, 
30 Sek. entsprechen, im Verhältnis stehen: 


re Bor 1:2,02:2,81.... für die schwache Beleuchtung, und 
te Aue 1:1,95:2,86.... für die intensive Beleuchtung. 


Di BR, Dieses Resultat kann man mit Rücksicht auf die sta- 
tistische Methode und die fallende Zahl der zur Verfiigung 
stehenden Beobachtungen (166 für 20 Sek. und 151 für 80 Sek.) 
als gut bezeichnen. Die Abweichungen betragen in der Reihen- 
folge 1, 6, 2,5, 5 Proz. Wenn in der Flüssigkeit nur eine Spur 
von Strömungen vorhanden wäre, so wäre auch die Bahn der 
Teilchen teilweise gestreckt und die Quadrate müßten dann 
enorm steigen. 

Für die Vergrößerung der Bewegung der Kohlenteilchen 
für die Zeiten von 20 und 30 Sek. wurde nach der vorigen 
Methode gefunden, daß sie 


Ben ' ‘ fiir die Zeit von 20 Sek. 87 Proz., 


" Aus vorigen Ausführungen folgt, daß die vorliegenden 
| ey verläßlich sind und daß auch die vergrößerte Be- 
wegung der Kohlenteilchen rein brownisch war. 


ns, 9. Die Folgerungen aus den Messungen sind: 4 


1. Die Brownsche Bewegung existiert auch ohne Beleuch- 
tung, in der optischen Dunkelheit. Weil die Bewegung der 
_ Mastixteilchen unverändert bleibt, wenn die Beleuchtung 
Er. einer sehr niederen zu einer mehr als tausendmal höheren 
Intensität gesteigert wird, so wäre jeder Zweifel, daß die Be- 
wegung plötzlich aufhört, wenn die letzte optische Strah- 
lung verschwindet, unlogisch. 


2. Das absorbierte Licht vergrößert die Brownsche Bewegung. 
Damit soll die allgemein akzeptierte Gouysche Anschauung 
(Abs. 2) korrigiert werden. Der Irrtum von Gouy ist leicht 
erklärlich. Unter den beiden Kurventeilen A und B (Fig. 3) 
kann man wohl keinen auffallenden Unterschied wahrnehmen, 

ie und trotzdem ist das mittlere Quadrat der Verschiebungen 

in B um 40 Proz. größer als in A. Und wenn wir bedenken, 
daß die eingezeichneten geraden Kurvenstrecken in Wirk- 

Ye 4 lichkeit wieder verwickelte Brownsche Linien sind, so sieht 
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man, daß es unmöglich ist, die behandelte Vergrößerung un- 
zahlenmäßig zu konstatieren. 


10. Die Vergrößerung der Brownschen Bewegung 
durch das Licht läßt sich im Rahmen der kinetischen 
Theorie nicht auslegen. 

a) Die einfachste Erklärung dieses Phänomens ist in der 
Erwärmung der das beleuchtete Teilchen umgebenden Flüssig- 
keitsschichten zu suchen. Die von den Teilchen absorbierte 
strahlende Energie wird als Wärme den nächsten Flüssigkeits- — 
schichten abgegeben und durch ein stationäres Temperatur- 
gefälle wird sie in die ganze Flüssigkeit geleitet. Demnach 


> 


befindet sich das Teilchen ständig in einer wärmeren Flissig- = : 
keit als das Thermoelement anzeigt. Wir wollen jetzt prüfen, 
ob das Teilchen und die nächsten Flüssigkeitsschichten ge- — 
nügend hohe Temperaturen erreichen konnten, daß die ge- Br, ' | 
messene Bewegungsvergrößerung eintreten möchte. 
Aus der Einsteinschen Formel folgt 
A3. 43 T, „7, 
= 


Versuchsweise finden wir, daß die nächste Flüssigkeitsum- 

gebung bei den Kohlenteilchen um 9,3 0C über die Tempe- — 

ratur des Thermoelementes erwärmt sein müßte, wenn die 

45 proz. Vergrößerung stattfinden soll. Mi 

Aus der Wärmeleitungstheorie folgt, daß in einem sta- 

tionären Zustande eine Kugelschale von den Radien r, und 
r, (r,<r,) eine Wärmemenge Q 

Q=4 ah 

in der Zeiteinheit passieren läßt, wenn k den Warmeleitungs- — 


koeffizienten, und T,, T,(T,>T,) die Temperaturen der 
inneren und äußeren Kugelfläche bedeuten. Das Teilchen 


den Radius der als kugelförmig betrachteten Kohlenteilchen 
bedeutet. Für den stationären Zustand ist Q=Q’. Für 
den Temperaturunterschied AT = T,—T,, welcher zwischen _ 


de} 


| 
| 
| 
t absorbiert aus dem Lichtstrom I eine Wärmemenge Q’ er 
2 
Q’ = no?la, 
n wo a = 0,95 den Absorptionskoeffizienten und 9 = 4:10? cm Br; 
it 
| 


=. 


der innersten Flüssigkeitsschicht mit Radius =o und eine 
sehr entfernten Schicht r, (r > o) besteht, bekommen wir 


AT = 7% 9 = 0,0003°C 


Die nötige Temperaturdifferenz beträgt 9,39 C, die wirk- 
lich mögliche 0,0003° C. Wir sehen, daß die Erklärung unserer 
Erscheinung durch die Erwärmung der Teilchenumgebung 
nicht gelingt. 

b) Wenn wir außer der obigen Erwärmung der Flüssig- 
keit noch einen möglichen Temperatursprung an der Grenz- 
fläche Teilchen-Flüssigkeit voraussetzen und so die Erklärung 
nur durch die erhöhte kinetische Energie des Teilchens suchen, 
so finden wir, daß dieser Sprung bei den Kohlenteilchen 132°C 
betragen müßte, wenn die gegebene Vergrößerung stattfinden 
sollte.2) Die Temperatur der Kohlenteilchen müßte dann 
einen unzuläßlichen Wert von 152° C erreichen. 

c) Die behandelte Erscheinung erinnert an den Crookes- 
schen Radiometer. Die unter dem Einfluß des Lichtes leb- 
hafter sich bewegenden Rußteilchen sind ein mikroskopischer 
Flüssigkeitsradiometer. Wenn wir dieselbe auf die übliche 
Weise durch den Rückstoß der erwärmten Flüssigkeitsmole- 
küle zu erklären versuchen, so stoßen wir auf folgendes Di- 
lemma. Entweder müssen die erwärmten Moleküle so hohe 
Temperatur besitzen, daß die gefundene Vergrößerung ein- 
treten kann, und wir finden wie im Falle sub b), daß sie bis 
4 auf 152°C erwärmt sein müßten, was der Wärmeleitungs- 
theorie widerspricht, oder wenn wir diese berücksichtigen, 
bekommen wir, daß die Moleküle wie in dem Falle sub a) nur 
so eine Erhöhung ihrer kinetischen Energie erfahren könnten 
welche einer Temperaturerhöhung von nur 0,0003°C ent- 
spricht und die Vergrößerung kann dann nicht eintreten. 


11. Versuch einer Erklärung der Erscheinung. 

a) Kann die Bewegungsvergrößerung photoelektrisch er- 
klärt werden? Setzen wir voraus, daß das von den Kohlen- 
teilchen absorbierte Licht teilweise zur Emission der Photo- 


1) Das Leitungsvermögen k des benutzten Gemisches Wasser- 
Glyzerin wurde zu 1,2-10-* berechnet. 
4 BR 2) Berechnet aus dem Verhältnis 
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elektronen dient. Die Rückstöße, welche diese Emission be- 
gleiten, möchten dann die Brownsche Bewegung der Teil- 
‘chen vergrößern. Diese Idee wollen wir jetzt auf unserem 
konkreten Falle durch folgende Umschlagrechnung prüfen. 


Das Kohlenteilchen machte in einer Sekunde die mittlere 


räumliche Verschiebung 4 1:10-* cm. Betrachten wir sie 


als eine geradlinige. Bei der intensiven Beleuchtung absor- — Br 


bierte das Teilchen die Energie W = 7,3-10-%erg/sec und — 


seine räumliche Verschiebung betrug 1,21-10- cm/sec. Die 
Bewegungsvergrößerung d 4 = 2,1-10-5 cm/sec. Dieser Zu- 


wachs dA stammt nicht aus der Molekularbewegung der 


Flüssigkeit und das Teilchen verbraucht nach dem Stokes- 


schen Gesetze die Energie 
Wy = 6%le(d 4)? =8,9-10-erg/seo 


um ihn zu vollziehen. Zur Berechnung der Bewegungsenergie — 
W,, welche das Teilchen durch die vorausgesetzte Elektronen- © 


emission bekommen kann, benutzen wir folgenden Satz, der BER 


unmittelbar aus der Impulsgleichheit folgt. 


Aus der Energie e eines Emissionsvorganges bekommt 


der Körper m, die Energie e, und der Körper m, die Energie e, 
so, daß das Produkt Energie x Masse konstant bleibt. 
em, = +&%=e. 

Weil die Masse des Kohlenteilchens m, = 4-10—1% g un- 
gemein ete | als die des Elektrons m, = 9-10-%8 g ist und 
setzen wir e = W und e, = W,, so bekommen wir einfach 
W, = Wm,/m, = 1,6-10—" erg/sec. Weil W,> W,, sehen wir, 
daß auch in dem Falle, wenn die gesamte absorbierte Energie 
photoelektrisch verbraucht und dabei alle Elektronen in der- 
selben Richtung emittiert würden, die Rückstoßkräfte nicht 


genügen die beobachtete Vergrößerung zu verursachen. — 


Man könnte trotzdem auf der photoelektrischen Erklärung 
unter der Annahme verharren, daß das emittierte Elektron 
bei der Zurückkehrung dem ionisierten Brownschen Par- 
tikel einen kleinen Teil seiner Bewegungsenergie zurückgibt.!) 


1) Z. B. wenn ein Elektron, welches mit dem Quant hy = 5-10-1 erg 
(berechnet für violettes Licht) emittiert wurde, nur !/, Proz. seiner Be- 
wegungsenergie dem Teilchen zurückgibt, würde schon die nötige Energie 
gedeckt sein. 
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b) Die molekulare und die Brownsche Bewegung und 
thre Vergrößerung durch das Lächt ist eine Folge der Strahlung, 


Betrachten wir die molekulare Bewegung als keine hypo-' 
 thetische, sondern als eine experimentelle Tatsache. 


Wenn 
wir die mechanische Weltanschauung verlassen, so handelt 
es sich darum, wie die Bewegung von einem energetischen 


Standpunkte aus zu erklären wäre. 


Jedes Massenteilchen, dessen Temperatur oberhalb des 
absoluten Nullpunktes liegt, steht durch unaufhörendes Strahlen 
mit der Umgebung im dynamischen Temperaturgleichgewichte, 


. Von dieser Strahlung setzen wir voraus, daß sie sich nicht 


strahlt (emittiert oder absorbiert) wird. 


Tu 


> 


ix 


fl 
Fin 


nur quantenhaft, sondern auch strahlenartig abspielt, das 
ist, daß ein Energiequantum & in einer festen Richtung ge- 
Jeder dieser Akte 
wird mit einem Impuls ¢/f(c) begleitet. f(c) bedeutet eine 
Funktion der Lichtgesch windigkeit.!) 

Da die Strahlung vollkommen unregelmäßig, sagen wir 
„brownisch‘‘ vor sich geht, so wird das Massenteilchen durch 
die Impulse hin und her bewegt. Auch ein einziges in einem 
Vakuum schwebendes Massenteilchen macht die Brownsche 
Bewegung. Eine geradlinige Molekularbewegung existiert 
praktisch nicht Die Umwandlung eines Bruchteiles der 
strahlenden Energie in die Bewegungsenergie des Teilchens 
und umgekehrt äußert sich durch den Dopplerschen Effekt. 
Wenn der Beweis gelingen wird, daß auch bei dieser An- 
schauungsweise die mittlere kinetische Energie eines Massen- 


 teilchens gleich 3RT/2N ist, so können wir uns auch von 


diesem Standpunkte aus die schönen Resultate der kinetischen 
Gastheorie aneignen. 

Die zwei konstatierten Tatsachen (Abs. 9) sind jetzt 
leicht qualitativ erklärbar. Die Brownsche Bewegung dauert 
auch in der optischen Dunkelheit, weil bei der gewöhnlichen 
Temperatur die Moleküle und die Teilchen im langwelligen, 
unsichtbaren Gebiete strahlen. Absorbiert das Teilchen aus 
dem optischen Energiestrom die Energie W, so wird seine 
Brownsche Bewegung durch diese Absorption, weil sie nur 
in einer Richtung geschieht, nicht beeinflußt, sondern erst 
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durch die langwellige, brownische Wiederausstrahlung der- _ 
selben wird sie gesteigert. = 


Im Mangel einer strengen Theorie versuchte ich durch — 
folgende Überlegung auch zu einem quantitativen Aufschlusse = 
darüber zu gelangen. Das fast vollkommen schwarze Kohlen- 
teilchen tauscht mit der Umgebung in einer Sekunde de 
Energie W* — 8x 020 T? aus, wo 420? die Teilchenoberfliche 
und o die Stefan-Boltzmannsche Konstante bedeutet. 
Während der intensiven Beleuchtung ist seine gesamte brownisch _ 
gestrahlte Energie gleich W* + W. j 


Die Energie W* bestehe aus n-Quanten, von deren Größe © 
wir nichts voraussetzen. Legen wir durch das Teilchen eine 
feste Achse x. Die Wahrscheinlichkeit, daß ein Quantum 
so strahlt, daß seine z-Komponente des Strahlungsdruckes — 
in die +a oder —z-Richtung fällt, ist gleich groß und ist Y,. 
Für diesen Fall folgt aus der Wahrscheinlichkeitsrechnung, 
daß aus den n-Quanten ein Überschuß q 


a 


+ 
in einen oder anderen Richtung (+2 oder zu er- 
warten ist. Dieser Überschuß q; maitiplisiert mit den zu- = 


die wahrscheinliche resultierende Kraft, welche das 
verschiebt. 
Die Verschiebung 4 ist also proportionell der Quadrat- — 


wurzel aus der gestrahlten Energie und die Quadrate 4? sind 
proportionell der gestrahlten Energiemenge W*, bzw. W* +W. 
So folgt für das Verhältnis 


4 A,?:4;2 = W*:(W* + W) = 80 T!:(86 74 + IA). 
Wenn wir in die letzte Gleichung die Zahlenwerte (6 = 
5,65-10-5 erg/sec- qem grad*, T = 293, I = 2,2/60 cal/sec- gem, 

A = 4,19-107 erg/cal) einsetzen, so berechnen wir für das 


Verhältnis der Verschiebungen der Kohlenteilchen 100:146, 
und das Resultat der Messungen war 100:145. : 


Von dieser interessanten Ubereinstimmung kann man 
wohl nicht behaupten, daß sie einen sicheren Beweis für aa 
oben gemachten Voraussetzungen darstellt, aber trotzdem bi 
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8 W. Pospisil. Vergrößerung der Brownschen Bewegung usw. 


kann sie die Vermutung berechtigen, daß in der Sache ein 
‚großes Stück der objektiven Wahrheit liegt.!) 

Ich hoffe, daß sich ein Theoretiker findet, der die hier 
skizzierte Frage lösen wird. 

Ich halte es für meine angenehme Pflicht, an dieser Stelle 
dem Hrn. Prof. Dr. V. Posejpal, welcher mir die Durch- 
führung dieser Arbeit ermöglichte, meinen besten Dank aus- 


zusprechen. 4 
Zusammenfassung a 


Es wurde experimentell quantitativ festgestellt 
a) Im Einklang mit der bisherigen Anschauung, daß 
die Brownsche Bewegung der Substanzen, welche das be- 
nutzte Licht nicht absorbieren (Mastix), ohne Rücksicht auf 
die Intensität der Beleuchtung unverändert bleibt und daß 
sie deshalb auch in der optischen Dunkelheit andauert. 
¢ b) Im Gegensatz zu der bisherigen Anschauung, nach 
welcher die benutzte intensive Beleuchtung die Bewegung 
nur in dem Maße vergrößern kann, soweit sie die Temperatur 
der Flüssigkeit erhöht und ihre Viskosität herabsetzt, wurde 
konstatiert, daß die Bewegung der das Licht absorbierenden 
Substanzen (Kohle, Gummigutt) durch eine intensive Be- 
 leuchtung um viel mehr vergrößert wird, als zu erwarten wäre. 
’ Die Bewegungsvergrößerung stellt eine verwandte Er- 
_ scheinung der Photophorese und des Radiometers vor und läßt 
sich nicht im Rahmen der kinetischen Theorie erklären. Sie 
könnte entweder photoelektrisch oder besser unter der An- 
nahme, daß die molekulare Bewegung eine Folge der schwarzen 
Strahlung ist, erklärt werden. 


> Gearbeitet 1921—28 im Physikalischen Institute der 
Karlova-Universität in Prag. 


1) Der Gedanke, daß der Strahlungsdruck die molekulare Bewegung 
verursachen kann, ist nicht neu, wie ich aus dem Buche Eucken, „Die 
Theorie der Strahlung und der Quanten“, erfahren hatte. Besonders 
interessant ist die Vermutung von Solvay (dortselbst S.3 und 369), 
mach der sich ein Teilchen nur dann bewegt, wenn es strahlt. Vgl. auch 
Einstein und Hopf, Ann. d. Phys. 88. S. 1105. 1910. Mit Anwendung 
des Strahlungsdruckes könnte man auch die Ehrenhaftsche Photo- 
_ phorese unter der Annahme erklären, daß die intensiv beleuchteten sub- 
re mikroskopischen Teilchen die absorbierte Energie gréBtenteils entweder 
_ auf der beleuchteten oder auf der unbeleuchteten Seite ausstrahlen. 
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